Robotics  Research 
l^hnical  Report 


II 


Polygon  Optimization  Problems 

by 

Jyun-Sheng  Chang 


Technical  Report  No.  240 

Robotics  Report  No.  78 

Aug,  1986 


\, 


c 

^    (Jl-U 

f^   C   (0 

I     <U    N 

«x:  -H 
Eh  w  e 

I     -H 
M    C  4J 
CJ   3    a    . 
W  >i  o  w 
Oj  t-a       g 

S         CO) 
O      -  O  i-H 

C  >i  o 

>;  x:  O  Q< 
z  u  cu 


New  York  University 
nt  Institute  of  Mathematical  Sciences 

Computer  Science  Division 

25 1  Mercer  Street   New  York,  N. Y  1 00 1 2 


Polygon  Optimization  Problems 

by 

Jyun-Sheng  Chang 


Technical  Report  No.  240 

Robotics  Report  No.  78 

Aug,  1986 


New  York  University 

Dept.  of  Computer  Science 

Courant  Institute  of  Mathematical  Sciences 

251  Mercer  Street 

New  York,  New  York    10012 


Work  on  this  paper  has  been  supported  by  Office  of  Naval  Research  Grant  N00014-82-K- 
0381,  National  Science  Foundation  CER  Grant  DCR-83-20085,  and  by  grants  from  the 
Digital  Equipment  Corporation  and  the  IBM  Corporation. 


POLYGON  OPTIMIZATION  PROBLEMS 

by 

Jyun-Sheng  Chang 

New  York  University 

Advisor:  Chee  K.  Yap 

ABSTRACT 

This  thesis  examines  polygon  optimization  problems  arising  from  the  stock- 
cutting  problem.  Two  types  of  problems  are  considered:  the  inclusion 
problems  and  the  enclosure  problems.  The  inclusion  (enclosure)  problems 
ask  for  a  maximum  polygonal  subset  (minimum  polygonal  superset)  of  a 
given  polygon,  satisfying  certain  conditions.  Both  the  area  and  perimeter 
metrics  on  the  polygons  can  be  used  as  the  measure  of  optimality. 

Various  geometric  properties  and  algorithms  for  these  problems  are  shown. 
The  mam  results  are:  (1)  An  O(n^)  time  iO(n^)  time)  algorithm  for  fmding 
a  maximum  area  (perimeter)  convex  subset.  (Only  exponential  time 
algorithms  existed  previously  for  the  problem.)  (2)  An  0(n'  log  n  log  k) 
time  algorithm  for  finding  a  minimum  area  enclosing  convex  k-gon.  (3)  An 
0{n~)  time  algorithm  for  finding  a  minimum  perimeter  enclosing  triangle. 
(4)  An  0(nk'^)  time  algorithm  for  finding  a  minimum  enclosing  k-gon  with  a 
fixed  shape. 
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CHAPTER   1 


Introduction 


1.1.   Computational  Geometry 

Geometric  objects  such  as  points,  liaes  and  polygons  are  the  basis  of  a  broad  variety  of 
important  applications  and  give  rise  to  an  interesting  set  of  problems  and  algorithms.  The 
name  geometry  reminds  us  of  its  earliest  use  -  for  the  measurement  of  land  and  materials.' 
Today,  computers  arc  being  used  more  and  more  to  solve  larger-scale  geometric  problems. 
For  example,  VLSI  circuit  design  requires  a  test  for  intersection  of  line  segments;  computer 
graphics  involves  hidden  line  removal.  Over  the  past  two  decades,  a  set  of  tools  and 
techniques  has  been  developed  that  takes  advantage  of  the  structure  provided  by  geometry. 
This  discipline  is  known  as  Computational  Geometry . 

Early  works  in  computational  geometry  have  followed  four  major  directions:  searching 
problems,  the  construction  of  convex  hulls,  proximity  problems,  and  intersection  problems 
[Shamos-Preparata  (1985), Brown  (1979),Dobkin-Snyder  (1979)]. 

In  the  following,  we  give  a  brief  review  of  these  problems.  For  a  more  comprehensive 
survey,  see  [Brown  (1979)]. 

1.1.1.    Searching  Problems 

Searching  is  a  fundamental  operation  intrinsic  to  many  computational  tasks;  it  is  very 
common  to  retrieve  some  particular  information  (to  query)  from  a  large  amount  of  previously 
stored  information  [the  file)     There  are  rwo  kinds  of  geometric  searching  problems: 

(1)  Location  problems,  where  the  file  is  a  set  of  regions  defined  by  some  geometric  objects 
such  as  line  segments  and  the  query  is  the  region  that  contains  a  given  point. 

(2)  Range  searching  problems,  where  the  file  is  a  set  of  geometric  objects  such  as  points  and 
the  query  is  the  subset  of  the  file  that  is  contained  in  a  given  geometric  shape  in  the 


'  Concrete  evidents  [Eves  (1972)]  exist,  indicating  that  the  Babylonians  and  Egyptians  had 
accumulated  a  lot  of  empirical  knowledge  of  geometric  shapes.  See  [Eves  (1972)]  for  a 
complete  survey  of  geometry.  See  also  [Shamos  (1978)]  for  the  history  of  geometry  from  the 
perspective  of  a  computer  scientist. 


space. 

Sec  [Shamos  (1978),Dobkin-Lipton  (1976),Lce-Preparata  (1977),  Lipton-Tarjan 
{1980),BcDtlcy-Saxc  (1979)]. 

1.1.2.  Construction  of  Convex  Hulls 

The  convex  hull  of  a  set  of  points  is  not  only  a  fundamental  geometric  structure,  but  is 
also  a  vehicle  for  the  solution  of  many  problems  in  the  literature  and  this  thesis. 

The  computation  of  convex  hull  has  been  studied  extensively.  See  [Gramham 
(1972),Jarvis  (1973),Prcparata-Hong  (1977).Bcntlcy-Shamos  (1976)]. 

1.1.3.  Intersection  Problems 

Intersection  problems  arise  in  such  areas  as  computer  graphics  (hidden  line  removal) 
and  printed  circuit  design.  Sec  [Shamos-Hocy  (1976),  Bentley-Ottmann  (1979),  Preparata- 
MuUer  (1979)]. 

1.1.4.  Proximity  Problems 

Proximity  problems  include  the  following  apparently  diverse  problems: 

(1)  The  closest  pair  of  a  planar  point  set. 

(2)  The  nearest  neighbor  of  each  point  in  a  planar  point  set. 

(3)  The  spanning  tree  of  a  planar  point  set  with  minimum  total  length 

(4)  The  nearest  neighbor  of  a  given  point  (in  a  point  set  given  beforehand). 

See  [Bentley-Shamos  (1976),  Shamos-Hoey  (1975),  Yao  (1975),Cheriton-Tarjan  (1976)]. 

1.2.    Optimization  Problems 

Our  goal  in  this  thesis  is  to  explore  a  fifth  direction  in  computational  geometry,  namely 
Optimization  Problems 

Optimization  problems  involve  finding  an  object  that  r^ptimizes  certain  objectives  while 
satisfying  various  constraints.  In  the  case  of  geometric  optimization  problems,  we  are  given 
a  set  of  geometric  constraints  and  arc  to  find  some  geometric  object  with  minima!  or  maximal 
geometric  properties  such  as  distance,  perimeter,  or  area 

In  this  general  sense,  proximity  problems  can  also  be  viewed  as  optimization  problems. 
However,  using  the  locus  approach ,'  the  loci  of  proximity  (the  Voronoi  digram)  arc  found; 
the  plane  is  partitioned  into  regions,  each  of  which  is  the  locus  of  the  points  closer  to  some 
point  than  to  any  other  point  in  the  file.  By  locating  the  query  point  in  a  region  in  this 
partition,  the  closest  point  is  found      In  other  words,  the  proximity  problems  are  solved  by 


"  There  is  also  a  dividc-and-conquer  approach  to  the  closest  pair  problem  [Bentley-Shamos 
(1976),Bcntley  (1980)]. 


searching  the  Voronoi  digram  of  the  file.    There  are  other  optimization  problems  that  do  not 
have  such  proximity  property  and  are  not  easily  solved  by  searching. 

In  this  thesis,  we  focus  on  optimization  problems  in  which  the  constraints  are  expressed 
as  containment  between  polygons.  The  general  framework  for  this  kind  of  optimization 
problems  can  be  posed  as  follows:  Let  P,  Q  be  families  of  polygons,  and  let  )x  be  a  real 
function  on  polygons  with  the  property  that  for  all  P  in  P,  Q  in  Q: 

P  is  contained  m  Q     =>     ^(/')  s  ^(C). 

Note  that  this  property  holds  if  p.  measures  the  area.  If  the  polygons  in  Q  are  convex, 
the  property  also  holds  when  ji  measures  the  perimeter.'  With  this  notation  the  classes  of 
(polygon)  inclusion  and  enclosure  problems*  can  now  be  defined: 

Inc(^P,  Q,  \i)  :  Given  P  i  P,  find  the  ji -largest  Q  ^  Q.  that  is  included  in  P. 

Enc{P,  Q,  ji)  :  Given  Pip,  find  the  ^.-smallest  Q  d  Q  that  encloses  P. 

Typical  examples  for  Pare  P^,„,  the  family  of  all  simple  polygons  and  f^^„,  the  family  of 
all  convex  polygons.  Q  can  typically  be  P^,  the  family  of  all  triangles  or  more  generally  P^, 
the  family  of  all  convex  Jt-gons. 

These  problems  are  posed  in  a  continuous  (as  opposed  to  a  discrete  or  combinatorial) 
setting.  Since  we  want  to  design  discrete  algorithms  to  run  on  digital  computers,  it  is 
essential  that  we  first  find  a  "finiteness  criterion"  for  the  solution  of  each  problem.  Once  the 
finitencss  criterion  is  found,  we  can  combine  it  with  various  problem  solving  paradigms  to 
develop  an  efficient  algorithm  for  the  problem  The  themes  of  this  thesis  will  be  the  study  of 
the  geometric  properties  of  the  polygon  optimization  problems,  finding  fmiteness  criteria  for 
these  problems,  determining  their  usefulness  for  an  algorithm,  and  the  actual  constructions 
and  analyses  of  algorithms. 

1.3.    Previous  Work  in  Geometric  Optimization  Problems 

We  review  some  of  the  inclusion  and  enclosure  problems  that  have  been  studied.  The 
enclosure  problem  is  closely  related  to  the  "stock-cutting  problems"  which  are  concerned  with 
cutting  a  sheet  of  material  into  smaller  subparts  under  various  constraints  (such  as  all 
subparts  are  congruent  to  a  given  shape)  and  satisfying  some  optimality  criterion.  See  [Dori- 
Ben-Bassat  (1983),  Albano-Orsini  (1979),  Christofidcs-Whitlock  (1977),Gilmore-Gomory 
(1967)].  As  pointed  out  in  [Dori-Ben-Bassat  (1983)],  the  enclosure  problem  is  a  key 
subproblem  in  the  more  general  stock-cuning  problems.    More  generally,  the  enclosure  and 


This  can  be  shown  quite  easily  by  an  induction  argument  and  the  triangle  inequality. 
Alternatively,  these  might  be  called  the  inscription  and  the  circumscription  problems. 


inclusion  problems  can  be  viewed  as  a  way  of  approximating  polygons.  For  example,  for  the 
purposes  of  detecting  collision  in  robotics,  we  typically  approximate  a  complicated  shape  by  a 
simpler  enclosing  body. 

1.3.1.    Inclusion  Problems 

The  inclusion  problem  Inc{^^^„,  Py,  area)  was  first  considered  by  Dobkin  and  Snyder. 
Their  algorithm  runs  in  linear  time  [Dobkin-Snyder  (1979)].  The  techniques  were  extended 
by  Boycc,  Dobkin,  Drysdale,  and  Guibas  to  the  problems  /nc(.P„„,  P^,  area)  and 
lnc{P^^„,  P^,  perimeter)  for  any  fixed  *  [Boyce-Dobkin-Drysdale-Guibas  (1982}].  Their 
algorithms  run  in  0(k  n  log-n)  time  Note  that  unlike  the  more  general  problems  described 
below,  the  finiteness  of  these  problems  is  easy  to  show.  It  follows  from  the  fact  that  the 
vertices  of  any  maximal  *-gon  must  be  a  subset  of  the  vertices  of  the  input  polygon.  So  the 
techniques  of  [Dobkin-Snyder  (1979),Boyce-Dobkin-Drysdalc-Guibas  (1982)]  are  not 
sufficient  for  the  more  general  problems  such  as  /nc(/'j^,  f^,  area),  since  they  rely  on  the 
convexity  of  the  input  polygon. 

The  potato-peeling  problem  was  first  posed  (in  a  more  general  form)  by  Goodman  who 
obtained  various  mathematical  properties  of  solutions  to  the  problem.  The  potato-peeling 
problem  in  our  notation  is  /nf(,r,„,  P^^„,  area).  The  only  partial  solution  in  [Goodman 
(1981)]  gives  the  criterion  for  a  finite  solution,  if  the  polygon  has  n  s  5  sides.  Quite 
independently,  Woo  studied  the  same  problem,  dubbing  it  the  "convex  skull"  problem  [Woo 
(1983)].  The  area  and  perimeter  versions  of  the  potato-peeling  problem  are  solved  in  0{n') 
and  O(n')  time  respectively  in  [Chang-Yap  (1984), Chang-Yap  (1986)].  Rectilinear  versions 
of  the  potato-peeling  problem  have  been  addressed  in  [Chazclle-Drysdalc-Lee 
(1984),McKenna-0'Rourke-Suri  (1985),  Wood-Yap  (1985)].' 

1.3.2.    Enclosure  Problems 

The  simplest  case  of  the  enclosure  problem.  Enc{P,,^.  P^,  area)  was  solved  by  Klce  and 
Laskowski  in  0{n  log-n)  time  [Klee-Laskowski  (1985)]  It  was  improved  by  O'Rourke, 
Aggarwal,  Maddila,  and  Baldwin  to  linear  time  which  is  optimal  [U'Rourke-Aggarwal- 
Maddila-Baldwin  (1984)].  For  the  problem  £nr(",  ,...  ^,.  area),  [Chang-Yap  (1984)]  provided 
an  0(n'log  *)  time  algorithm.  By  further  refincmeni.  [Aggarwal-Chang-Yap  (1985)]  solved 
the  same  problem  in  C>(n-log  n  log  k)  time.  Wc  remark  that  Dori  and  Ben-Bassat  claimed  to 
have  a  linear  time  solution  to  this  problem  [DoriBen-Bassat  (1983)].  However,  their 
optimality  proof  is  faulty;  indeed  [O'Rourke  (1984b)]  has  provided  some  counterexamples. 

Other  researchers  have  restricted  the  solution  to  sub-families  of  convex  polygons  with 
the  intent  of  obtaining  faster  algorithms.  The  restricium  may  be  imposed  on  the  angles  of 
the  polygon,  on  the  lengths  of  its  edges,  or  on  both 


^  [Cha7clle-Drysda!e-Lce  (1984)]  formulates  the  problem  slightly  differently  but  [McKenna- 
O'Rourke-Suri  (1985)]  shows  the  connection  to  pot.no  peeling. 


As  an  example  of  angular  restriction,  we  can  restrict  the  family  f^  to  rectangles. 
Freeman  and  Shapira  proposed  an  O(n^)  time  algorithm  [Frccman-Shapira  (1975)]  for  this 
problem.    Toussaint  improved  it  to  linear  time  [Toussaint  (1979)]. 

In   general,   let  0  =  (B, 9^)   be   any   sequence   of   angles   with   each   6,  <  ti    and 

k 

21  6,  =  (k-2)-n.  Let  Pq  denote  the  family  of  convex  t-gons  whose  interior  angles  are  given 
/- 1 

by  the  sequence  6.  Thus  we  have  the  problem  of  finding  smallest  polygon  from  Pq  enclosing 
a  given  convex  polygon:  Enc{P^^^,  Pq,  ^).  Depano  and  Aggarwal  solved  the  problem  in 
0{kn)  time  [DcPano-Aggarwal  (1984)].  For  the  three  dimensional  versions  of  these 
problems,  O'Rourke  described  an  O(n^)  time  algorithm  for  the  smallest  rectangular  box 
enclosing  a  polyhedron  [O'Rourke  (1984a)]. 

We  can  also  impose  a  restriction  that  involves  both  angles  and  the  lengths  of  the  edges: 
restricting  the  solution  polygon  to  regular  polygons.  An  0{nT{k))  time  solution  was  given  in 
[DePano-Aggarwal  (1984)]  for  finding  smallest  regular  t-gons  enclosing  a  given  polygon 
where  T{k)  is  an  exponential  function  of  k. 

In  general,  let  P  be  any  polygon  and  shape{P)  be  the  family  of  polygons  obtained  by 
the  transformations  of  scaling,  rotation,  translation  and  reflection  of  P .  Then  we  have  the 
shape  problems,  lnc{P^^„,  shape{P),  p.)  and  Enc{P^,^,  shape{P),  \i.).  We  will  describe  an 
0{nk^)  time  solution  for  the  enclosure  shape  problem  (*  is  the  number  of  edges  in  P). 
Obviously,  for  this  particular  problem,  whether  p.  measures  area  or  perimeter  makes  no 
difference. 

A  special  situation  in  the  general  stock-cutting  problem  is  the  construction  of  the 
congruent  replications  of  a  given  shape  in  a  board  so  that  waste  is  minimized.  This  problem 
has  a  practical  application  in  mass-production  systems  where  the  board  area  is  relatively  large 
compared  to  the  given  figure. 

One  approach  [Dori-Ben-Bassat  (1984)]  often  taken  to  solve  the  problem,  is  the 
following:  First  enclose  the  given  figure  with  a  smallest  possible  polygon  that  can  pave  the 
plane,  or  in  other  word,  cover  the  plane  by  replications  of  the  same  figure  without  any  gap  or 
overlap.  One  kind  of  such  polygons^  [Kershner  [1968]]  is  the  family,  Wpo,,,,  of  hexagons. 
For  each  hexagon  P  in  WL^,.  there  exists  a  pair  of  edges  e  and  e   \t\  P  such  that 

(1)  e  and  e   are  two  edges  apart, 

(2)  e  is  parallel  to  e  ,  and 

(3)  1^1  =   \e  \\  e  and  e   equal  in  length. 

Then  the  problem  in  question  reduces  to  £nr(P,,^,  //:    ^.,  area)  -   the  problem  of  finding  the 
smallest  hexagon  in  WL,,,  enclosing  a  given  polygon. 


1.3.3.    Smallest  Enclosing  Circles  and  Ellipses 

Optimization  problems  involving  geometric  shapes  other  than  polygons  are  also  of 
interest.    Minimum  area  enclosure  with  circles  and  ellipses  has  been  studied. 

As  pointed  out  by  Post,  an  ellipse  enclosing  a  convex  polygon  provides  an  analytic 
approximation  for  the  polygon  which  may  be  useful  in  certain  situations  [Post  (1981)].  For 
instance,  approximation  algorithms  for  point  location  can  be  developed  which  run  in  constant 
time.  Consider  the  problem  of  determining  whether  a  given  point  is  contained  in  the 
polygon.  There  is  an  exact  algorithm  [Shamos  (1978)]  that  runs  in  0{n  log  n)  time. 
However,  determining  whether  a  point  is  in  an  ellipse  is  just  a  matter  of  substituting  the 
coordinates  of  the  point  into  the  formula  that  determines  the  ellipse.  If  a  large  set  of  points 
needs  to  be  checked  for  inclusion  in  a  fixed  polygon,  we  can  use  the  minimum  area  enclosing 
ellipse  to  quickly  throw  out  all  the  points  that  are  not  in  the  ellipse  and  check  the  points  that 
are  in  the  ellipse  using  the  exact  algorithm.  Post  described  an  0(n'^)  time  algorithm  for 
determining  the  minimum  area  ellipse  enclosing  an  n-sided  polygon  [Post  (1981)].* 

Finding  the  smallest  circle  that  encloses  a  given  set  of  planar  points  is  discussed  in 
[Rademacher-Toeplitz  (1957)]  It  also  appears  in  operations  research  as  the  minimax 
facilities  location  problem  in  which  we  seek  a  point  whose  greatest  distance  to  any  point  of 
the  set  is  a  minimum.  Elzinga  and  Hearn  gave  a  quadratic  algorithm  for  the  problem 
[Elzinga-Hearn  (1972)).  Mcgiddo  improved  the  running  time  to  0(n)  which  is  optimal  for 
the  problem  [.Megiddo  [1982]] 

In  this  thesis,  we  will  concentrate  on  polygonal  optimization  problems. 

1.4.   Thesis  Outline 

Chapter  2  establishes  the  context  and  direction  of  this  thesis  with  a  simple  example  - 
the  minimum  area  corner  containing  a  set  of  given  points.  But  before  the  problem  is  looked 
at  and  an  algorithm  given,  we  first  settle  such  basic  issues  as  problem  specification, 
representation  of  geometric  objects,  and  computational  model.  We  show  the  process  of 
finding  a  finiteness  criterion  for  the  solution  and  the  process  of  developing  an  algorithm 
based  on  the  finiteness  criterion. 

Chapter  3  is  a  stud\  of  minimum  area  enclosure.  We  look  at  the  general  case  where  the 
enclosing  polygons  have  k  sides,  k  2  3.  We  present  a  finiteness  criterion  and  develop  an 
algorithm  for  the  minimum  area  enclosing  k-gon  The  algorithm  employs  the  techniques  of 
divide-and-conquer  and  dynamic  programming  and  trades  storage  for  fast  running  time.  The 
analysis  of  the  running  time  and  storage  requirement  of  the  algorithm  is  also  given  This  and 
all  succeeding  chapters  conclude  with  a  list  of  unsolved  problems. 

Chapter  4  is  concerned  with  finding  minimum  perimeter  enclosing  triangles  V*  c 
present  the  necessary  finiteness  criterion  and  an  algorithm  that  applies  the  dividc-and-conquer 


The  running  time  has  been  improved  to  0(n-]  recently  by  the  same  author. 


strategy  in  aa  interesting  way. 

In  Chapter  5,  we  introduce  the  basic  geometric  properties  of  the  inclusion  problem  and 
solve  a  geometric  problem  that  arises  from  the  inclusion  problem.  Finiteness  criteria  and 
algorithms  for  the  inclusion  problem  (both  area  and  perimeter  metrics)  arc  presented  and 
analyzed  in  Chapter  6  and  7. 

In  Chapter  8,  wc  cover  several  miscellaneous  problems  that  are  the  restricted  versions 
of  the  problems  discussed  in  the  previous  chapters. 

Chapter  9  summarizes  the  thesis  and  points  out  directions  for  further  work. 


CHAPTER  2 


An  Example:  The  Smallest  Enclosing  Corner  Triangle 


In  this  chapter,  wc  explore  some  of  the  basic  issues  that  arise  as  we  attempt  to  solve 
geometric  problems  on  a  computer.  To  illustrate  these  issues,  Wc  will  use  the  simple 
example  of  finding  a  smallest  triangle  at  the  corner  of  the  first  quadrant,  enclosing  a  set  of 
given  points.    In  solving  the  problem,  we  face  the  following  questions: 

(1)  A    precise    problem    specification;    how    the    input    and    output    of    the    problem    are 
represented. 

(2)  A  model  of  computation  and  an  appropriate  cost  measure  to  measure  the  complexity. 

(3)  Upper  bounds  and  lower  bounds  of  the  problem. 

2.1.    The  Model  of  Computation 

Wc  will  use  a  mixture  of  English-like  statements,  basic  control  structures  (sequence, 
branch,  repetition,  and  recursion),  and  traditional  mathematical  notation  to  present  our 
solution  to  a  problem.  However,  to  be  able  to  analyze  the  computational  complexity  of  an 
algorithm,  wc  need  to  relate  the  notation  we  use  to  a  formal  computational  model. 

A  model  of  computation  defines  the  basic  operations  that  algorithms  can  perform  and 
the  cost  of  each  operation.  It  is  essential  to  have  a  precisely  specified  model  of  computation 
in  order  to  talk  about  such  things  as  how  much  time  is  sufficient  or  necessary  to  solve  a 
certain  problem  -  upper  bounds  and  lower  bounds.  The  choice  of  a  computation  model  is  a 
compromise  between  realism  and  mathematical  tractability .  It  is  to  be  as  close  to  actual 
computers  as  possible  and  yet  still  permit  thorough  analysis 

The  choice  depends  on  the  kind  of  problems  one  is  dealing  with  For  instance,  in 
searching  problems,  we  are  given  a  collection  of  objects  and  are  to  choose  an  element  (or  a 
subset)  from  the  collection  satisfying  some  given  criterion.  The  solution  of  these  problems 
consists  of  objects  given  and  wc  do  not  have  to  create  any  new  objects.  Often  such  problems 
can  be  solved  by  using  only  comparisons  and  the  four  arithmetic  functions.  On  the  other 
hand,  with  problems  such  as  finding  the  smallest  triangle  containing  a  given  planar  point  set, 
we  need  to  compute  the  geometric  parameters  of  the  triangle  For  these  problems,  we  often 
require  more  general  functions. 


The  real  RAM  seems  to  be  a  reasonable  choice  for  much  of  the  work  in  computational 
geometry.  It  is  a  random-access  machine  based  on  the  one  described  in  [Aho-Hopcroft- 
Ullman  (1974)],  but  each  storage  location  holds  a  real  number  instead  of  an  integer.  The 
following  is  a  list  of  operations  available  at  unit  cost' 

(1)  The  arithmetic  operations  +  ,  -  ,  x  ,  /. 

(2)  Comparisons  between  two  real  numbers  <,s,=,<>,s,>. 

(3)  Functions  such  as  trigonometric,  exponential  and  logarithmic  functions. 

(4)  Finding  roots  of  a  fixed  degree  polynomial  equation. 

(5)  Indirect  addressing  of  memory. 

This   model   basically   embodies   the   capabilities   provided    by    any   high-level   programming 
language. 

2.2.  Upper  Bounds  and  Lower  Bounds 

We  will  use  a  notation  suggested  by  Knuth  (1968)  to  describe  asymptotic  upper  bounds 
and  lower  bounds  on  complexity  of  an  algorithm  or  a  problem:  0(J{n))  denotes  the  set  of  all 
functions  gin)  such  that  there  exist  positive  constants  c  and  n^  satisfying  \g{n)\  s  cf(n)  for 
all  n  ^  riQ.  n(/(n))  denotes  the  set  of  all  functions  g{n)  such  that  there  exist  positive 
constants  c  and  n^  satisfying  |^(n)|  ^  cf{n)  for  all  n  ^  n^. 

0(J(n))  represents  functions  that  are  at  most  as  large  as  /(n)  within  some  constant 
factor.  So,  an  upper  bound  of  0{J{n))  on  the  running  time  of  a  problem  is  established  by  an 
algorithm  that  solves  the  problem  of  size  n  in  cf{n)  time,  for  some  c  >  0.  On  the  other 
hand,  n(/(n))  represents  functions  that  are  at  least  as  large  as  /(i).  So  it  is  precisely  the 
notation  for  lower  bounds  of  the  running  time  of  a  certain  problem  or  algorithm. 

Vi'e  can  talk  alternatively  about  the  storage  needed  for  an  algorithm  to  solve  some 
problem.    The  notation  is  the  same  as  for  the  time  complexity. 

2.3.  The  Representation  of  the  Problem  and  Its  Solution 

Consider  the  following  problem: 

Problem  2.1.  Given  a  set  5  of  n  points  (in  the  first  quadrant),  find  the  smallest 
corner  triangle  Aoab  such  that  Aoab  D  5,  o  is  the  origin,  and  a  and  b  are  on  the 
positive  X-  and  y-axes  respectively. 


The  logarithmic  cost  criterion  for  an  integer  RAM  assigns  a  greater  cost  to  manipulations  for 
larger  integers  than  for  smaller  integers.  But  for  a  real  RAM  it  makes  more  sense  to  use  the 
uniform  cost  criterion. 


10 


Figure  2.1. 


Wc  can  represent  the  points  as  a  list  of  real-number  x-  and  y-  coordinates: 

S  =  {p,  =   ix„  >,),  X,  ^  0,    and  >,  ^  0,  0  s  ,  <  n  }. 

The  solution  can  be  represented  by  the  real  coordinates  of  a  and  b: 

a  =   (r,  0)     and     b  =   (0,  h),    z,  w  >  0 

Sec  Figure  2  1.    The  size  of  ^oab  depends  on  the  metric  used.    The  r»o  metrics  of  interest 
arc 


(1)      Area  measure:  \^oab\  = 


zw 


(2)     Perimeter  tneasurc:  |Aojii  =  ; -►  «■  ^  ^z--^w 


V7. 
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'n-l 


Figure  2.2. 

More  complicated  geometric  objects  will  be  represented  in  terms  of  points.  For 
example,  a  line  segment  will  be  represented  as  a  pair  of  points.    A  simple  polygon  P  will  be 

represented  by  a  list  of  its  corners  (vq,  v, *'n-i)  i°  counterclockwise  order,  with  no 

three  consecutive  corners  coUinear.  The  list  begins  with  the  corner  that  has  least  y- 
coordinate.-    See  Figure  2.2. 

2.4.    Preliminaries 

The  smallest  corner  triangle  problem  as  stated  is  a  continuous  one  in  the  sense  that 
there  are  infinite  number  of  candidates  for  the  solution.  A  finiieness  criterion  limits  the 
solution  of  the  problem  to  a  finite  set  of  objects,  so  a  finite  algorithm  can  be  constructed  to 
find  the  solution  from  the  finite  set.  In  order  to  solve  the  problem,  we  need  to  find  such  a 
finite  criterion.    Let  us  start  with  a  definition  of  "smallness:" 

Definition.  A  comer  triangle  C  is  a  local  minimum,  if  there  exists  e  >  0  such  that  the  size  of 
C  is  always  greater  than  the  size  of  C.  for  each  corner  triangle  C  that  contains  S  and  lies  at 
Hausdorff  distance  <  t  away  from  C. 

This  definition  is  not  immediately  useful  in  developing  an  algorithm  for  finding  the 
smallest  corner  triangle,  because  it  does  not  limit  the  solution  to  a  finite  set.  In  order  to 
proceed,  we  need  to  find  a  finite  criterion.    We  first  introduce  the  concept  of  the  convex  hull 


*  If  two  corners  have  the  same  least  v-coordinate,  the  list  begins  with  the  one  that  has  smaller 
x-coordinate.  This  is  the  same  as  the  standard  form  described  in  [Shamos  (1978)].  It  should 
also  be  noted  that  a  simple  polygon  in  other  representations  can  be  converted  to  the  standard 
form  in  linear  time. 
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of  a  set  of  points: 

Definition.  The  convex  hull  CH{S)  of  a  set  of  points  5  is  the  intersection  of  all  convex  sets 
containing  it. 

The  convex  hull  of  a  planar  point  set  is  a  convex  polygon  [McMullen-Shepard  (1971)] 
which  can  be  found  in  0{n  log  n)  time  [Preparata-Hong  (1977)].    Let  e,  be  the  directed  edge 

of  P  =  (vp,  v, V,.,)  from  v,  to  v,,^^rno<i  n)'  0  s  J  <  n.    Each  edge  <r,  has  an  angle  ^e, 

which  is  the  angle  measured  counterclockwise  from  the  positive  x-axis  to  ?,.  A  polygon  is 
convex  if  and  only  if  in  standard  form  0  s  ^ie^  <  ^e^,  .  .  .  ,  <  ^e„_^  <  Itx  [Shamos 
(1978)]. 

Lemma  2.1.  ab  overlaps  a  tangent  to  CH{S).  Furthermore,  if  v,  is  the  corner  of  CH{S) 
with  the  largest  x-coordinate ,  v,  with  the  largest  y -coordinate,^  then  any  corner  of  CH{S)  on 
ab  is  counterclockwise  between  v,  and  v,. 

Proof.  First  of  all,  by  the  definition  of  convex  hull,  any  convex  set  that  contains  5  also 
contains  CH{S).  Thus  the  corner  triangle  Aoab  contains  CH{S),  since  it  is  convex  and 
contains  S.  Furthermore,  since  we  are  concerned  with  finding  the  smallest  corner  triangle, 
we  can  assume  that  the  line  segment  ab  is  a  tangent  to  CH{S).  In  other  words,  ab  passes 
through  a  corner  or  two  adjacent  corners. 


Figure  2.3. 


'  If  there  arc  two  corners  with  the  same  largest  x-coordinatc  (v-coordinaie),  then  we  take  the 
one  with  larger  y-coordinate  (i-coordinate). 
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Refer  to  Figure  2.3  for  discussion.  Let  v,  be  a  corner  of  CH(S)  which  is  also  on  ab. 
From  the  convexity  of  CH(_S),  if  v,  is  not  between  v,  and  v„  then  clearly  v,  is  dominated  by 
cither  v^  or  v,.  (A  point  is  dominated  by  the  other,  if  both  of  its  coordinates  are  strictly 
smaller  than  the  other's.)  Since  w,  and  v,  are  both  on  or  under  ab,  the  domination  implies 
that  V,  is  under  ab ,  a  contradiction.  Thus,  any  corner  of  CH{S)  on  ab  must  be 
counterclockwise  between  v^  and  v,.    D 

For  the  rest  of  the  chapter,  we  assume  that  the  convex  hull  of  5  is  the  polygon 
''  ~  (^0'  ^1'  ■  •  •  *'m-  i)'  ^/  ~  (*.'  y)  ^°^  '  ~  0,  .  .  .  ,  m  -  1,  m  s.  n ,  and  v^  (resp.  v,)  is  the 
corner  of  P  with  the  largest  x-  (resp.  y-)  coordinate. 

Notation.  For  each  edge  e^  =  v^v,,,  in  P,  r  ^  /  <  r,  let  Z.  be  the  line  passing  through  v,  and 
V,.,,  and  c,  d  be  the  intersections  of  L  with  x-  and  y-axes  respectively.  We  use  £,  to  denote 
cd. 

2.5.    Finiteness  Criteria  for  the  Smallest  Corner  Triangle  Problem 

2.5.1.    The  Area  Metric 

Definition.    For  any  tangent  cd  of  P  such  that  c  is  on  the  positive  T-axis  and  d  is  on  the 

positive  y-axis,  we  say  that  cd  is  area  low  (high)  if  cd  HP  C  cm  (md),  where  m  is  the  center 
of  cd.  Otherwise  {  m  }  Q  cdDP  and  cd  is  said  to  be  area  balanced.  We  will  just  say  low, 
high,  and  balanced  for  short,  when  the  metric  is  understood  to  be  the  area.  In  particular,  for 
r  ^  1  s  r,  if  m,  is  the  center  of  £,,  then  £,  is  low  {high)  if  v,.,  (v,)  is  below  (above)  m,. 
Otherwise  E,  is  balanced. 
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Figure  2.4. 


Lemma  2.2.  ab  is  balanced 

Proof.  If  ab  is  high  (rcsp.  low),  then  by  perturbing  the  line  segment  clockwise 
(counterclockwise)  around  p  the  lowest  (highest)  point  of  intersection  with  P,  we  can  reduce 
the  area  of  Aoab,  contradicting  the  assumption  that  ^oab  has  the  smallest  area.  (See  Figure 
2.4.)    Wc  conclude  that  ab  must  be  balanced     D 

Lemma  2.2  gives  us  a  finiteness  criterion  from  which  we  can  derive  an  algorithm:  P, 
the  convex  hull  of  5,  will  be  computed  first,  ab  may  contain  an  edge  of  P  or  it  may  pass 
through  a  corner  of  P.  Ihcre  arc  linear  number  of  edges  and  corners  in  P.  If  the  balanced 
tangent  in  each  situation  can  be  verified  or  found  in  constant  time,  then  we  can  solve  the 
smallest  corner  triangle  problem  in  0(n)  time,  by  an  exhaustive  search  for  the  balanced 
tangent  The  next  lemma  gives  an  additional  property  of  tangents  of  the  convex  hull  of  5 
which  enables  us  to  use  a  binary  search  instead  of  the  exhaustive  search  for  the  balanced 
tangent: 

Lemma  2.3.  The  sequencr  of  flush  tangents  [E,.  .E     ,]  of  P  consists  of  a  sequence  of 

low   tangents,  followed  by  at  most  one   balanced  tangent,  followed  by  a  sequence   of  high 
tangents 

Proof    The  lemma  is  the  consequence  of  the  following  facts: 

(1)  If  £,  is  high  or  balanced  then  £,.  .  is  high 

(2)  If  £,  is  low  or  balanced  then  £,_  ,  is  low.  D 
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Lemma  2.4.   We  can  determine  whether  a  flush  tangent  is  high,   low,  or  balanced  in  0(1) 
time . 

Proof.    Consider  a  tangent  £,  of  P,  r  s  i  s  r.    £_  intersects  the  v-axis  at  (0,  y)  where 


y  = 


'.-1    *, 


Then,  £,  is 

(1)  low,  if  >|-|  is  smaller  than  -r, 

(2)  high,  if  >•,  is  greater  than  ■^,  and 

(3)  balanced,  otherwise. 

It  is  clear  that  0(1)  time  is  needed  for  the  computation.   D 


Lemma  2.5.  //£,_|  is  low  and  £,  is  high,  then  there  exists  a  balanced  tangent  that  passes 
through  V|. 

Proof.    By  a  continuity  argument.    D 

2.5.2.    The  Perimeter  Metric 

For    the    perimeter    metric,    the    smallest   corner   triangle    is   characterized    by    a   more 
complicated  condition  than  the  bisecting  property  for  the  area  case 


Figure  2.5. 
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Definition.    Let  cd  be  a  tangent  of  P  such  that  c  and  d  arc  on  the  x-  and  y-&xe.s  respectively. 
Let  4i  be  the  angle  ^dco  and  >  the  angle  ^odc.    Let  m  be  the  unique  point  such  that 


cm  cot  ^^ 
2_ 

md  cot  -^ 
2 


=  1. 


Wc  will  call  m  the  {perimeter)  center  of  cd.  And  wc  say  that  cd  is  perimeter  low  (high)  if 
cdnP  Q  cm  (md).  Otherwise  cd  is  said  to  be  perimeter  balanced.  We  will  just  say  low, 
high,  and  balanced  for  short,  when  the  metric  is  understood  to  be  the  perimeter.  See  Figure 
2.5.  In  particular,  for  a  flush  tangent  £,  of  P,  if  v,.,  (rcsp.  v,)  is  below  (rcsp.  above)  c  ,  the 
center  of  £,,  then  £,  is  low  (resp.  /i/^/i).    Otherwise  £,  is  balanced. 

Lemma  2.6.  at  is  balanced. 

Proof.  Assume  for  the  contradiction  that  ab  is  cither  low  or  high.  First  consider  the  case 
that  ab  IS  low  and  let  the  lowest  point  of  abHS  be  p .    See  Figure  2.6  . 


b 
b 


a       a 


Figure  2.6. 

Let  cd  be  a  line  seqment  such  that  p  i  cd  and  c  and  d  axe  on  the  positive  x-   and  v-   axes 
respectively.    Let 

J  =  pc,     g  =  pd. 


<t>  =    ^dco,     -y  =    ^  ode . 
If  L  denote  the  perimeter  of  Aacd,  then  we  have 


/cot^ 
dL    ^ 2_ 

'        /?  cot  -f- 


17 


(sec  Lemma  4.1  for  the  derivation).    Since  ab  is  low,  it  is  easy  to  see  that 

—         ^ 

ap  cot  -^-- 


d-i  .  y' 

>  "  >         pfc  cot  -^ 
^  2 

where  y'  =  ■^oba  and  <t>*  =   ^bao. 

So  by  a  infinitesimal  counterclockwise  rotation  of  ab  around  p ,  we  can  decrease  the 
perimeter  of  Aoab  while  maintaining  the  property  that  A^oab  D  S.  Similar  argument  applies 
to  the  case  where  ab  is  high.    Hence,  we  conclude  that  ab  must  be  balanced.    D 

Lemma  2.7.  The  sequence  of  edges  in  the  standard  form  of  P  consists  of  a  sequence  of  low 
edges,  followed  by  at  most  one  balanced  edge,  followed  by  a  sequence  of  high  edges. 

Proof.    The  lemma  follows  from  the  following  facts: 

(1)  If  £,  is  low  or  balanced  then  E,_,  is  low. 

(2)  If  £j  is  high  or  balanced  then  £,.,  is  high. 
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Figure  2.7. 

To  sec  (1),  draw  a  line  segment  cd  through  the  center  m,  of  £,  parallel  to  £,_,  such  that 
c  and  d  are  on  the  jt-  and  >-axes  respectively.  See  Figure  2  7  Let  m  be  the  center  of  cd.  It 
is  easy  to  see  that  (a)  m  is  above  m  .  and  (b)  m  ,  the  center  of  £_,,  equals  nrnOE,  , .  The 
point  m,_|  is  above  m,  since  m,.,  is  above  m  and  m  is  above  m  Hence,  £_  ,  is  low.  A 
similar  argument  applies  to  (2).    o 


Lemma  2.8.  li  lakes  0(1)  lime  to  determine  whether  E,  is  low    high,  or  balanced. 


Proof     From    the   definition,   it   takes   only   a   few    arithmciic   operations   and   trigonometric 
function  evaluations.    D 
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Lemma  2.9.  If  E^_^  is  low  and  £,  is  high  then  there  exists  a  balanced  tangent  that  passes 
through  V|. 

Proof.    By  a  continuity  argument.    □ 

2.6.    Aa  Algorithm  for  Both  Area  and  Perimeter  Metrics 

It  is  very  nice  that  wc  get  essentially  the  same  finitcncss  criterion  for  both  metrics,  area 
and  perimeter/  except  for  the  difference  in  the  definitions  of  "low,"  "high,"  and  "balanced." 
So  the  two  problems  can  be  solved  with  the  same  algorithm. 

To  solve  the  problem,  we  first  compute  P,  the  convex  hull  of  S.  Then  we  proceed  to 
find  the  two  comers  v,  and  v,  of  P.  That  can  be  done  by  some  kind  of  binary  search.  By 
comparing  the  y-coordinates  of  two  consecutive  corners  v,,  v,^.  of  P,  i  =  m/2  (m  is  the 
number  of  comers  in  P),  wc  can  determine  whether  v^  is  before  or  after  v,.  After  v  and  v 
arc  found,  another  binary  search  is  then  performed  to  find  the  balanced  tangent; 

Algorithm  2.1.  The  Smallest  Corner  Triangle 

Step  1       Find  the  convex  hull  P  of  the  planar  point  set  S . 
Step  2       Use  a  binary  search  to  find  the  corners  v^  and  v,  of  P. 

Step  3       Perform  a  binary  search  on  the  edges  that  are  counterclockwise  between  e,  and 
f,_i  to  find  the  balanced  edge,  if  it  exists.  Otherwise,  find  the  last  low  edge. 

Step  4       If  the  balanced  edge  exists,  say  <•,  =  v\v,. ,,  then  ab  overlaps  e ^.    Thus 

-«'Xv,-i->',) 


X, ,  0 

>',-i-.v, 


and  b 


0,y,- 


-t,-i--f. 


Step  5       Otherwise,  there  is  no  balanced  edge.  Let  *,_,  be  the  last  low  edge  and  e^  be  the 
first  high  edge. 

(1)  The  minimum  area  corner  triangle  is  determined  by  a  tangent  ab  of  P  which 
passes  through  v,  and  is  bisected  by  v,,  where 

a  =   (2x,,  0)  and  b  =   (0,  2y). 

(2)  The    minimum    perimeter    corner    triangle    is    determined    by    a    perimeter 
balanced  tangent  ab  through  \\  where^ 


But  as  we  will  sec  in  the  later  chapters,  that  is  not  always  the  case. 
The  formula  is  taken  from  [Oman  (1976)]. 
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yXx,*i2x.yy-)  ^,0,  ^  (2x,y,)^ -) 

a  =  (x  -^ : — ,  0)  and  b  =    (0,  y  +  : — .  D 

'         yM2x,yy-  x,^i2x,y,)^- 

In  conclusion,  wc  have 

Theorem  2.10.     The    smallrsi    corner    trianglf    containing    a    set    of  n    points    in    the  first 
quadrant  tn  the  plane  can  be  found  in  0{n  log  n)  time  and  0{n)  space. 

Proof.  The  computation  of  convex  bull  in  Step  1  takes  0{n  log  n)  time.  It  takes  0(1)  time  to 
determine  whether  an  edge  is  low,  balanced,  or  high,  so  the  binary  searches  in  Steps  2&3 
take  0(Iog  n)  time.  Steps  4&5  each  takes  only  constant  time.  Thus  the  problem  can  be 
solved  in  0{n  log  n)  time  Note  that  if  a  convex  polygon  is  given  instead  of  a  point  set,  then 
the  problem  can  be  solved  in  O(log  n)  time.    □ 

Note  that  if  we  have  the  convex  hull  already,  then  we  can  solve  the  problem  in  0(log  n) 
time. 


CHAPTER  3 


Minimum  Area  Enclosing  Polygons 


In  Chapter  2,  wc  examined  a  simple  case  of  enclosure  problems.  Now  we  consider  an 
enclosure  problem  that  has  more  practical  implications: 

Problem  3.1.  Given  a  convex  n-gon  P  and  a  fixed  i,  3  <  *  <  n,  find  a  convex  k- 
gon  Q  that  encloses  P  and  has  the  minimum  area  among  all  convex  polygons  that 
enclose  P . 

For    the    rest    of    the    chapter,    we    will    assume    that    P  =  (vq,v, v„_,),    and 

Q  -  (wq.h-i,  .  .  .  ,Wj.,)  (the  corners  are  listed  counterclockwise).  Let  *,  denote  the  directed 
edge  of  P  from  v,  to  v,.,,  (0  s  /  <  n).  Similarly,  let  j,  denote  the  directed  edge  of  Q  from 
w,  to  H,.  |,  (0  s  I  <  V).    The  subscripts  of  v,  and  h,  are  taken  mod  n  and  mod  *  respectively. 

What  characterizes  the  solution  polygon  2?  Since  we  are  concerned  with  a  minimum 
area  enclosing  t-gon,  we  can  assume  without  loss  of  generality  that  every  edge  of  Q  comes  in 
contact  with  P  Each  edge  of  C  is  a  tangent  of  P .  Define  a  tangent  7  of  P  to  be  extremal  or 
flush  if  it  passes  through  two  adjacent  corners  of  P  (thus  overlapping  an  edge  of  P).  We 
write  £,  for  the  tangent  that  overlaps  edge  f,.  We  will  assume  that  every  T  is  directed  so  that 
P  is  to  the  left  of  T.  Let  T'  denote  the  closed  half-plane  to  the  left  of  T.  It  is  clear  that  a 
minimum  area  convex  t-gon  Q  enclosing  P  must  be  the  intersection  of  k  half-planes 
determined  by  k  distinct  tangents  T^,  T.,      .  .  ,T^o{  P .    The  t-gon  Q  is  given  by 

* 

n  T- 


3.1.    A  Finiteness  Criterion  and  a  Nai\c  Algorithm 

In  this  section,  we  will  give  a  finiteness  criterion  for  a  minimum  area  enclosing  k-gon  Q 
and  discuss  its  usefulness  for  an  algorithm.  First  we  give  some  definitions  needed  for  stating 
the  finiteness  criterion. 
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an  -4 -pair 


a  V-pair 


Figure  3.1. 

DeTinition.  For  an  ordered  pair  of  directed  lines  (Z, ,  M),  if  the  angle  counterclockwise  from 
i.  to  Af  is  greater  than  -n,  then  we  say  {L,  M)  is  V-shaped.  Otherwise  we  say  (/. ,  M)  is  A- 
shaped.  Sec  Figure  3.1.  An  edge  j,  \n  Q  \%  &  V-edge  (resp.  an  ;4-edgc)  if  its  two  neighboring 
edges  (i,_,,  J,,,)  are  V-shapcd  (resp.  A-shapcd).  We  call  (5,_,,  J,,,)  supporting  lines  of  s,. 
We  call  (j_|,x,.,)  an  A-pair  or  a  \'-pair  accordingly.  (Remark;  'V  and  'A'  are  chosen  for 
the  shapes  of  these  letters,  being  mnemonic  for  the  orientations  of  the  supporting  lines  -  if 
we  imagine  the  line  s^  as  horizontal  and  directed  to  the  left  and  P  as  sitting  below  s^.) 


Figure  3.2.  Tangent  7  and  edge  cc   embedded  in  an  A-butterfly 


Dernition.  Consider  an  edge  5  of  ^ .  Let  m  be  the  center  of  i,  =  h^v,  .  ,.  If  P  Omlv,. ,  =  0 
(m»v,.i  denotes  the  closed  interval  between  m  and  «  .  ;).  then  we  say  that  j  is  early. 
Similarly,  if  /Timlv,  =  0,  then  it  is  late  Otherwise  m  (  s^r\P  and  we  say  that  s  is 
balanced.    Let  T  be  the  tangent  of  P  that  overlaps  s  .    We  say  that  j,  is  extremal  or  flush,  if  T 
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is  extremal.    (Note  that  "extremal"  or  "flush"  is  defined  with  respect  to  P .) 

Definition.  Consider  an  edge  sofQ.  Let  T  be  a  tangent  of  P  that  overlaps  5,  and  passes 
through  a  corner  of  P,  say  v^.  Let  S  and  U  be  the  two  tangents  of  P  that  overlap  j,.  ,  and 
J,.  I  respectively.  Let  S  intersect  Ej  and  f^.,  at  a  and  b  respectively  and  U  intersect  Ej  and 
£  .  I  at  a  and  b  respectively.  We  call  the  region  B,  consisting  of  the  rwo  triangles  AabVj  and 
A  a  b  Vj,  a  butterfly.  (See  Figure  3.2.)  The  edge  j,  of  Q  is  said  to  be  embedded  in  the 
buncrfly  B^.  Furthermore,  B^  is  an  A-butterfly  if  j,  is  an  M-edge  and  a  V-butterfly  if  j,  is  a  V- 
cdge.    Note  that  the  butterflies  are  defined  with  respect  to  the  polygons  P  and  Q. 

In  the  following,  we  give  three  lemmas  that  characterize  the  solution  polygon  Q: 

Lemma  3.1.  All  edges  in  a  minimum  area  enclosing  k-gon  Q  are  balanced.    Furthermore, 
every  A -edge  is  flush. 

Proof.  Consider  an  edge  s,  of  Q  embedded  in  a  butterfly,  say  fl^,  0  ^  1  <  t,  0  s  ;  <  n.  The 
edge  J,  extends  to  become  a  tangent  T .  Since  Q  has  minimum  area,  the  area  of  BJ  =  7"  ns, 
as  a  function  of  7  is  at  a  minimum.  Refer  to  Figure  3.2  for  discussion.  BJ^  is  the  union  of 
triangles  Aabv  and  ^abv  .  Consider  the  case  where  B,  is  a  V-butterfly.  First  assume  for 
the  sake  of  contradiction  that  j,  is  not  balanced.  So  j,  is  cither  early  or  late.  If  it  is  early, 
then  clearly  a  counterclockwise  rotation  of  j,  would  result  in  a  smaller  area  for  Q, 
contradicting  the  optimality  ol  Q .  A  similar  argument  applies  in  the  case  of  a  late  edge 
Thus  a  V-edge  in  Q  must  be  balanced. 

Next,  we  consider  the  case  where  fl,  is  an  A-butterfly.  If  j,  is  early,  then  we  can  rotate 
5,  counterclockwise  around  the  center  of  s,.  The  area  of  Q  will  be  reduced  (or  kept 
unchanged  when  the  rwo  neighboring  edges  of  5,  arc  parallel)  by  the  rotation.  In  any  event, 
now  s^  does  not  intersect  P  any  more,  contradicting  the  optimality  of  Q.  Similar  arguments 
apply  for  the  case  of  a  late  edge  and  the  case  of  a  balanced  but  non-flush  edge.  So  it  is  clear 
that  J,  must  be  flush  and  balanced  Thus  every  A-tdgt  of  Q  must  be  both  balanced  and  flush 
D 

Lemma  3.2.    There   are    at   most  Iko   V -edges   in   a   convex  polygon    with    more    than    three 
edges.    Furthermore,  when  there  are  exactly  two  V -edges,  they  are  adjacent. 

Proof.  Let  e,  and  e  be  two  distinct  V-edges  in  a  convex  n-gon  P .  From  our  definition  of  V- 
edges,  we  know  that  the  angle  a  counterclockwise  between  f  _,  and  ?,.,  is  greater  than  r. 
and  so  is  the  angle  3  between  e  ,  and  ?,.,.  But  since  P  is  convex,  the  orientations  of  the 
edges  of  Q 

are  in  order,  so  a  and  3  must  overlap.    Thus  either;  =  1  -  l(mod  n)  or  /  =  _;  -  1  (mod  n),  i  c 
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e,  and  e   arc  adjacent.    It  is  also  clear  that  there  can  not  be  more  than  two  V-cdgcs.    D 

From  the  last  lemma,  we  know  that 

(1)  There  are  at  least  k-2  ^ -edges,  and  thus  k-2  flush  edges  in  Q . 

(2)  If  there  are  two  V-cdges  in  Q  then  they  must  be  next  to  each  other. 

Property  (2)  allows  us  to  strengthen  the  criterion  by  showing  that  there  are  in  fact  k-  1  flush 
edges  in  Q: 

Lemma  3.3.  For  k  ^  3,  there  exists  a  minimum  area  kgon  Q  that  has  at  least  k-1  edges 
flush  with  P. 

Proof.    Assume  that  the  two  adjacent  V-edges  of  Q,  say  j,_[  and  s^,  are  not  flush. 


^,-i('*',) 


Figure  3.3. 


Refer  to  Figure  3.3  for  discussion.  Let  5.  ,(m  )  (resp.  •S|.|(*v,))  denote  the  line  that 
passes  through  h,  and  has  the  same  direction  as  j  ,  (rcsp.  J,.,).  Let  us  perturb 
simultaneously  5,. ,  and  s  by  rotating  them  around  their  respective  centers  c,_ ,  and  c,,  in  such 
a  way  that  j,_ ,  becomes  j  _ ,,  j    becomes  i, ,  and 


-•,  ^  !,-:(»■,)    nj,.,(HO    , 


n.v 


A  "gain-and-loss  analysis'  shows  that  after  the  perturbation,  the  area  of  Q  is  decreased  by  an 
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amount  equal  to  the  area  of  the  diamond-shaped  region  bounded  by  j,_|,  j  ,  s,_-i{w),  and 
^i-i(^i)-  (The  decrease  is  essentially  due  to  the  fact  that  (j,.^,  5,,,)  is  V-shaped  and  i,,,  and 
s,  arc  balanced.)  This  contradicts  the  optimality  of  Q  and  proves  that  at  least  one  of  the  two 
V-edgcs  must  be  a  flush  edge.  Thus  we  have  at  least  k-  1  flush  edges  in  Q.  This  lemma  also 
appears  independently  in  [DcPano  (1984)].    D 

In  summary,   a  minimum  area  i-gon  Q   enclosing  a  given  n-gon  P   has  the  following 
properties: 

(1)  At  least  A-  1  edges  of  Q  are  flush  with  P. 

(2)  All  midpoints  of  edges  of  Q  are  in  P . 

From  this  finiteness  criterion,  we  can  device  a  brute-force  optimization  algorithm  as 
follows:   The  algorithm  will  guess 

(1)  the  t -  1  edges  of  P  that  are  flush  with  Q ,  and 

(2)  from  the  t-  1  edges,  two  neighboring  edges  to  put  the  last  (possibly  non-flush)  edge  in 
between. 

Each  guess  represents  a  possible  "configuration"  of  the  solution.  For  each 
configuration,  the  algorithm  computes  the  last  edge  that  minimizes  the  area  of  Q'  and  verifies 
that  the  optimality  criterion  above  is  satisfied.  This  is  done  in  0(n)  time  for  each  guess. 
This  naive  algorithm  runs  in  exponential  time  while  using  only  linear  storage.  Can  we  do 
better?  Can  we  trade  more  storage  for  a  faster  running  time?  In  the  following  sections,  we 
will  answer  the  question  positively.  By  exploiting  a  dynamic  programming  technique  and 
finding  more  geometric  structure  in  the  minimum  area  enclosing  polygons,  we  can  develop  a 
much  faster,  polynomial-time  algorithm. 

3.2.    The  Dynamic  Programming  Approach  to  Optimization 

TTie  fault  of  the  naive  algorithm  described  in  the  previous  section  is  that  it  exhausts  all 
possible  configurations  of  the  solution  and  solves  the  problem  as  a  whole.  One  configuration 
may  be  just  slightly  different  from  the  other,  yet  they  arc  all  recomputed  from  scratch.  So 
there  is  a  wasteful  amount  of  redundant  work.  Fortunately,  dynamic  programming  [Bellman 
(1957)]  comes  to  the  rescue.  In  essence,  a  dynamic  programming  algorithm  calculates  the 
solutions  to  all  subproblems;  proceeding  from  smaller  subproblcms  to  the  larger  problem. 
All  the  partial  solutions  to  subproblcms  are  stored  in  tables.  The  advantage  of  dynamic 
programming  lies  in  the  fact  that  once  a  subproblem  is  solved,  the  solution  is  stored  and  can 
be  used  later,  without  recalculation,  in  solving  other  subproblcms. 

Specifically,  for  our  problem,  instead  of  computing  the  minimum  area  enclosing  k-gon 
as  a  whole,  we  will  compute  sequences  of  consecutive  edges  which  arc  part  of  the  minimum 
area  enclosing  *-gon  Q.    Before  we  define  these  sequences  precisely,  we  need  the  following 

'  This  is  essentially  the  same  as  computing  V|(/,  ;')  and  ©,(/,  ;)  described  in  the  next  section. 
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notatioD: 


NoUtion.    For  n  ^  3,  and  0  s  a,,  a,, 


.  .a,  <  n,  3  s  I  <  n,  wc  write 


a^  <  a^  < <  a, (mod  n) 

if  there  exists  a  number  /,  0  s  /  <  n,  such  that 

a|  +  /(naod  n)  <  3^  +  /(mod  n)  <  ,  .  .  .  ,  <  a,-^  /(mod  n). 

NoUlion.    Let  E(i)  denote  the  set  of  all  tangents  of  P  that  pass  through  v,  but  not  v.,.    Let 
E'  =  U  =(()•    It  is  clear  that  every  tangent  belongs  to  a  unique  E((),  0  s  /  <  n. 

1-0,  .1-1 

For     instance     £,  =  v.-i  ^  c(i)-      For     a     tangent     7  €  c(<),     the     direction     of     T     is 

counterclockwise  between  those  of  £,_,  and  £,.    For  a  sequence  of  tangents  (7",,  7", T„) 

of  P,  if  their  angles  are  in  order  in  a  counterclockwise  traversal  of  the  unit  circle,  then  we 
say 

T^■s  T.^ ^  T„(mod  P). 

Furthermore,  if  they  are  all  distinct,  we  write 

7",  <  r,  < <  r„(mod  P). 


r:^ 


3  , 


Figure  3.4.  An  (i,  ;)-chain  of  length  3 


Dennition.    Let  e^  and  r,  be  two  distinct  edges  of  P.    For  0  £  h  s,  i-3,  we  define  an  h-sided 
(i,  j)-chain  (or  simply,  a  chain,  if  h,  i,  and;  arc  understood)  to  be  any  sequence  of  tangents 

C  =  (7',,  7, 7^)       such       that       for       every       m.       1  ^  m  s  /i,       7^  €  5(u^)       and 

/  <  i<i  <  u,  < <  u^  <  ;(mod  n).     N^'c   define    the    'krea   of  such   a   chain   C   to   be   the 

area  of  the  polygon 
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£*  (lEj'  DCjC,   n 


',6,^-)' 


where  CjC,  is  the  line  directed  from  c^,  the  center  of  e^,  to  c,,  the  center  of  «,.  See  Figure  3.4. 
The  chain  C  is  flush  if  all  its  tangents  are  flush  with  P . 

Notation.  Let  W^{i,j)  denote  the  set  of  all  A-sided  (j,  y)-chains  and  let  VV^((,  ;)  denote  the 
set  all  /i-sided  flush  (i,  y)-cbains.  Wl{i,  j)  Q  W>,(<,  ;)  and  Wf,(i,  j)  is  empty  if  there  are  less 
than  h  corners  couterclockwise  between  v,  and  v,  (exclusively).  We  use  C^(i,  j)  to  denote  any 
chain  with  minimum  area  among  all  chains  in  Wi^{i,j),  and  call  it  an  optimal  chain  for 
W^{i,j).  Its  area  is  denoted  by  a^(i,j).  Similarly,  we  use  Cl{i,j)  and  a^(i,  ;)  to  denote 
any  optimal  flush  chain  for  W^(i,  j)  and  its  area,  respectively. 


/\P'.J 


Figure  3.5.  a^,(i,  J)  =  Area{^c,CJp ,  ,) 


Let  /3,  denote  the  intersection  of  £,  and  Ej.  Then  the  area  of  to(''  -')  (  ~  '^o(''  j))  '^ 
determined  by  that  of  Ac^c^p,  ,: 

ao(',  ;)   =   ao('.  ;■)  =  Area(Ac,c.p,  ), 

which  can  be  computed  in  0(1)  time.    See  Figure  3.5.    Obviously,  Qo(',  ;)  can  be  computed 
in  0{n')  time  for  all  /',  j. 

By  Lemma  3,3,  we  know  thai  there  are  *-  1  flush  edges  in  Q .  So  we  can  find  two  flush 
edges  s  and  s  in  Q  such  that  s  (resp  s  )  overlaps  E  (resp.  £,)  and  s  and  5  decompose  Q  into 
an  optimal  one-sided  (i,  y)-chain  and  an  optimal  (t-3)-sided  flush  (J,  i)-chain.  See  Figure 
3.6  for  an  example. 
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optimal  one-sided  flush  (2,  6)-chain 


optima!  one-sided  (6,  2)-chain 


Figure  3.6. 

So  Problem  3.1  can  be  reduced  to  the  r*o  following  subproblems: 

Problem  3.2.    For  0  ^  t,  j  <  n,  i  ^  j,   find  the  optimal  one-sided  (i,  7)-chain  for 

w'i('.;). 

Problem  3.3.    For  0  ^  i.j  <  n,  i  i=  j  and  0  <  /i  s  i-3,  find  the  optimal  /i-sided 
flush  ((,  7)-chain  for  H'^(i,  j). 

Hence,  if  we  compute  ir,(/,  ;)  and  Cl.i(i,  j)  and  their  areas,  for  all  i  *  j,  then  the  area 
of  Q  can  be  obtained  by 

min  {  ai-30.  0^  a|(/,  ;)  }. 

Consequently,  if  two  problems  mentioned  above  are  solved,  then  Q  can  be  computed 
with  0{n-)  additional  steps.  As  it  turns  out,  all  our  solutions  to  both  subproblems  require 
fl(n-)  time,  so  the  running  time  of  all  our  algorithms  is  determined  by  that  of  the  algorithms 
for  Problem  3.2  and  Problem  3.3. 


3.3.    More  Geometric  Structure 

In  order  to  solve  Problem  3.2  and  Problem  3  3  cfficicntiv,  we  need  to  know  more  about 
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the  structure  of  the  subproblems. 

Dennilion.      For    C     =     (T", T,)     i     W  ,(i ,  j)    and    C'     =     (r, TJ     €     1V,(;,;), 

i  <  j  <  J  (mod  n),  we  say  that  C  intersperses  C    if 

7",  :s  7",  s  r,  s  7,  s  ,  .  .  .  ,   s  7"^  s  r^(mod  P). 
In  particular,  for  C  €  ^'^(i,  ;)  and  C    €  W^*(i,  ;),  if  T,  =  £„  and  7",  =  £     for  1  s  /  s  A  and 

a,  s  a,  s  a,  s  a,  s  ,  .  .  .  ,   s  a^  s  fl^(mod  n), 
then  C  intersperses  C  ■ 

Lemma  3.4.  (Interspersing    Lemma    I)    //  /' s  y  s  ^' (mod  n),    then    Cl(i,j)    intersperses 
CliiJ). 

Proof.    We  will  show  the  interspersing  by  induction  on  h: 

Basis  {h=  1):  Let  C,  P  ^  W,((,  j)  be  (E,J  and  {E„)  respectively,  i  <  m'  <  m  <  ;(mod  n). 
Let  C  ,  P  €  ^^1(1,7)  be  (£„)  and  (£„  )  respectively,  /<  m'  <  /n,m'  <  _/(mod  n).  Refer  to 
Figure  3.7  for  discussion.  We  will  show  that  if  Area(C)  s  Area{P),  then 
Area{C  )  <  Area[D  ).'    Now,  it  is  easy  to  see  that  if 

Area{C  )  =  Area{C)*  y 

Area{V)  =  Area{r)  +  y^  5. 


'  In  a  recent  paper  Aggarwal,  Klawe,  .Moran,  Shor,  and  Wilber  proposed  an  algorithm  that 
takes  advantage  of  this  stronger  property  (rather  than  the  interspersing  property)  and  finds 
optimal  flush  chains  with  a  saving  factor  of  0(log  n).  Thus  they  were  able  to  find  the 
minimum  area  enclosing  i-gon  in  0{n-\og  k)  time  [Aggarw  al-Klawe-Moran-Shor-Wilber 
(1986)]. 
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Figure  3.7.  The  dotted  lines  rcpreseiit  edges  of  P 


Then  wc  have 

Area{C  )  =  Area{')  -  7 

S  Area{D)~  ^ 

<  Area(V)'  7^5 

=   Areai::). 

So  AreaiC  )  <  Area{7).  Lei  C](i.j)  and  C]{i,j)  be  {E ^)  and  {E ^.)  respectively  If 
i  <  a  s,  b  <  j  <  ]  (mod  n)  then  C'(i,  ;)  intersperses  C'(',  ;  ).  Assume  for  the  purpose  of 
contradiction  that  i  <  b  <  a  <  j  <  j  (mod  n).  Let  l)^  denote  the  chain  that  belongs  to 
W|((,_/)  and  shares  the  same  tangents  with  C'(i.j):  T,  =  (f.J  Similarly,  Pf.  denotes  the 
chain  that  belongs  to  W^[(i,  ;)  and  shares  the  same  tangents  with  C\(i.  j  ),  Df.  =  {£,).  NN'e  say 
that  Tfl  and  P^  arc  the  result  of  "edge  exchange"  of  C\(i ,  j)  and  C]ii,  j  ).  From  the  optimality 
of  C](i.  j),  we  have 

Area{C;{i.j))  ^  Arfa{V,) 

Since  C]U,j)  =  (£j).  Pt>  ^  '**'i('.  y)>  and  1  <  b  <  a  <  j  <  j  (mod  n).  wc  have 

Area{PJ  <  Arfa(C\0  ,  J  )), 

contradicting  the  optimality  of  C](i.j)-     We   call   the   argument  used  here   to  establish  the 
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basis,  the  "edge  exchange  argument." 

Induction  (>i>l):  Let  E„  be  the  last  tangent  of  C.  Similarly,  let  £„  be  the  last  tangent  of  C  . 
If  m  s  m'  <  /(mod  n)  then  by  the  induction  hypothesis,  the  (/i  -  l)sided  subchain  of  C  from 
i  to  m  intersperses  the  (/i-l)-sided  subchain  of  C  from  i  to  m ' .  Hence  C  intersperses  C  as 
desired.  Otherwise,  m'  <  m  <  j  (mod  n).  A  similar  edge  exchange  argument  can  be  used  to 
complete  the  proof.    D 

Observe  that  the  proof  for  the  last  lemma  does  not  depend  on  the  fact  that  the  chains 
involved  are  flush  chains.  Thus  it  is  valid  even  when  C;,((,  ;)  and  C^(/, ;  )  are  used  instead  of 
C'hU,  j)  and  Cl(i,j): 

Lemma  3.5.  (Interspersing    Lemma    II)    Let    i  s  y  s  _/  (mod  n).    If    C  =  Cf,ii,j).    and 
C    =   ^a('>  .'  )  ''""  ^  intersperses  C  ■ 

In  fact,  this  property  holds  for  optimal  chains  in  any  subsets  W ^{i,  j)  Q  ^/,(',  ;)  and 
'^^('.j)  C  ^'a('. /)  which  are  closed  under  "edge  exchange."  Perhaps,  it  is  possible  to 
exploit  this  to  obtain  a  more  efficient  algorithm.  However,  since  we  will  not  use  this 
observation,  we  omit  its  formulation.' 

What  the  interspersing  lemmas  mean  for  our  problems  is  that  if  we  compute  some 
optimal  chain  first,  then  other  optimal  chains  can  be  computed  more  efficiently  using  the 
information  provided  by  the  first  optimal  chain.  The  interspersing  lemmas  enable  us  to  limit 
the  range  of  search  for  the  edges  of  optimal  chains,  thus  speed  up  the  dynamic  programming 
process.* 

3.4.    Ad  Algorithm 

In  this  section,  we  exploit  the  interspersing  lemmas  to  solve  the  two  subproblems: 
Problem  3.2  and  Problem  3.3. 

3.4.1.    Computing  One-Sided  Optimal  Chains 

Recall     that     Problem     3.2     calls     for     us     to     find     C^{i,j]     and     a|(/,;)     for     all 

I,  J  =  0 n  -  1,  (  't  j.    The  value  of  Q|((  ,  ;)  is  represented  by  a  real  number.    We  choose 

to  represent  l,(/,  ;)  by  the  index  of  the  corner  of  P  that  it  passes  through  and  by  its  angle; 
V|((,  j)  and  e,((,  ;).  More  precisely,  if  t,((,  ;)  is  a  tangent  T  of  P ,  then  7  i  5(V,(i,  ;))  and 
0i('.j)  =  i^T.  If  W',(j,;)  =  0,  then  t ,(',;)  is  undefined  and  we  conventionally  let 
^\(',J)  =   '+1  and  a|(/,;)  =  ==.    We  can  view  Problem  3.2  as  a  problem  of  computing  the 


'  In  O'Rourke  et.  al.  [O'Rourke-Aggarw  al-Maddila-Baldwin  (1984)]  there  is  another  form  of 
interspersing  lemma  for  the  smallest  triangles  enclosing  a  given  polygon 

"  See  [Yao  (1980)]  for  another  way  of  speeding  up  dynamic  programming. 
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entries  of  three  n  x  n  matrices  a,,  V,,  and  0,. 

Algorithm  3.1.  Optimal  One-Sided  Chains 

Step  1       Let  1  =  0. 

Step  2       Let     ;    =  /+l(modn)     and     j    =   /-l(modn),      and      let      V'|(i,^)  =   i-l      and 

<*i('.y)  =  "•    Compute  V,((,  7  ),  and  a,((,)  )  by  using  a  slightly  modified  version 

of  Algorithm  2.1.' 
Step  3       Let;,  =    [(j  +y  )/2J. 
Step  4       If  (?,,  f^  )  is  an  ;4 -pair,  then  1 1(/,  ^i)  =    C[(i,  7,).    Then  a  ,((,  y',)  is  given  by 

min  (  aj,(',  *)  +  aQ(k,j^)->-AreaiAc,c^c,  )  } 

And    ^1(1,  _/,)  =  k    and    6|(i,  j,)  =  ■^f*    where    t    is    the    index    that    determines 

Step  5       Otherwise  (i.y,)  is  a  V-pair,  and  V,(/,  7,),  a^(i,j^)  can  be  computed  by  using  the 
modified  version  of  Algorithm  2.1. 

Step  6       Compute 

{  V, (/,;■),  a,(,-,y)  i;  =  /+l y,-l  }  and 

by  settmg  7    to  the  value  of  7,  and  applying  Steps  3-5  recursively. 
Step  7       Compute 

{  V, (/,;).  a, (,,7)  I;  =  7,-! 7' -  1  } 

by  setting  7    to  the  value  of  7,  and  applying  Steps  3-5  recursively. 

Step  8       Repeat  Steps  2-7  for  i  =    1 n  -  1.    D 

We  now  determine  the  running  time  and  storage  requirement  of  Algorithm  3.1.  Let 
^0  .7  .  '■)  denote  the  time  spent  on  computing  the  i-th  row  entries  from  column  7  to  column 
7  ,  where  r  is  the  number  of  edges  bcrween  V|(/,  7  )  and  V|(i,  7  ).    The  recursion  formula  for 

r(j  ,J  ,r)^   T(j  .7,.  r,)-r(7,.;   .  r-r.)-O(r), 


We  have  here  the  same  problem  as  Problem  2.1,  except  that  the  rv*  o  supporting  lines  e^  and 
e  arc  not  the  x-axis  and  v-axis  Using  Algorithm  2  1  and  taking  into  account  the 
orientations  of  e ^  and  e  ,  V ^{i ,  j  )  and  Q  \,i .  j  )  can  be  computed  in  logarithmic  time,  since  the 
convex  hull  of  (  v,. ,,  .  .  .  ,  w  ._    }  is  readily  available  to  us. 
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where  r,  is  the  number  of  corners  between  V^{i,j)  and  V^(i,j^).  If  there  are  s  corners  from 
*'■.]  '°  *'_i'  '^^°  ''  '^  ^^^y  ^°  verify  that  T(J  ,  j  ,  r)  =  0{r  log  s).  Since  r  and  5  are  both 
less  than  n,  this  proves  that  0(n  log  n)  time  is  sufficient  to  compute  all  the  entries  of  the  i-th 
rows  of  a,,  V|,  and  0,.  There  are  n  rows  in  each  matrix.  Hence  Algorithm  3.1  runs  in 
0{n-\og  n)  time  and  requires  0{n-)  space. 

3.4.2.    Computing  Optimal  Flush  Chains 

In  this  subsection,  we  proceed  to  give  an  algorithm  for  Problem  3.3;   computing  optimal 

flush  chains,  C^(i,  y)'s  and  their  areas.    Let  l^(<,  7)  =  [T^ T;,].    We  define  V^(i,  j)  to 

be  the  "index"  of  the  middle  tangent  of  Clii,j),  i.e.,  Vl{i,  j)  =  a  if  7"(aij-|  =  E^.  Again, 
we  regard  V^(i,  _/)  as  an  n  x  n  matrix.  V^((,  ;)  =  i+l  and  a^((,y',)  =  00  when  there  are  at 
most    A-1     edges    in    P     from    «,    to    e^    (exclusively).       Observe    that    if    the    matrices 

Vp  Vj V^  are  available,  then  the  optimal  chains,  Cl{i,  j),  can  be  easily  reconstructed  in 

0(h)  time  for  each  /',;,  0  s  i  ^  j  <  n.  Hence,  Problem  3.3  can  be  reduced  to  the  task  of 
computing  the  entries  of  matrices  Vl  and  q^,  for  all  h  between  1  and  k-3. 

Algorithm  3.2.  Optimal  Flush  Chains 

Step  1       Let  J  =  0  and  h  =  \. 

Step  2       Let  a  =    l-^!^^J  and  i  =    f-^1. 

Step  3       Let  j    =  i+l(modn)   and  j    =  i-l(modn),   and   Vl(i,j)  =  i+\,   a^(;, ;)  =  oc. 
°^'h{',J)  is  given  by 

min  a'((,  i)  +  a^(*,  ;  )  +  Area{^c^c^^c  ■). 

)'  <  k  <  i'(mod  11 

And  V^((,  7,)  is  given  by  the  number  k  that  provides  the  minimal  value. 
Step  4       Let;,  =    [0  'j)f2\. 
Step  5       a^(;,  _/|)  can  be  obtained  from 

min       ^  (  a*(/,  *)+ a^(*,;,)-/lrtfa(Ar,rjC,)  }. 

*  =  i'J!(i..;')       .y^i I.  ■  I 
And  V^'((,  _/,)  is  given  by  the  number  k  that  provides  the  minimal  value. 
Step  6       Compute 

{  v;(',y),a;(/,7)  I;  =  ;  -1 ;•;-!} 

by  setting;    to  the  value  of;,  and  applying  Steps  4-5  recursively. 
Step  7       Compute 
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{v;{i,j),a;u,j)  |y  =  ;,  +  ! ;  -u 

by  setting  ;    to  the  value  of  y,  and  applying  Steps  4  5  recursively. 

Sifp  8       Repeal  Steps  2-7  for  i  =    1 n  -  1  and  /i  =  2,  ...,*- 3.    □ 

We  now  show  Algorithm  3.2  runs  in  0(Jtn-log  n)  time.  Let  TiJ  ,  j  ,  r)  denote  the  time 
spent  in  computing  the  i-th  row  entries  from  column  y  to  column;  where  r  is  the  number  of 
edges  between  Vj|(i, ;)  and  Vl(i,j  ).  An  analysis  similar  to  the  one  in  the  last  subsection 
shows  that  T(j  ,  j  ,  r)  =  0(r  log  .r)  where  s  is  the  number  of  corners  from  v^  to  v 
counterclockwise.  Since  r  and  s  are  both  smaller  than  n,  this  means  the  /th  row  of  each 
matrix  can  be  filled  in  0{n  log  n)  time.  There  are  n  rows  in  each  matrix,  and  we  have  i-3 
sets  of  matrices,  so  Algorithm  3.2  runs  in  C>(tn-log  n)  time  and  requires  0{kn')  storage. 

As  we  mentioned  before,  the  cost  of  computing  optimal  chains  dominates  the  cost  of  the 
whole  problem,  so  we  conclude  that  the  minimum  area  enclosing  k-goa  can  be  found  in 
O(i:n-log  n)  time  by  using  0{kn')  storage. 

In  order  to  find  the  minimum  area  enclosing  *-gon,  optimal  chains  of  length  1  and 
optimal  flush  chains  of  lengths  from  1  to  t-3  are  computed.  But  are  all  these  chains 
absolutely  necessary?  In  the  next  section,  we  describe  an  improvement  that  does  without 
some  of  them.    Thus  less  time  and  storage  is  required. 

3.5.    Improvement 

Next  we  will  improve  the  0(*n-log  n)  time  to  0{n-\og  k  log  n)  by  applying  an 
observation  that  was  first  used  by  [Chang-Yap  (1984)].  For  all  h  between  1  and  k-2,  the 
above  procedure  computed  all  /j-sided  optimal  flush  chains  However,  we  now  show  that  to 
compute  the  matrices  V'.j  and  q*.  ,,  we  only  need  to  compute  V;,  and  a,  for  at  most  2  log  k 
values  of  h  To  show  this,  it  is  sufficient  to  prove  that  for  any  h  Si  2,  we  can  compute  the 
matrices 

in  0(n'log  n)  time  from  the  entries  of 

v:,v:.,.a;.a:.,  o) 

where  m  ((/i-l)/2J.  The  method  of  the  procedure  above  clearly  shows  that  V'^  and  al 
can  be  obtained  from  (1).  To  see  that  V^.,  and  a^.,  can  also  be  obtained  form  (1),  we 
observe  that 

(    1>./2J.|>,  2)   }  C   (m.m^n. 
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Wc  conclude  that 

Theorem  3.6.    For  an\  given  n-gon  P,  the  mmimum  area  enclosing  kgon  can  he  computed 
in  O(n-log  k  log  n)  time  and  O(n-log  k)  storage. 

3.6.    Unsolved  Problems 

The  minimum  area  enclosing  t-gon  problem  has  a  trivial  lower  bound  i\{n).  And 
indeed  in  the  case  *  =  3,  O'Rourke,  Aggarwal,  Maddila,  and  Baldwin  have  shown  a  linear 
algorithm  thus  closing  the  gap  between  the  upper  bounds  and  lower  bounds  [O'Rourke- 
Aggarwal-Maddila-Baidwin  (1984)).  But  for  the  genera!  case  i  >  3,  there  is  still  a  gap 
berween  the  upper  bounds  and  lower  bounds  and  the  optimality  of  our  algorithm  remains  to 
be  determined.  In  fact,  it  seems  that  to  achieve  a  subquadratic  algorithm  even  for  a 
minimum  area  quadrilateral,  some  new  geometric  properties  would  be  required. 


CHAPTER  4 


Minimum  Perimeter  Enclosing  Triangles 


4.1.  Statement  of  the  Problem 

In  this  chapter  wc  consider  the  same  smallest  i-gon  problem  as  in  Chapter  3,  except  that 
the  size  of  a  polygon  is  measured  by  its  perimeter  instead  of  area: 

Problem  4.1.    Given  a  convex  polygon  P,  find  a  Jt-gon  Q  which  encloses  P  and  has 
the  minimum  perimeter  among  all  A-gons  enclosing  P ,  k  ^  3. 

A  special  case  of  Problem  4.1  is: 

Problem  4.2.    Given  a  convex  polygon  P,  find  a  triangle  Q  which  encloses  P  and 
has  the  minimum  perimeter  among  all  triangles  enclosing  P. 

For     the     rest     of     the     chapter,     we     assume     that     P  =  (vq,  v, >■,_,),     and 

Q  =  (hv,  w..  .  .  .  ,H'^_i)  (the  corners  are  listed  counterclockwise).  Let  e^  denote  the  directed 
edge  of  P  from  v,  to  v,.,,  (0  s  i  <  n).  Similarly,  let  j,  denote  the  directed  edge  of  Q  from 
w,  to  w,.  |,  (0  s  /  <  it). 

In  this  chapter,  wc  show  a  property  of  the  solution  polygon  Q  in  the  general  case. 
Problem  4.1.  But,  we  only  give  a  finitcness  criterion  and  an  algorithm  for  the  special  case. 
Problem  4.2. 

4.2.  Criteria  for  the  Minimum  Perimeter  Enclosing  Polygon 

As  it  turns  out,  the  edges  in  the  minimum  perimeter  enclosing  *-gon  are  characterized 
by  a  more  complicated  condition  than  the  bisecting  property  for  the  area  case.  To  state  the 
finitcness  criterion,  we  need  the  following: 
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-A 


'•'2 


Figure  4.1. 

Notation.  Consider  an  edge  j,  of  Q.  (See  Figure  4.1.)  Let  <i>  be  the  interior  angle 
prccceding  s,  and  -y  the  angle  following  s,  in  a  counterclockwise  traversal  of  Q, 
<J>  =   •^H,_,H',w,_,, -y  =   ^H',.-,H',.  ]M',.    For  a  point  p  €  J,,  let  y5(j,,  p,  C)  denote 

/  cot  -^ 


g  cot 


where  f  =  py>\  and  g  =  p  v*, . ,. 

Dermition.    Consider  an  edge  5  of  C      From  the  previous  chapter,  we  know  s  passes  through 
either  one  or  two  consecutive  corners  of  P.    So  let  v^  and  v     be  the  corner(s)  that  s  passes 

through,  j   =  j  OT  j+  l(mod  n).    We  say  that  s  is 

(1)  high,  if  Z{s,  Vj,  Q)>  I, 

(2)  low,  if  ;i{s,  V  ,  Q)  <  1,  and 

(3)  balanced,  otherwise. 
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Figure  4.2. 


Lemma  4.1.  (Butterfly  Lemma  II)  Every  edge  in  a  minimum  perimeter  enclosing  k-gon  is 
balanced. 

Proof  Refer  to  Figure  4.2  for  discussion.  First,  consider  a  non-flush  edge  s  of  the  minimum 
enclosing  i-gon  Q.  Recall  the  definition  of  butterflies  from  Chapter  3.  Let  s  be  embedded 
in  a  butterfly  B  centered  at  v  (which  we  will  call  the  pivot  of  s)  with  two  tips  ab  and  a  b  . 
Let  L  be  the  length  of  the  edge  plus  the  length  of  the  part  of  the  wings  that  are  in  7",  where 
T  is  the  tangent  determined  by  s.  Without  loss  of  generality,  assume  L  =  ac  '  cc  ~  c  b' . 
For  the  optimalit)  of  Q,  L  must  be  at  a  minimum.  Let  us  perturb  s  counterclockwise  by  an 
angle  B,  resulting  in  a  new  edge  s    =  dd  .    Let 

/    =  dVj,     g    =  d  v^,     h  =  cd,     h    =  c  d  , 
^    =   ■i^v  da  =   <i>  -  b ,     y    =    ■^  b  d  V    =   y  +  b . 


Then 


,    _      /  sin  <|> 

sin  (4)-  5) 


g 


=      R  sin  II 


h  = 


J  sin  6 

sin  (4>  -  5) 

g  sin  5 


sin  (-y  ~  5)  sin  {y  ^  h] 

dL  =  (J  *g*h  -J-g-h) 


dy  =   b 


dy  6-0  6 
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[jjjj  /  (sin  (^  ->-  sin  6  -  sin  (4)  -  S))  _|_  g  (sin  >  -  sin  S-  sin  (-y  +  S)) 
i-o  6sin(<|>-6)  6  sin  (7+6) 

..      /  (sin  4)  -  sin  4)  cos  6  +  sin  8  +  sin  S  cos  ^) 

6-0  6  sin  (4)  -6) 

g  (sin  -y  -  sin  "y  cos  5  -  sin  5  -  sin  8  cos  y) 
5  sin  (7  +  6) 


=  lim 


sin  6 


6-0       6 


/(l-Hcos4>)       J?  (1  + cos  7) 
sin  4>  sin  7 


/  (1^  cos  4))       g  (l  +  cos  7) 
sin  4>  sin  7 


=  /  cot  ^  -  g  cot 


Since  C  is  an  edge  of  a  iBinimum  perimeter  enclosing  polygon,  =  0.    Hence 


f  cot 


g  cot 


1. 
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Figure  4.3.  edge  cc   embedded  in  an  ^4 -butterfly  centered  at  v, 

This  extremum  actually  has  a  meaningful  geometric  interpretation  (see  also  [Chakerian- 
Langc  (1971)]  ):  Construct  the  larger  circle  S  such  that  S  is  tangent  to  cc  ,  ah,  and  a  b  ,  ]i  B 
is  a  V-buttcrfly.  (If  B  is  an  A -butterfly,  then  5  is  the  smaller  circle,  having  the  same 
properties.)  Sec  Figure  4.3.  Lei  d  =^  SHah  and  d  =  SHa  b  .  Then  L  =  cc  -^  b  c  *  ca  = 
ad  +  b  d  .  Let  5*  be  the  unique  such  circle  which  contains  v..  Then  L  is  minimized  if  5  =  S' . 
since  any  other  line  through  w,  will  cut  5*  and  have  a  corresponding  larger  "tangent"  circle 
associated  with  it  f(ir  which  the  sum  ad  ~  b  d   is  larger. 

For  flush  edges  of  Q.  wc  argue  that  if  they  are  not  balanced  then  there  is  a  way  of 
perturbing  them  such  thai  the  perimeter  of  Q  is  decreased  The  argument  is  exactly  the  same 
as  the  one  given  in  Lemma  2.6 

We  conclude  thai  e\cry  edge  in  a  minimum  perimeter  enclosing  i-gon  is  balanced.    D 

In  addition  to  the  ahovc  lemma,  we  have  the  following  propcrts  for  the  solution  of  the 
special  case,  Problem  4.2: 

Lemma  4.2.  Givrn  a  pol\gon  P  and  a  trtanglr  T.  lei  Q  br  the  minimum  area  i perimeter i 
triangle  thai  encloses  P  and  is  congruent  to  T  Then  one  oj  the  three  edges  of  Q  is  flush 
with  P 


Proof      I  his  lemma   is  due   to  Aggarwal  and   DePamv     1  he   basic  uk-al   of  their   proof  is  the 
following:    Consider  perturbing  Q   by  rotating  each  of  its  edges  b\   some  angle      It  turns  mil 
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that  there  is  always  a  direction  of  rotation  for  which  the  area  of  Q  decreases,  contradicting 
the  optimality  of  Q.    See  [DcPano-Aggarwal  (1984)J  for  more  details.    D 

Lemma  4.3.  There  is  at  least  one  flush  edge  in  the  minimum  enclosing  triangle  Q. 

Proof.  Assume  for  the  purpose  of  contradiction  that  there  is  no  flush  edge  in  Q.  Then, 
according  to  Lemma  4.2,  we  can  rotate  the  triangle  Q  in  one  way  or  the  other  while  keeping 
the  same  intersection  points  with  P  and  maintaining  the  same  shape  so  that  the  area  of  Q  is 
decreased  (sec  also  the  shape  problem  in  Chapter  8).  Since  the  triangles  have  the  same 
shape,  a  decrease  in  area  also  results  in  a  decrease  in  perimeter,  contradicting  the  optimality 
of  g.    D 

In  summary,  every  edge  in  the  minimum  perimeter  A-gon  is  balanced.  For  the 
triangular  case,  we  have  the  additional  property  that  one  of  the  three  edges  is  flush.  Thzx 
amounts  to  a  finiteness  criterion  that  will  enable  us  to  develop  an  0{n^)  time  algorithm  for 
Problem  4.2,  provided  that  when  given  a  flush  edge  and  two  pivots,  we  know  how  to  find 
two  balanced  edges  through  the  pivots. 

In  the  next  section,  we  explore  more  geometric  properties  of  the  minimum  perimeter 
enclosing  triangles  which  make  an  O(n-)  time  algorithm  possible. 

4.3.    More  Geometric  Properties 

If  an  edge  j  of  C  is  flush  with  an  edge  e  of  P,  we  say  Q  is  j-anchorcd  or  ("-anchored. 
From  Lemma  4.3  in  the  previous  section,  we  know  that  we  can  restrict  our  search  for 
minimum  perimeter  triangles  to  ^-anchored  triangles  for  every  edge  e  of  P .  Furthermore, 
for  each  e  in  P,  we  only  have  to  consider  pivots  that  will  give  rise  to  balanced  tangents. 
Clearly,  there  are  0{n-)  candidates  for  such  pairs  of  pivots.  However,  we  will  show  in  this 
section  that  there  are  at  most  a  linear  number  of  such  pairs  and  that  they  can  be  found  m  just 
O(log-n)  "overhead"  time. 

For  the  rest  of  this  chapter,  we  fix  the  flush  edge  e  to  be  e^,  =  VqV,.  Let  A(v,)  be  the 
distance  from  v,  to  the  flush  edge  e.  If  ^  is  the  set  of  all  farthest  points  in  P  from  e,  then 
there  is  an  index  i  such  that  either 

e    =   {  V,  }  andO<  hiv^)  < <  h{y,)  >  /i(v,.,)  > >  /i(v„_,)  =  0, 

or 

e    =   ^.,^  and  0  <  /i(v,)  < <  /i(v,)  =  /i(v,.,)  > >  /ifv^.,)  =  0. 

In  the  following,  we  assume  e  =  {  v,)  for  simplicity.  It  is  not  difficult  to  modify  the 
arguments  for  the  other  case.  The  index  i  can  be  found  in  O(log  n)  time  b\  a  binary  search. 
Then  we  can  separate  the  edges  of  P  other  than  e  and  e    into  two  sequences, 
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IF,\1 


t-\]  =    [*|,  «2,    .   .   .   ,*,-,], 


G  =  [G,  I  1  S  I  ^  n-/]  =   [<„  «,. 


,e 


n-lJ- 


If  wc  consider  e  as  pointing  horizontally  to  the  right,  then  clearly  an  ^-anchored  triangle 
enclosing  P  must  have  its  right  side  pivoting  at  a  corner  v^  and  its  left  side  at  Vj,, 
\  <  a  <  r(mod  n),  t  <  b  -^  n-  l(mod  n).  For  all  the  ^-anchored  triangles  enclosing  P ,  we 
can  define  an  equivalence  relation  for  their  apexes  depending  on  the  r»o  pivots.  If  two 
triangles  have  the  same  left  pivot  and  right  pivot,  then  their  apexes  arc  equivalent.  We  can 
then  partition  the  plane  into  equivalence  classes  determined  by  this  relation.  All  the  triangles 
determined  by  the  sequence  of  tangents  (Eg.  •5',  7").  where 

s  €  H(<),  r  €  E(n-;)u{  £„_,,,  }-{£„.,}. 
have  their  apexes  in  the  region  bounded  by 

£,_,,E„£„.^  and£„.^_,. 
We  denote  this  region  by  Region{2i,  Ij).    See  Figure  4.4. 


I  \ 


/.-;-! 


-1-1 


Figure  4.4.  Region(2i,  2j) 


Dennition.  For  1  s  /  <  r  and  1  s  _;  s  n-t,  let  j,  and  .5,  be  the  respective  edges  overlapping 
F  and  G  in  a  triangle  determined  by  £p,  £,,  and  C  .  Let  M  be  an  r  x  /  matrix  over  the  set 
{  LL.  LH,  HL,  HH  ]  such  that  M(i,  j)  is 

(1)  LL,  if  J,  is  either  balanced  or  low  and  s.  is  either  balanced  or  high,  or  if  thc\   do  not 
intersect  each  other  on  the  half-plane  £q  , 

(2)  LH ,  if  s.  is  either  balanced  or  low  and  j,  is  low, 

(3)  HL,  if  j|  is  high  and  j-,  is  either  balanced  or  high,  and 
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(4)     HH ,  if  j|  is  high  and  j,  is  low. 
Wc  call  ht  the  ratio  matrix  for  F  and  G . 

Dennition.    If  M{i,j)  =  LL  and  there  exists  {k,  I),  k  >  i  and  /  >  j,  such  that  M{k,  I)  =  LL, 
then  wc  say  that  M(/,  j)  is  dominated  by  A/(J1:,  /). 


a  =  a 


Figure  4.5. 


Lemma  4.4.  Let  s^  =  be  arid  Sj  =  ca  be  the  respective  edges  overlapping  L  and  R  in  the 
triangle,  Aabc,  determined  by  £q,  L,  and  R,  tangents  of  P .  Similarly,  let  s,  =  be  and 
Jj  =  c  a  be  the  respective  edges  overlapping  L  and  R  ,  in  a  triangle.  Aa  b  c  ,  determined 
by  Eq,  L,  and  R  .    If  R  <  R    <  L(vnod  P)  and  p   =  s^Hs^  and  q  i  S2,  then 

(1)  Ji{s^,  p,  Aabc)  <  Ji{s\,  p  ,  ^ah  c  ). 

(2)  5(j;,  q,  Aabc)  >  >C(5,,  <?,  Aa  b  c). 

Proof.  Refer  to  Figure  4.5  for  discussion.  First  of  all,  since  R  <  R  <  L(mod  P),  c  lies 
between  a  and  c,  and  b  lies  between  a  and  b  .  Property  (2)  follows  from  (i)  cq  >  Fq ,  and 
fii)  -^  acb  <   ^a  c  b  . 

In  the  following,  we  show  that  Property  (1)  is  also  true.    Let  Ji{s^,p,  Aabc)  =   y,  and 
R{s^,p,Aabe)  =   h'.    Let 


^b,     p    =   ^b\     y  =  ^c,     y 
b  p  sin  3 


^c    (  =    ^acb), 


bp  cot  -^ 

^  2 


cp  cot 


1. 


.     bp 


sin  3 


cp 


-    c  P  ^in  "y 
sin  -y 


After  substituting  for  bp  and  cp ,  we  get 
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—                    .  cot  -^ 
_    b  p  sin  -y  sin  P 2_ 

JT  sin  y   sin  0  j-q;  X 
2 

Since  0<p<p<'TT,0<-y<7    <  -n ,  ac  >  ac  ,  and  ai)    >  at,  we  have 

0  3 

cot  -'^  cot  -'^ 

— ^<  ^. 


cot  -^f-         cot  ^ 
2  2 


and  by  the  law  of  sines 


sin  -y  sin  g     _    ab  ac  ^ 

sin  7   sin  P  ac  ab 


b  p  cot  ^ 
2 
y  <  =  >*' I     y  <  w. 

c  p  cot  -^ 

^  2 

Wc  have  shown 

Z{s^,  p,  t^abc)  <  R{s\,  p  ,  ^abc).a 


Lemma  4.5.  Each  column  of  M  can  be  characterized  as  a  sequence  of  LL  entries  followed 
by  a  sequence  of  either  LH  or  HL  er:tries  (possibly  empry)  which  in  turn  is  followed  by  a 
sequence  of  HH  entries 

Proof.    The  lemma  follows  from  the  following  two  properties  of  M: 

(1)  If  «(/,  ;-l)  is////  then  ^f{i,j)  =  HH. 

(2)  If  M(i,  ;^  1)  IS  LL  then  A/(/,  j)  =  LL  . 

Both  properties  (1)  and  (2)  are  the  consequences  of  a  symmetrical  version  of  Lemma  4.4.    D 

Lemma  4.6.  The  necessary  condition  for  Region(i,j)  to  contain  a  minimum  perimeter 
triangle  is  that  M (i ,  j)  €  LL  and  M(i^  1,  j+  !)€////. 

Proof.    Follows  from  Lemma  4.4.    D 
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Lemrni  4.7.  If  Regional,  j)  contains  a  local  minimum,  then  M(i,  j)  €  LL  and  M  (i ,  j)  is  not 
dominated  by  another  entry  that  is  in  LL . 

Proof.  If  M(i,j)  is  dominated,  then  M(i+l,j+l)  is  LL .  By  Lemma  4.6,  Region{i,  j)  can 
not  contain  a  minimum  perimeter  triangle.    □ 

Lemma  4.8.  There  are  at  most  0(n)  regions  that  contain  local  minima. 

Proof.  For  Region{i,  j)  to  contain  a  local  minimum,  M{i,j)  must  belong  to  LL  and  not  be 
dominated.  Otherwise,  by  Lemma  4.5,  M{i+  1,  j+  1)  €  LL.  We  now  show  by  induction  that 
there  are  at  most  r+l  undominated  LL  entries  in  an  rx  /  matrix.  Let  M{r,  i)  be  the  last  LL 
entry  in  the  rth  row  and  M(j,  I)  be  the  last  LL  entry  in  the  /th  column.  By  definition  all  the 
entries  that  fall  between  row  1  and  row  j-l  or  between  column  1  and  column  i-l  are 
dominated.  By  the  induction  hypothesis,  there  are  at  most  r- ;  + /-y  undominated  LZ,  entries 
that  fall  bcrween  row  /  and  row  r-l  and  between  column  y  and  column  /-I.  So  in  total,  the 
number  of  regions  that  arc  determined  by  a  pair  of  diagonal  LL/HH  entries  in  the  matrix  M  is 

i  +  j  +  r-i+l-j  =  l+r  =  0{n).  0 

The  amount  of  work  spent  on  finding  local  minima  in  a  region  is  independent  of  n,  the 
size  of  the  problem.    Thus 

Lemma  4.9.  It  takes  0(1)  time  to  determine  the  local  minima  in  a  given  region. 

4.4.    A  quadratic  Algorithm 

We  now  give  an  algorithm  which  exploits  all  the  properties  mentioned  above: 

Algorithm  4.1.  Minimum  Perimeter  Triangles 

Step  1       Let  e  be  e^,  of  P . 

Step  2       Divide  the  rest  of  P  into  two  sequences  of  tangents  F  =  {  F,  |  /  =   1 r  }  and 

C  =  {  G,  I  /■  =    1,  ...,/)  Let  M   be  the  ratio  matrix  for  F  and  G.    Let  r,„„   =    1, 

Step  3        Let  I  =    (r,^,  -  r^,)  div  2. 

Step  4       Perform  two  binary  searches  to  find  the  last  LL  entry  M{i,a)   and  the  first  HH 
entry  M(i,b)  m  the  i-th  row  of  M . 

Step  5       For    ;  -  /^,^ ,  .  .  .  ,a,     if     M((*l,y+1)  €  HH     then     compute     the     minima     in 

Regior\(i.  j)  and  register  the  minima. 

Step  6       If  /  <  r^_  and  a  ^  /^^^ ,  let  r,„„   =  (+1  and  /^,  =  a-^1.    Then  recursively  do  Steps 
3-5. 
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Step  7      If  /  >  r,„^,  and  b  <  />„-l.   let   r^,  =  i-l    and   /,„,   =6-1.     Then   recursively   do 
Steps  3-5. 

Steps       For?  =  «,,*2 e„_^,  repeal  Steps  2-7.    o 

Note  that  Steps  5&6  arc  justified  by  the  fact  that  all  the  entries  in 
A^ftoH-'""  1.  '■/<,»  •■''"21'  are  in  LL  \J  LH  and  those  in  M[i-  l../^,,  a  *  2.  r^,]  are  in  HL  UHH  ,' 
all  of  which  can  be  ignored,  since  they  are  not  capable  of  producing  diagonal  LL/HH  pairs. 

4.5.  Analysis 

Let  T{n,  m)  be  the  number  of  steps  needed  to  seek  out  all  the  regions  that  can  possibly 
contain  local  minima  for  all  triangles  anchored  on  an  edge.  Since  Step  4  takes  only 
logarithmic  time,  we  have 

7(n,  m)  =  log  m  +  r(y,  m,)  +  7(y,  m-m^-^2). 

That  is 

T(n,  m)  =  O(log-n). 
So  the  total  number  of  steps  needed  is 

rt  X  max(0(log-n),  0(n))  =   0(n=). 
Wc  have  proved 

Theorem  4.10.    Minimum  perimeter   triangles  enclosing  a  given  convex  polygon  P   can  be 
found  in  0{n-)  time 

4.6.  Unsolved  Problems 

The  perimeter  versions  of  the  enclosure  problems  are  obviousl>  more  difficult  than  the 
area  versions.  In  this  chapter,  we  solve  the  triangular  enclosure  problem  with  the  perimeter 
measure.  We  do  not  know  whether  our  algorithm  is  optimal,  since  we  do  not  have  any  lower 
bound  for  the  triangular  case  or  the  general  case  except  for  the  trivial  linear  lower  bound 
Furthermore,  we  do  not  know  whether  this  method  applies  to  the  general  case  where  the 
enclosing  polygons  have  more  than  3  sides: 

Problem  4.3.    Given  a  convex  polygon  P  and  a  fixed  number  i  >  3,  find  a  *-gon 
enclosing  P  with  minimum  perimeter  among  all  *gons  enclosing  P . 


^  M[a    b .  c.  d]  denotes  the  matri;i  that  is  obtained  by  taking  the  slice  of  Si  from  row  a  to  row 
b  and  column  c  \o  column  d. 

^  This  follows  from  Lemma  4.5. 
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And  also  unsolved  is  the  generalization  from  2  dimensions  to  3  dimensions: 

Problem  4.4.     Given    a    convex    polytope    P,    find    tetrahcdrals    containing    P    with 
minimum  surface  area  among  all  tetrahcdrals  that  contain  P . 

Without  doubt,  the  problem  becomes  even  more  difficult  in  higher  dimensions. 


CHAPTER  5 


The  Balanced  Chain  Problem 


In  Chapter  2,  wc  dealt  with  the  smallest  corner  triangle  problem  which  asks  for  a  line 
segment  such  that  (i)  its  two  end  points  are  on  the  two  given  lines  respectively,  (ii)  the  three 
lines  form  a  triangle  that  encloses  a  given  set  of  points,  and  (iii)  the  area  of  the  triangle  is 
minimized.  Now,  we  generalize  the  problem  by  asking  for  a  convex  polygonal  path 
originating  from  one  of  the  given  lines  and  ending  at  the  other  line.  We  are  not  only 
interested  in  the  problem  itself,  but  also  in  its  implication  for  another  problem  to  be 
discussed  in  the  next  chapter,  namely  the  potato-peeling  problem. 

5.1.    The  Smtliest  Corner  Problem 

The  following  definitions  are  needed  for  the  statement  of  the  problem: 

Dennition.  For  a  directed  line  L  ,  let  Z. "  denote  the  closed  half-plane  to  the  left  of  L ,  and  L  '  , 
the  closed  half-plane  to  the  right  of  L.  For  a  pair  of  directed  lines  L  and  M,  if  Z.n,V/'  is 
directed'  toward  LDM ,  then  we  say  that  L  precedes  M .  See  Figure  5.1 


Figure  5.1.  Z.  precedes  M 


Here,  we  assume  that  L  O M     has  the  same  direction  as  L . 
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Definition.  Let  ^  be  a  set  of  points  in  the  plane.  For  a  pair  of  directed  lines  (Lq,  Z.,)  such 
that  Lq  precedes  L,,  let  an  ({Lq,  Z,,,  R))-chord  (or  just  chord  for  short)  denote  a  directed  line 
segment  contained  in  Lq  DZ,,"  passing  through  at  least  one  point  of  R  and  with  its  two 
cndpoints  in  Lq  and  Z,,  respectively.  The  direction  of  a  chord  is  from  the  endpoint  in  L^  to 
the  one  in  L,.  The  points  of  /?  in  a  chord  are  called  the  pivots  of  the  chord.  We  say  that  a 
sequence  of  (Z-q,  Z,,,  /?)-chords  C  =  (C,,  Cj C,„)  is  upward  convex,  if 

(1)  For  1  =  0 m,  C ^  intersects  C,.  ,  at  x,  (Cq  =  Lq  and  C„.  ,  =  Z.,). 

(2)  For  I  =  0,  .  .  .  ,  m  -  1,  C,  precedes  C,^ ,. 
Let  p(C)  denote  the  simple  polygon  of 


Lo-nzLf  n(nc,-) 


1- 1 


A  sequence  C  of  (Z,q,  Z,,,  /?)  chords  is  said  to  enclose  R  if  every  point  of  R  which  is  in 
Lq  r\L^ ,  is  also  in  p(C).  See  Figure  5.2.  If  a  sequence  C  =  (^€^,€2,  ■  ■  ■  ,C„)  of 
(Lq,  Z,,,  /J)-chord  is  upward  convex  and  encloses  /? ,  then  we  say  C  is  an  (Lq,  L^,  R)-chain  of 
length  m.  We  call  it  an  (Z-q,  L^)-chain  when  Z?  is  understood,  or  simply  a  chain  when  (Lq,  Z,,) 
and  R  arc  understood.  x,_|.x,  is  an  edge  of  C  and  also  known  as  the  truncated  version  of  the 
chord  C,  («'  =  I.  .  .  .  ,m)  and  j,  is  a  node  of  the  chain.  Z,q  and  L,  are  called  the  backward  and 
forward  supporting  lines  respectively. 


LA  =  C, 


^oi=Co) 


Figure  5.2.  An  (Z.^,,  Z.,,  /?)-chain,  (C,,  C,,  C3); 
the  points  of  R  are  indicated  by  'x  ' 


Now  we  can  state  the  smallest  corner  problem  formally: 

Problem  5.1.  Given  a  pair  of  directed  lines  (Z.^,  Z.,)  in  the  plane  {Lq  precedes  Z.,) 
and  a  set  of  points  R  in  Z.q'  HZ.,"  ,  find  an  (Lq,  Z.,,  Z;)-chain  C  such  that  the  area  of 
p(C)  is  minimized. 

For  the  rest  of  the  chapter,  wc  fix  an  arbitrary  pair  of  directed  lines  (Lq,  Z.,)  and  a  set  R 
of  n  points  R  in  Z.q  nz,,"  .  For  simplicity,  assume  that  there  are  no  three  coUinear  points  in 
R. 
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5.2.    The  Balanced  Chain  Problem 

The  solution  to  the  smallest  corner  problem  exhibits  a  balancing  condition  similar  to 
that  of  the  simplified  case  described  in  Chapter  2.  So,  first  we  generalize  the  notion  of 
balancing  in  order  to  characterize  the  solution. 

Notation.  For  two  distinct  points  p  and  <?,  let  /?  |(7  denote  the  reflection  of  p  about  ^,  i.e.  p  |^ 
is  a  point  on  the  line  pq  such  that  q  is  the  midpoint  between  p  and  p  \q . 

Definition.       Let      C     be     an      (Lq,  L^,  R) -chain     with      nodes      (xq- ^i ^m)      ^"^      'c' 

•)t,_,x,n/?  =  {  p^,  q^  }  (possibly  p,  =  q,  and  p,  is  between  jt,. ,  and  q^.  Sec  Figure  5.3.  The 
edge  x,_,i:,  is  said  to  be  area  balanced  (or  just  balanced  when  the  metric  used  is  understood 
to  be  the  area),  if  the  four  points 

X,.,,    x\q^,    x,.,|p,,    X, 

occur  in  that  order.    Note  that  the  four  points  need  not  to  be  distinct.    We  allow  the  cases 

(1)  x,|<7,  =  x,_|Lp,,  and 

(2)  P,  =  q,  (so,  X,.,  =  xj^,  and  x,_,|p,  =  x,). 

In  case  (2),  we  have  p^,_|  =  p^,  which  is  our  original  notion  of  balancing  in  Chapter  2. 


9, 

-X— 


P, 

-X— 


>, 


-f,W, 


Figure  5.3.  A  balanced  double-pivot  edge  x,_|X, 


If  p,  ^  q  ,  we  call  {p  ,  q  )  a  double-pivot .  In  this  case,  C,  is  double-pivot  chord  and  x,_  |X,  is  a 
double-pivot  edge  Otherwise,  p^  =  q  and  we  call  p,  a  single-pivot.  Similarly,  C,  is  called  a 
single-pivot  chord  and  x  _  |X,  is  called  a  single-pivot  edge.  A  chain  C  is  balanced  if  every  edge 
in  L  is  balanced  and  a  simple  chain  is  a  chain  that  consists  of  singlc-pivot  edges  only.  Thus,  a 
simple  balanced  chain  is  a  chain  that  consists  of  single-pivot  balanced  edges. 


Let  c_  be  the  center  of  edge  x,.|X,  and  let  x,  and  x,_|  denote  x  ,|p^  and  xj^, 
respectively.  It  is  clear  that  if  x,.|X,  is  balanced  then  r,  is  between  p.  and  q  Otherwise 
x,.|X,  is  not  balanced  and  one  of  the  following  three  conditions  is  true: 


(i)      p,  and  q^  arc  in  x,_,r,, 
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(ii)     p,  and  q,  arc  in  c,x,; 

(iii)    x,_,|p,  appears  before  x}q,  in  a  traversal  of  C  from  jtq  to  x^. 


We  call  -t,_|i:,  an  early  edge,  a  /aff  edge,  or  a  splittable  edge  according  to  which  of  (i),  (ii), 
and  (iii)  is  true.    Sec  Figure  5.4  . 


(a)  an  early  edge 


(b)  a  late  edge 


P,     9. 
~)C  ■><    - 


S-l 


9,     P, 

-X-K  - 


S-1 


(c)  a  spittable  edge 


Figure  5.4. 


Lemma  5.1.  An  {Lq,  L  ^,  R ) -chain  C  has  the  minimum  area  among  all  (Z-q,  Z.  |  ,  R)-chains,  if 
it  is  balanced. 


Proof.  Let  C  be  an  {L^,,  Z.,,  /?)-chaiD  with  minimum  area.  Let  x,_|jr,  be  an  edge  of  C  with  the 
double-pivot  (p,,  <?,)  and  assume  for  the  sake  of  contradiction  that  x,_,x,  is  not  balanced. 
First,  we  consider  the  case  that  x,_,j:,  is  an  early  (resp.  late)  edge.  (Refer  to  Figure  5.4  for 
discussion.)  It  is  easy  to  see  that  by  perturbing  x,_,x,  infinitesimally  clockwise  around  q^ 
(resp  counterclockwise  around  p,),  we  can  decrease  the  area  of  J,  contradicting  the 
optimality  of  t.  (Sec  Figure  5.4(a)  and  (b).)  In  the  case  of  a  splittable  edge,  we  can  break 
the  edge  x,_  |X,  at  the  midpoint  c  between  J:,_  Jp,  and  x}q^  into  two  edges  x,.  ^c  and  rx,  (thus 
the  term  "splittable").  It  is  clear  that  x,_|C  is  an  early  edge  and  ex,  is  a  late  edge.  Then  by 
an  infinitesimal  clockwise  perturbation  of  x,_,r  around  p,  and  a  countcrclock\i  isc 
perturbation  of  ex,  around  q  .  we  can  decrease  the  area  of  p(C),  a  contradiction.  (See  Figure 
5.4(c).)  Note  that  after  the  perturbations  mentioned  above,  the  chain  C  remains  upward 
convex  and  still  encloses  R.    We  conclude  that  in  all  cases  every  edge  of  a  minimum  area 
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chain  C  must  be  balanced.    □ 

By  Lemma  5.1,  the  smallest  corner  problem  reduces  to  the  following  problem  of  finding 
balanced  chains: 

Problem  5.2.    For  a  pair  of  directed  lines  (Lq,  L,)  in  the  plane  (Lq  precedes  L^), 
and  let  ^  be  a  set  of  points  in  Lq'  flLf  .    Fmd  balanced  (Lq,  L^,  ^)-chains. 

Balanced  chains  have  some  interesting  properties  which  arc  also  useful  for  an  algorithm. 
In  the  two  following  sections  we  show  these  properties.  An  O(n')  time  algorithm  will  be 
given  in  the  last  section. 

5.3.    The  Existence  and  Uniqueness  of  Balanced  Chains 

Balanced  chains  exhibit  some  very  surprising  properties.  It  turns  out  that  of  all  the 
upward  convex  chains  supported  by  a  pair  of  directed  lines  that  enclose  a  given  set  of  planar 
points,  there  is  one  and  only  one  balanced  chain.  In  this  section  we  give  a  proof  for  the 
existence  and  uniqueness  of  the  balanced  chain.  The  proof  itself  is  interesting  and  gives 
some  intuitive  idea  of  the  difficulties  involved  in  finding  such  a  balanced  chain  and  hints  at 
the  process  of  decomposition  of  simple  balanced  chains  to  be  discussed  in  the  next  section. 
Furthermore,  the  existence  proof  is  a  constructive  one;  it  is  an  algorithm  for  finding  balanced 
chains,  although  we  do  not  have  an  upper  bound  for  the  algorithm. 

Theorem  5.2.    Let  (Lp,  L^)  be  a  pair  of  directed  lines  and  R  be  a  set  of  points  on  the  plane. 
Then  there  is  one  and  only  one  balanced  (Z.q,  L^-chain  with  respect  to  R. 

We  show  a  useful  lemma  along  the  way  to  proving  this  theorem.  The  following 
notations  arc  needed  for  this  lemma: 

Notation.    Let  C  be  a  balanced  (Z,p,  Z.,,  ^)-chain  and  let  p    and  <?,  be  the  pivots  of  the  chord 

C,  in  C  (possibly  p,  =  q^)  and  let  (xq,  x, j„)  denote  the  nodes  of  C.    If  p,  *  q    then 

assume  q^  lies  berween  p ^  and  x,.  Recall  that  for  any  point  p  ,  L^{p)  denotes  tfic  line  through 
p  and  parallel  to  Z.,.  Let  5,(C)  denote  the  strip  berween  the  parallel  lines  L|(x,_|)  and  L^{x) 
where  Z.|(x,  _ ,)  is  excluded  from  S  but  i.|(x,)  is  included.  Let  /1,(l)  denote  the  strip  berween 
the  parallel  lines  Z.|(x,_|)  and  L^(p).  It  is  important  to  note  that  we  exclude  the  line  L|(j,_|) 
from  A^  but  include  Z.,(p,)  in  A,.  Let  fl,(J)  denote  the  strip  berween  Z.,(p,)  and  Z.|(x,_  ||p,); 
again  Z.,(p,)  is  excluded  but  L  {x ^  ^\p )  is  included.  Similarly,  D,(C)  (resp.  £,(l))  denote  the 
strip  and  between  the  lines  L^(x\q)  and  L^(q)  (resp.  /.,(<?)  and  Z.,(j:,))  where  L^{x\q)  (resp. 
L^(q))  is  excluded  but  /.,(£?,)  (resp  Z.|(.r,))  is  included  \N  e  will  just  say  5,,  /t,,  fl,,  C„  and  D 
for  shon,  when  the  chain  in  question  is  understood     .Note  that  ,4,nfi    =  0,  D,C\E,  =  0  and 
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5,iC)  =  A,UB,UD,UE,. 

Definition.  Recall  that  the  chain  C  divides  the  quadrant  Lq"^  nz.*  into  a  finite  region  p(C)  and 
an  infinite  region.  It  is  convenient  to  consider  p(t)  to  be  above  the  chain  and  the  infinite 
region  to  be  below  it.  Thus,  each  of  the  strips  5,,  A,,  fl,,  etc.,  will  be  broken  up  into  two 
half-strips ,  above  and  below  the  chain,  respectively.  See  Figure  5.5  and  Figure  5.6  for  these 
regions. 


E     ,    D. 


B. 


ii(',)  ^,0',)  ^,(x,-,) 


'   P,  =  ?, 

B^  =  E, 

A,=  D, 

V-l 


Figure  5.5.  The  regions  A,,  D,,  fl,  and  £,,  two  cases 

Notation.  Let  o,(C)  =  A,UD,  and  P,(^)  ~  B,^^r  Note  that,  except  in  the  extreme  case  of 
'i-\\p,  ~  ^,k,.  "i^P;  is  non-empty  because  of  the  bracketing  property.  Also  p,  =  q^  if  and 
only  if  A,  =  D,  if  and  only  if  B,  =  E,. 


B,nD, 


Figure  5.6.  The  shaded  areas  are  forbidden  for  nodes  of  C  ,  two  cases 
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Lemma  5.3.  Let  C  be  a  balanced  (Lq,  Z.,,  R)-chain  and  let  L^  be  a  directed  line  preceding 

i.|.  Lq  *  Lq.    Let  C    =   (Cp  Cj C„  )  a  balanced  (L^,  L,,  R)chain  and  let  x/j.  p^s 

be  the  nodes  and  pivots  of  C  ■    Then  x    can  not  lie  above  C  in  0|(C)  for  any  i  and  j. 

Proof.  Assume  for  the  sake  of  contradiction  that  Xj  lie  above  C  in  a,  for  some  ;  and  i.  If 
j  =  m  ,  by  definition  x„  lies  in  £,„  and  hence  is  not  in  any  a,.  Choose  j  to  be  the  largest 
index  such  that  x  (J  <  m  )  is  above  the  chain  in  o,  for  some  i.  Hence  x  ,  and  x  j  arc 
either  below  C  or  not  in  a^  for  some  *.  We  will  prove  that  x^,,  or  x^.j  ''^^  above  the  chain 
in  some  a^  (k  >  /),  thereby  contradicting  our  choice  of  j.    There  are  two  cases  to  consider. 

(a)  Consider  the  case  where  x  ,  is  in  5^,  for  some  k  >  i.  If  x^. ,  lies  below  C ,  then  p  , 
would  be  below  C-  This  contradicts  the  assumption  that  there  arc  no  points  below  C.  Thus, 
Xj,  I  lies  above  C  in  5^-  a^. 


L,{x)     L,{j>:)    Z.,(x,_,) 


^iK) 


Figure  5.7. 


Suppose  X|.|  lies  above  the  chain  in  fc",  (the  proof  for  B,-a^  is  similar).  C,.,  must 
intersect  L^(q^)  below  the  chain  (otherwise  q,  would  be  in  ^tiC  ),  contradicting  the  emptiness 
of  C  )  Let  r  denote  this  intersection  point  Similarly  C,. ,  intersects  Z. ,(;?,)  below  C  at  some 
point  y  (sec  Figure  5.7).  Then  the  midpoint  of  the  segment  x  .  x,.  ,  lies  in  v7  Note  that 
Pj.  I  must  lie  to  the  left  of  >•  and  i?, .  ,  to  the  right  of  z  to  satisfy  the  bracketing  property;  this 
implies  both  p^.  ,  and  q^ .  ,  lie  below  the  chain  l  ,  contradiction 

(b)  Consider  the  case  when  x  ,  stays  in  5,.  If  x  . ,  is  below  C.  then  an  argument 
similar  to  part  (a)  shows  that  x  . ,  is  above  '  in  a,  for  some  k  •  i  Otherwise  x  .  ,  is  above  J 
in  3,  and  we  have  3  possibilities: 

(1)     x^  (  -4,  andx^.,   (   0  , 
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(2)  x]  i  D,  and  Xj.^  €  £,. 

(3)  both  J    and  Xy.  I  arc  in  fl, no,. 

Now,  (1)  implies  p,  ^  p(C  )  and  (2)  implies  q,  q  p{C  )  contradicting  the  fact  that  C  encloses 
R.  If  (3)  holds,  then  this  contradicts  our  choice  of  j.  So,  in  all  cases,  we  have  a 
contradiction.    D 


^i(-«o1Pi) 


L,{p,) 


^i(-»;o)       ^c 


Figure  5.8. 


Proof  of  Theorem  5.2    There  are  two  parts  to  this  theorem:  (i)  (uniqueness)  there  cannot  be 
more  than  one  balanced  chain  and  (ii)  (existence)  there  exists  at  least  one  balanced  chain. 


(Uniqueness)      Let     C  =   (C,,  C;. 


,  C.)     be     a     balanced     chain        V>  e     derive 


contradiction  by  assuming  the  existence  of  another  balanced  chain  l     =   (C,,  C; C„  ). 

Let  p/s,  jt/s  and  p/s,  j:,'s  be  the  pivots  and  nodes  of  C  and  C  respectively.  To  apply  the 
previous  lemma  we  only  have  to  show  that  x  lies  in  a..,  for  some  i.j  Initially  assume 
*o  *  Xq.  Suppose  Jq  is  below  Xf,  on  L^,.  (See  Figure  5.8.)  Observe  that  i,  must  lie  above  C, 
otherwise  p ,  is  below  C  ■  U  x^  lies  in  any  a,  wc  are  done.  So  let  x ,  lie  in  P  -  a ,.  Note  that 
X,  is  to  the  right  of  Z.,(xo1p  |).  If  x,  is  in  fi,  then  p,  is  below  i ,  contradiction.  If  x,  is  in  f ,, 
let  7  be  the  intersection  of  C,  with  L^{q^).  Then  :  must  be  below  J.  Note  that  the  midpoint 
of  C,  hcs  to  left  of  z  and  hence  p,  lies  to  the  left  of  z.  This  implies  p,  is  below  J, 
contradiction.  Therefore  x^  cannot  lie  below  x  ,.  If  x^  is  above  X(,  on  L,,  then  a  symmetrical 
argument  applies  by  exchanging  the  roles  of  C  and  C  ■ 

It  remains  to  consider  the  possibility  x,,  -  x^^.  Since  C  *  J  .  let  x  .  ,  be  the  first  node 
such  that  X.,  *  x.,.  Suppose  x^.  ,  does  not  lie  in  the  line  through  xTT.,.  An  analysis 
similar  to  the  above  shows  that  x     ,  lies  in  q,  for  some  k.    Finally,  suppose  x,.,  lies  in  the 
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line  through  xJT^j.^.  If  it  lies  in  the  segment  i^7~*,-i  then  x^.j  must  lie  in  some  a^  by  the 
same  argument  as  the  previous  case.  Otherwise,  x^.,  clearly  is  in  A^.j  ^^'^  again  we  have  a 
contradiction. 

(Existence)  The  existence  proof  can  be  regarded  as  an  algorithm  although  we  do  not 
know  of  a  polynomial  time  bound  on  its  complexity. 

We  now  give  a  'scan  line'  algorithm  for  computing  the  balanced  (Z.q,  L|)-chain.  A 
sequence  of  (LqC')-  Z.,)-chains  are  computed  where  LqC')  '*  ^  ^^^^  parallel  to  Lq  at  a  distance  i 
to  the  left  of  Lq.  Imagine  the  line  Z-qCO  moving  from  the  infinite  right  toward  Lq  as  t 
approach  0  from  ».  As  r  -  0,  Lq(i)  -  Lq  and  the  (Z,o(r),  Z.,)-chain  becomes  the  {L^,  Z.,)-chain. 
During  the  process,  there  are  events,  t^,  ;,,...,;„  that  divide  the  scanning  process  into 
intervals  where  changes  in  the  (/.o(0.  ^i)-chain  are  smooth  within  each  interval.  More 
specifically,  between  events  the  pivots  remain  the  same  and  the  slope  of  each  chord  changes 
at  a  smooth  rate  (with  respect  to  t).  Initially  the  (Lq{'*>),  L,)-chain  consists  of  an  infinite 
chord  parallel  to  L^  through  the  point  p,  in /fnz.,"  that  is  farthest  from  L,.  Then  this  chord 
turns  continuously  clockwise  about  p,  as  Lq(i)  gradually  moves  toward  Lq  until  the  chord  hits 
a  new  point  q,  in  R.  This  is  the  first  critical  moment  /q  and  it  has  a  corresponding  balanced 
chain,  C  =  (xqX,). 

As  L^it)  continues  to  move  closer  to  Lq  and  Xq  moves  toward  x,  (as  usual,  x,  denotes  the 
i-th  node  of  the  balanced  chain,  C,  is  the  i-th  chord,  etc.),  Xq|p,  moves  toward  both  xj^,  and 
Pi  and  the  midpoint  r,  of  C,  moves  toward  q^.  Two  things  can  happen  at  the  next  critical 
moment  i  =  r,:  Either  (a)  XqIp,  meets  xjij,  and  the  double-pivot  chord  'splits'  into  two 
single-pivot  chords  balanced  at  p^  and  q^,  or  else  (b)  r,  meets  <?,  and  the  chord  starts  turning 
clockwise  around  q^  'leaving'  p,  behind.  If  (a),  the  splitting  event  occurs  first,  we  will  have  a 
chain  of  size  2.  And  as  Z.q(0  moves  further  in,  in  order  to  balance  C,  and  C;,  Xq  and  x, 
move  in  the  direction  of  /. ,  and  x,  moves  in  the  opposite  direction  of  /.,.  This  continues 
until  t  '-  r,  when  a  point  in  R  is  hit  by  either  C,  or  C,  and  turns  a  single-pivot  chord  into  a 
double-pivot  chord.  If  (b),  the  leaving  event  occurs  first,  the  double-pivot  chord  turns  into  a 
singlc-pivot  chord  while  continuing  to  turn  counterclockwise  unul  t  =  j,  when  it  hits  another 
point  of  R  and  turns  itself  into  a  double-pivot  chord.  So  this  process  of  forming  double- 
pivots,  splitting,  and  leaving,  continues  until  i  =  t^  when  L^it)  reaches  L^. 

In    general,    consider    the    (L^(t),  L|)-chain,    (C,,  C, C„).     Let   C      (for   some  j^, 

\  ^  Jq  s.  m-t-1)    denote    the    rightmost    double-pivot   chord.     We   choose  j^   to    be    m-1    and 
C     -  i.|,  when  there  arc  no  double-pivot  chords     As  L^it)  moves,  the  even-numbered  nodes 

(j,,,  2i  <  j^)   move    in   the   direction  of  C   ,   while   odd-numbered   nodes   (j:-,.,.  2i-1  <  Jq) 

move  in  the  opposite  direction  of  C   .    Observe  that  C  ,  for  i  =    1 ;.  -  1,  turns  clockwise 

if  I  is  even  and  counterclockwise  otherwise.    This  implies  that  the  length  of  the  truncated 
version  of  C     decreases,  if  y^,  is  even.    The  rest  of  the  (L^(t),  Z.|)-chain  remains  unchanged. 
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Wc  can  classify  the  possible  events  into  five  categories: 

(1)  (flanening)  Two  consecutive  chords  C-,,  and  Cj.-i  could  flanen  out  and  become  a 
double-pivot  chord.    It  should  be  noted  that  C,,-,  and  C2,  cannot  flatten  out 

(2)  (hitting)  One  of  the  single-pivot  chord  could  hit  a  point  in  /? ,  thus  becomes  a  double- 
pivot  chord. 

(3)  (splitting)    If  y'o  '*  odd,  the  truncated  version  of  C^   could  shorten  to  such  an  extent  that 

X,     ,|p     meets  x,  \q ,  .    Then  C,    loses  the  bracketing  property  and  splits  into  two  single- 

H         ■'0  .'0    •'0  -'0 

pivot  chords. 

(4)  (leaving  p   )    When  C.    is  shortening,  another  situation  could  also  arise.    The  midpoint 

r^   could  meet  q^    before  splitting  occurs.    Subsequently,  Cj   will  turn  counterclockwise 

around  a    leaving  p,    behind.    In  effect,  C,    turns  into  a  single-pivot  chord  balanced  at 
Vo  h  •'0 

(5)  (leaving  q^)  Finally,  if  ;o  is  even,  the  chord  C^  is  lengthening  and  the  only  event  that 
can  happen  is  the  midpoint  r  moving  right  and  meeting  Pj  .  Subsequently,  Cj  will 
pivot  clockwise  around  p.  ,  leaving  q.   behind. 

The  algorithm  first  establishes  a  balanced  chain  for  (Z,(rg),  L,)   and   then  repeats  the 
following  steps  until  Z,o(')  reaches  L^: 

Step  1       Find  the  rightmost  double-pivot  chord,  C^  ,  on  the  current  (i.Q(f),  L,)-chain.    In  the 
case  of  a  chain  consisting  of  just  single-pivot  chords,  take  i.,  for  Cj  . 

Step  2       For  each  chord  C,,  compute  the  value  of  t  when  an  event  involving  C   occurs: 

(1)  It  is  not  too  hard  to  compute  the  values  of  t  when  flattening,  leaving  or 
splitting  occurs  at  each  C  ,  since  those  are  determined  b)'  the  chain  alone. 

(2)  The  hitting  event  at  each  C,  is  more  difficult  to  compute  since  it  involves 
points  not  on  the  chain  But  it  is  clear  that  the  point  hit  by  C  has  to  be  on 
the  convex  hull  of  some  subset  of  R  containing  p..    So  we  partition  the  plane 

into  7o  strips  and  two  half-planes  by  the  lines  L^  (x,),  i  =   0 ^q.  where 

L  [x)  is  the  line  through  x  parallel  to  C  Let  H  be  the  con^■ex  hui!  of  those 
points  of  R  in  the  strip  between  /.,(x,_|)  and  Z,,(i,).  The  point  hit  by  an 
even-numbered  chord  C^,  is  the  next  point  on  //,,  clockwise  from  p-,  ;  the 
point  hit  by  C-,,. ,  is  the  next  point  on  //-  .  ,  counterclockwise  from  P;  -  i- 

Step  3       Find  the  smallest  value  for  i,  which  determines  the  next  event. 
Step  4       Update  (Z.p(/),  Z.|)-chain  accordingly. 


58 


Wc  conclude  that  there  exists  a  balanced  (Z.q,  L,,  /?)-chain.    D 

5.4.    The  Decomposition  of  Simple  Balanced  Chains 

From  the  proof  of  the  last  section,  wc  get  some  idea  of  the  difficulty  involved  in 
computing  the  balanced  chains.  We  know  of  no  simple  way  of  computing  balanced  chains 
directly.  Essentially,  wc  will  break  the  problem  up  into  many  smaller  problems,  find  the 
solution  for  each  of  them,  and  finally  merge  these  partial  solutions  to  form  the  solution  for 
the  original  problem. 

Consider    a    balanced    (Lq,  Z.,,  ^)-chain,    (€^,€2 CJ.     The    middle    chord    C„-, 

definitely  divides  the  chain  into  two  smaller  balanced  {Lq,  C„  2.  ^  )"  and  (C^-,,  L,,  R  )-chaiDS 
where  R  and  R  are  subsets  of  R.  However,  for  an  algorithm  to  work,  C^^  must  belong  to 
some  a  priori  finite  set  of  chords.  But,  in  general,  such  a  condition  is  not  met,  except  when 
C  ^  is  a  double-pivot  chord.  Our  idea  is  to  transform  a  simple  balanced  chain  in  some  way 
so  that  one  of  chords  becomes  a  double-pivot  chord.  Before  showing  how  that  is  done,  wc 
introduce  some  notation. 

Notation.     Let  R    be   a  set  of  n    points,   and  {pp/'j P„}    C   l^      For   a   line   L,   an 

(L,  P|,  /J;.  .  .  .  ,p^-path  is  a  polygonal  path  of  the  form 

n    =    (-To.-r, -^m), 


where  x^  i  L  and  for  i  =   1 m,  x,.  ,x,  passes  through  p^  and  the  center  of  x,.  |X,  is  p ,. 


L{x,) 


Figure  5.9.  An  (Z.,P|,p;,p,)-path 


For  an)    pom:  x   and  a  line  /. ,   let  L(x)  denote   the   line  which  passes  through  x   and 
parallel  to  L     See  Figure  5  9.    Note  that  the  line  L(i  )  is  a  function  of  L  and  ;> , ,  p,.  .p 
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only    (i.e.    L(x)    is     independent    of    />,.,,  .  .  -  ,p„    and     the     particular    choice    of    11). 

Furthermore,     for     any     /  =  0 m,     any     point    x     on     Z,(x,)     determines     a     unique 

(L,p^,P2 p„)-path  and  vice  versa.    In  particular,  the  choice  of  any  x  i  L  determines  a 

path.  Wc  can  think  of  a  path  as  a  configuration  in  a  system  of  "interconnecting  levers;"  each 
x,_|i,  is  a  strctchable  lever  that  pivots  about  the  fixed  axis  p,  and  nodes  x,_,  and  x,  arc 
constrained  to  slide  along  parallel  slots  L{x,_^)  and  Z.(x,)  respectively.  We  call  the  point 
Xq  (.  L  critical  with  respect  to  R,  if  the  (Z.,^,,^^,  .  .  .  ,p„)-path  determined  by  x^  has  the 
property  that  for  some  /,  cither 


(a)  x,_|X|  passes  through  some  corner  in  R  -  {  p  ^,  P2,  ■  ■  .  .p„},  oi 

(b)  X|_,,  X,,  X,.]  are  collinear. 

A  minimal  interval  I  ^  L  bounded  by  two  critical  points  is  called  a  critical  interval.  An 
(L,P|,p-,,  .  .  .  ,p^)-path  is  said  to  belong  to  a  critical  interval  /  if  its  first  node  Xq  is  in  /. 
The  concept  of  paths  and  critical  points  is  crucial  in  our  discussion  of  composition  and 
decomposition  of  balanced  chains. 

Lemma  5.4.  Given  a  line  L,  a  set  R  of  n  points,  and  a  sequence  (p^,  pj,  .  .  .  ,p„)  of  points 

from    R,    we    can    determine    all    the    critical    intervals    of  {L,p^,p2,  ■  ■  ■  iP m)'P'^'^^    with 
respect  to  R  in  0{n  log  n)  time. 

Proof.       First,     we     pick     an     arbitrary     point     Xq     on     L      and     compute     the     unique 

(Z.,p,,p,,  .  .  .  ,p^)-path,  (xq,  X| x„),  in  0(m)  time.    Then,  the  parallel  lines  Lo(x,), 

/■  =  0,  .  .  .  ,m,  can  be  easily  determined.    These  parallel  lines  divide  the  plane  into  stinps. 

For  each  r  £  R  -  {  p  ^,  pj p„  },  wc  determine  in  O(log  m)(  =  0(log  n))  time,  the  index 

/,  1  s  I  s  m,  such  that  r  and  p,  are  in  the  same  strip.  Then,  in  0(1)  time*  we  can  determine 
the  critical  pointy  €  Lq  corresponding  to  the  (Z.p,  p,,  P;'  •  ■  -Pm)'?^^^  ^^^^  passes  through  r. 
Also  for  each  /  =  1,  .  .  .  ,m-l,  in  0(1)  time  we  can  determine  the  critical  point 
corresponding  to  the  (Z.g,p,,p-,,  .  .  .  ,p„)-path  with  (x,. ,,  x,,  x,. ,)  being  collinear.  Since 
each  critical  point  takes  at  most  0(log  n)  time  to  compute,  0(n  log  n)  time  is  needed  in  total. 
D 

Definition.  Two  points  _v,z  on  the  line  Z.(x,),  for  any  /,  are  (L,P|,p-,,  .  .  .  ,p ^-equivalent , 
if  the  rwo  paths  through  y  and  i  both  belong  to  the  same  critical  interval     Let  t   be  a  simple 

balanced  (Lq,  L,)-chain  with  nodes  (xq,  x,,  .  .  .  .x,„)  and  pivots  (p,,P2 p„)-    The  chain 

C  determines  two  paths,  the  (Lq,  p  ^,  p  -, p^)-path  and  the  (Z.|,p„.,p„_|,  .  .  .  ,P|)-path. 

We  call  {L(^^,  p  ^,  p^,  ■  ■  .  ,p„)-path  the  forward  path  of  C  and  (Z.,,  p^,  p„_ ,,  .  .  .  ,p  ,)-path 
the  backward  path  of  C.  The  critical  interval  of  the  (Z,q,  p,,  p,,  .  .  .  ,p^)-path  containing  the 
first  node  x,^  is  called  the  forward  tolerance  of  C-    Similarly,  the  backward  tolerance  of  C  is 


This  is  due  to  the  fact  that  ^yx^p ,  is  similar  to  Arp,j  where  s  -  x,_|Xnz.(r). 
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the  critical  interval  of  the  (L,,  p „,  p„-^ p,)-path  containing  i„. 

Next,  we  describe  a  decomposition  process  that  transforms  a  path  determined  by  a 
simple  balanced  chain  until  the  path  becomes  critical,  or  in  other  words,  a  double-pivot  chord 
appears.  The  double-pivot  chord  is  then  used  to  divide  the  simple  balanced  chain  into  two 
simple  balanced  chains.    The  process  goes  as  follows: 

Let     C  =  (C| C„)      be      a     simple      balanced  (Lq,  L,,  /?)-chain      with     pivots 

(;7,,p, p^)     and     nodes     (xq,  x^ x„),     m  >  0.  The     chain     C     determines     an 

{Lq,p^,p., p„)-path,   (xj),  X,,  .  .  .  ,x„).    We  consider  moving  the  point  x^  along  Lq  in 

cither  direction  until  it  reaches  the  first  critical  point  Xq. 


Loi^k) 


L,    , 


Figure  5.10.  x^x,.,  passes  through  r 


Lei  (Xf),  X| x^)  be  the  (/,„,  p,,  p. p„)-path  at  x^^..    Then,  either 


(a)  X,  _  I JT,  passes  through  a  point  r  in  ^  -  (  p,,  pj p  „]  (see  Figure  5.10),  or 

(b)  X;_  I ,  x,,  ij .  I  are  collinear 

for  some  i,  0  <  i  s  m      We  assume  that  k  is  unique;  only  one  of  the  above  critical  rvenn 


occurs 


First    consider    possibility    (a).      Let    Z,,    be    the    directed    line     from    pj    to    r       For 
j=   1 *- 1 ,  let  C,  denote  the  chord  overlapping  X,.  I X, .    Observe  that 
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C  =  (c, c,.,) 

is  a  balanced  {Lq,  Z,2)-chain.    However 

(C,-, C„) 

does  not  represent  a  balanced  (Z.^,  L,)-chain,  because  while  the  pivot  p ^  is  at  the  center  of 
Jt„-iJ:„,  J„  is  not  on  Z.,.    To  obtain  a  balanced  (Lj,  Lj-chain,  we  continue  as  follows: 
Fori  =  k m,  let  jt,'  =  L^ix^CM^ix]). 


h  =  hi'k) 


hi't'O 


L,  =  Z.,(xJ  L^ixJ 


Figure  5.11. 


Let  C,  be  the  chord  overlapping  x,_  I  Jt, .    Then,  we  note  that 

c"  =  (c,., CJ 

is  a  balanced  (Z.,,  Z.|)-chain.  This  follows  from  the  fact  that  the  set  of  triangles  (see  Figure 
5.11) 

A  X,  X,  X,  (i  =  k.  .  .  .  ,  m) 

are  congrueni  Consider  the  (Z.-,,p^.|,  .  .  .  .p„)-palh  corresponding  to  J  :  It  is  imponan:  to 
see  from  our  construction  that  x^  and  x^  lies  in  the  same  critical  interval  with  respect  to 
(Z.,.  p, .  ,       .     .p^),  i.e.  as  we  move  from  i,  to  x^,  the  (Z.-,.  p^.  |.  .p,J-paths  encountered 
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along  the  way  arc  non-critical. 


Note  that  if  /■  €  ^t-jP*  then  the  truncated  version  of  L;  becomes  an  early  edge  of  C. 
Otherwise,  r  (  pjj:^  and  the  truncated  Z,;  ends  up  being  a  late  edge  of  C. 

Next,  we  consider  possibility  (b).    We  proceed  essentially  in  the  same  way,  except  this 
time  wc  set  i.-,  to  be  the  directed  line  from  p^  to  p^,  ,: 


For  I  =   1,  .  .  .  ,*,  let  C,  denote  the  chord  overlapping  -t,_,jt, .    Wc  immediately  get  a 
balanced  (Lq,  L2)-chain 

C    =  (C, C,_,). 

For  /  =  t+  1,  .  .  .  ,  m,  let 

x~,    =  Z.,(jr,)n/.,(x;) 


and  C,  be  the  chord  determined  by  x,_  |X, .    Then 

c"  =  (c;.2 cj 

is    a    balanced    (Lj,  Z.,)-chain,    if    x^^ ,     and    x^,,     lie    in    the    same    critical    interval    of 

(^:.  Pk-i Pj- 

Note  that  Lj  becomes  a  splittable  chord  after  the  above  adjustment  process. 

This  completes  our  process  of  decomposing  a  simple  balanced  chain  C-    The  resulting 
pair  of  simple  balanced  chains  C   and  C    is  called  the  decomposition  of  C.    It  should  be  noted 
that,  however,  each  of  c     and  C    does  not  in  general  enclose  the  set  R      Rather,  C    and  C 
enclose /f(Z,f^,  L,)  and  R{L-^,  Z.,)  respectively,  where 

r{Lq,  z.;)  =  /?nz.j  nz-f  and  R{L.,  L^)  =  rhl;  dl:  . 

We  have  shown  constructively  thai  every  simple  balanced  chain  of  length  m  >  0  can  be 
decomposed  into  two  simple  balanced  chams  of  length  <  m.  We  have  assumed  that  i  is 
unique  in  (a)  and  (b)  above.  It  is  not  hard  to  provide  the  necessary  modification  for  the  case 
where  more  than  one  critical  events,  occur.  When  that  happens,  C  is  decomposed  into  rw  o  (or 
more)  chains. 

5.5.    The  Composition  of  Simple  Balanced  Chain<) 

We  now  consider  how  the  process  in  the  previous  section  may  be  reversed,  i.e.  given 
two  balanced  chains  C    and  C  .  finding  out  whether  a  chain  C  can  be  composed  from  J"    and 
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<?      P 

-X X- 


Figure  5.12.  Composing  an  (Lq,  L-^)-ch&m  and  an  (Lj,  Z.,)-chain 

Let  C  be  the  balanced  (Lq,  Lj,  R  )-chain  and  C  be  the  balanced  (1,2,  Z,,,  R  )-chain.  The 
backward  supporting  line  of  C  coincides  with  the  forward  supporting  line  of  C  .  Let  R  and 
R  be  subsets  of  some  point  set  R.  L,  is  a  directed  line  from  p  to  q  where  p  and  q  are  points 
from  R  and  Lq  precedes  Lj  ^nd  L,  precedes  L,.    See  Figure  5.12.    Let  m   be  the  length  of  C  , 

m    be  the  length  of  C  ,  and  let  x, ,  /  =  0 m,  and  x, ,  /  =  0,  .  .  .  ,  m    be  the  nodes  of  C 

and  C    respectively. 


^i(-^-i)     ^o(-fm) 


Figure  5.13.  Cuuiposition  through  an  early  edge 


If  C  and  C  form  the  decomposition  of  some  simple  balanced  (Lq,  L,, /?)-chain,  then 
that  can  only  come  about  by  the  decomposition  process  in  one  of  the  following  three  ways 
corresponding  to  possibilities  (a)  and  (b)  described  in  the  previous  section: 


(A)    x^  Xq  is  an  early  edge.    (Sec  Figure  5.13. 
Let  m  =  m  *  m  *  1  and  set 
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P,  = 


p, 

U  i  =   1,  .  .  . 

9 

if  1  =  m  +  1 

P.-m- 

if  1  =  m  +2, 

Let  j1[  =  x'^lq  and  jt^  =  Z-o^'m  )'^^i('-i)-  '*  '*  clear  from  the  decomposition  process 
that  C   is  the  composition  of  t     and  C  ,  if  (i)  x„  Xq  is  an  early  edge,  (ii)  x^    and  x„    are 

(i-o.Pi.P: p„  )-cquivaIcnt,    and    (iii)    x„     and    xl,    arc    (L,,  p^,  p„  _,,...  ,p, ,  ,7)- 

equivalent. 

As  a  consequence  of  the  balancing  of  the  edges  of  C  ,  the  critical  interval  of  the 
(Lq,  PpPi,  •  .p„  )-path  containmg  x^  can  be  obtained  by  a  translation  of  the  forward 
tolerance  of  C  that  maps  Xg  to  x„  .  Thus,  (ii)  can  be  verified  in  0(1)  time,  if  the  forward 
tolerance  of  C  is  available.  Similarly,  we  can  do  the  same  with  (iii),  using  the  backward 
tolerance'  of  C  ■ 

h  'm-l      *0    1    P  'm- 


Figure  5.14    Composition  through  a  late  edge 


(B)     x^  Xq  is  a  late  edge     (See  Figure  5.14.) 
Let  m  =  m  ■>-  m   -^  I  and  set 


P,  = 


P, 

P 

P 


/  -  m  -  I 


if  /  =    1 m 

if  «■  =   m  ^  \ 

if  /  =   m  +  2,  .  .  .  ,  m. 


Let  x^  .,  =  x„  |p  and  x„  .,  =  Z.,(xp)  nLQ(x„  .,).  Then  C  is  the  composition  of  C  and 
C  ,  if  the  following  requirements  arc  met:  (i)  x„  Xq  is  a  late  edge  (ii)  x„  .  ,  and  x^  .  .  are 
(Z-Q.  p,,P;.  .        ,p^  .  p  )cquivaleni,    and    (iii)   x„  . ,    and   x^,    arc    (L^,  p^  ,  p^  _^ P,)- 


'    Notice    that   the    critical   interval   on   /.|(x_|)   containing  x_,    can    be   obtained    from    that 
containing  x^,  by  a  translation  that  maps  X(,  to  x_ , . 
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equivalent. 


^o('m-i)      ^i(->^'-i) 


■I     '■m   '\ 


Figure  5.15.  Composition  by  splitting 


(C)    x^Xq  is  a  splittable  edge.    (See  Figure  5.15.) 
Let  m  =  m  +  m  +2  and  set 


P 
<? 
P. 


if  1  =   1,  ....  m 
if  /  =  m  +  1 
if  I  =  m  +  2 
if  /  >  m  +  2 


Let  x^  .,   =  x,„  \p    and  x,,  =  x^\q ,   and  let  x^  . ,  =  Z.^(x.  ,)ni.,(x,^^  ,).    Then,  as  the 
third  possibility,  C  is  the  composition  of  C    and  C  ,  if  (i)  x,„  Xq  is  a  splittable  edge,  (ii)  x^ 


and 


and     x,^.,      are      (L^,  p , ,  p,,  .  .  .  ,p„., /5)-equivalent,      and      (iii)     x„ 
(I.,,  /^^  .  p„  _, P|,  <?)-equivalent. 

To  summarize,  we  note  that  if  one  of  the  three  conditions  above  is  satisfied,  then  J  and 
C  form  the  decomposition  of  C-  Otherwise,  there  is  no  way  in  which  we  can  compose  a 
balanced  chain  from  C    and  C  ■    Hence, 


Lemma  5.5.  Given  the  simple  balanced  (Lq,  Z,-,,  R  )-chain  C  ,  the  simple  balanced 
(Z.-,,  /.,,/?  )-chain  C  and  their  forv.ard  and  backward  tolerances,  it  takes  0(1)  time  to 
verify  whether  C  and  C  can  be  composed  to  create  the  simple  balanced  (Lq,  L,,  R \)R  )• 
chain  C- 
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5.6.    An  Algorithm  for  Finding  the  Balanced  Chain 

Now  wc  arc  ready  to  describe  an  0{n')  time  algorithm  for  finding  the  balanced  chain 
that  is  supported  by  a  pair  of  given  directed  lines  (Lq,  L^)  and  encloses  a  set  of  given  points 
R. 

Notation.  Let  R-  denote  directed  lines  L  defined  by  an  ordered  pair  of  points  in  R  such  that 
Lq  precedes  L  and  L  precedes  Z,,.  Wc  only  deal  with  lines  with  directions  between  Lq  and  Z,,, 
since  only  they  arc  relevant  to  an  upward  convex  (Lq,  Z.,,  >?)-chain. 

The  algorithm  is  a  two-step  process.  First  we  compute  all  the  simple  balanced  chains 
supporting  a  pair  of  lines  from  R'U{  L^,  L^  }.  If  there  is  a  simple  balanced  {L^,  Z.|)-chain 
then  we  are  done.  Otherwise  the  balanced  (Lq,  Z.,)-chain  is  not  simple.  In  that  case  there  is 
more  than  one  way  of  computing  the  balanced  chain.  Wc  describe  one  that  is  based  on  graph 
search.  From  the  simple  balanced  chains  we  have  computed,  we  can  define  a  digraph  G 
whose  vertices  correspond  to  the  lines  in  {  L^,  Z.,  }U/?-.  Let  there  be  an  edge  in  C  between 
the  two  vertices  representing  lines  C  and  C  ,  if  there  exists  a  simple  balanced  (C,  C  )-chain. 
A  graph  search  on  G  is  then  conducted  to  find  the  path  connecting  the  nodes  representing  Lq 
and  Z.,  which  corresponds  to  the  balanced  (Z,q,  L|)-chain. 

Algorithm  5.1.  The  Balanced  Chain 

Step  J  For  each  pair  of  lines  (C,  C  )  in  /?-x/?',  find  /?(C,  C  ),  the  set  of  points  in  R  that 
is  in  C"nC".  If/f(C,C)=0,thenwe  have  a  simple  balanced  chain  of  length 
0.  Otherwise,  Z?(C,  C  )  is  not  empty  and  we  have  a  smallest  corner  triangle 
problem.  Let  the  smallest  corner  triangle  be  determined  by  a  chord  C  .  If  C  is 
single-pivot  chord  then  it  is  balanced  with  respect  to  (C ,  C  ).  Surely  we  have 
found  a  simple  balanced  chain  of  length  1  Otherwise  C  is  a  double-pivot  chord, 
and  it  may  or  may  not  be  balanced  according  to  our  new  definition  of  balancing." 
If  C  is  not  balanced,  it  is  clear  that  the  balanced  (C,  C  )-chain  is  either  not  simple 
or  is  of  length  2  or  more. 

Step  2       Repeat  Step  .?  until  no  new   chain  is  generated. 

Step  3  For  each  pair  of  extremal  chords  (C,  C  )  €  (  Z.^,,  Z.,  }U>f-,  if  the  simple  balanced 
(C,  C  )-chain  has  not  been  generated,  wc  will  attempt  to  generate  it  by  composing 
simple  balanced  chains  computed  previously  For  that,  we  iterate  th'ough  all 
chords  C  €  /?-  (C  precedes  C  and  is  preceded  by  C),  checking  whether  the 
(C,  C  )-  and  (C  ,  C  )-chains  have  been  computed  and  whether  they  form  a 
decomposition  of  the  (C .  C  )-chain.  If  so,  the  (C,  C  )-chain  is  constructed.  In 
order  to  facilitate  subsequent  composition  of  chains,  we  compute  the  forward  and 


It  is  always  "balanced"  in  the  old  sense  that  its  tenter  is  on  an  edge  of  the  convex  hull  of 
/f(C.  C  ). 
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backward   tolerances   for   the   (C,  C  )-chain   with   respect  to   the   enlarged   set  of 
points,  R(C,  C  ),  right  after  it  has  been  constructed. 

Step  4  After  we  have  computed  all  simple  balanced  chains,  if  there  is  a  simple  balanced 
chain  supported  by  (Lq.  ^i)  t^en  return  the  chain. 

Step  5  Otherwise,  the  balanced  chain  is  not  simple.  Construct  a  graph  G  with  nodes 
corresponding  to  {  Lq,  L,  ]UR'.  If  there  exists  a  simple  balanced  chain  supported 
by  a  pair  of  lines,  then  put  an  edge  between  their  corresponding  nodes.  Then 
conduct  a  depth-first  search  on  G  to  find  all  the  paths  between  the  two  nodes  g 
and  g  corresponding  to  Lq  and  Z,,  respectively.  The  verification  of  the  balancing 
condition  of  a  double-pivot  chord  can  be  done  when  post-visiting  its  corresponding 
node.  According  to  our  uniqueness  and  existence  theorem,  one  and  only  one 
balanced  chain  will  be  found,   n 

To  analyze  the  complexity  of  this  procedure,  we  divide  the  cost  into  three  parts: 

(1)  The  cost  of  verifying  the  possibility  of  composition. 

(2)  The  cost  of  actually  composing  the  chains  and  computing  the  forward  and  backward 
tolerances. 

(3)  The  cost  of  searching  the  graph. 

First  let  us  count  the  cost  of  part  (1).  Step  1  takes  0{n^\og  n)  time,  since  there  are 
O^n")  pairs  of  lines  to  go  through  and  each  pair  takes  0{n  log  n)  time  to  process.  Each  time 
Step  2  IS  executed,  we  go  through  all  triples  (C,  C  ,  C  )  verifying  whether  it  is  possible  to 
compose  the  (C,  C  )-  and  (C  ,  C  )-chains.  So,  there  arc  O(n^)  instances  of  testing  for  a 
possible  composition  with  each  test  costing  0(1)  time  (from  Lemma  5.5).  Since  we  repeat 
Step  2  at  most  n  times,  part  (1)  takes  O(n^)  time  in  total. 

Because  of  the  uniqueness  of  the  balanced  chains,  there  are  at  most  O(n')  chains  that 
are  composed  during  the  entire  procedure.  The  cost  of  part  (2)  is  C)(n-log  n),  since  the 
computation  of  tolerances  for  each  chain  takes  0{n  log  n)  time  (from  Lemma  5.4)  and  the 
construction  of  a  chain  from  two  smaller  chains  takes  only  linear  time. 

The  cost  of  graph  search  and  balancing  condition  verification  of  part  (3)  is  O(n'),  since 
the  size  of  graph  is  only  0(n").  As  for  the  storage  requirement,  since  there  are  O(n^) 
balanced  chains  each  of  which  takes  at  most  0{n)  space  to  keep  track  of,  we  need  0{n^) 
storage  space.    Thus 

Theorem  5.6.    The  balanced  chain  problem  can  be  solved  in  0(n'')  time  using  0{n^)  space. 

In  the  next  chapter,  we  will  deal  with  the  potato-peeling  problem,  the  problem  of 
finding  a  minimum  area  convex  polygon  contained  in  a  given  non-convex  polygon  There  are 
some  fundamental  facts  about  the  solution  of  this  problem.  One  of  them  is  that  all  the  edges 
of  the  minimum  area  included  polygon  are  balanced.    Finding  these  balanced  edges  lies  right 
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at  the  heart  of  the  solution  to  the  potato-pccling  problem.  The  computation  we  described  in 
this  chapter  for  finding  balanced  chains  will  be  used  with  slight  modification  in  solving  the 
potato-peeling  problem. 

5.7.    Unsolved  Problems 

We  can  generalize  the  smallest  corner  triangle  problem,  in  another  waj;  instead  of 
asking  for  a  chain  of  any  number  of  chords,  we  insist  that  the  chain  consist  of  k  chords  k  >  1. 
Then  we  have  the  following  unsolved  problem: 

Problem  5.3.  Given  a  pair  of  directed  lines  (Lq,  Z,,)  (Z,q  precedes  Z,,),  a  set  of 
points  R  in  L^  DZ,,"  ,  and  a  fixed  number  i  >  1,  find  an  (Lq,  Z,,,  /?)-chain  of  length 
k  with  the  minimum  area. 

It  would  also  be  interesting  to  know  some  lower  bounds  for  the  balanced  chain 
problem. 


CHAPTER  6 


The  Potato-Peeling  Problem 


As  J.  Goodman  observed  in  [Goodman  (1981)],  the  usual  way  in  which  a  swivel-type 
potato  peeler  is  used  to  produce  a  convex  peeled  potato  suggests  the  problem  of  maximizing 
the  volume  of  a  convex  body  contained  in  a  given  non-convex  body.  In  this  chapter,  we  look 
at  the  planar  case  of  finding  a  maximum  area  convex  subset  contained  in  a  planar  n-sided 
simple  polygon. 

6.1.   Statement  of  the  Problem 

Given  a  non-convex  polygon  P ,  we  are  to  find  any  maximum  convex  subset  Q  contained 
in  /•;  it  is  intuitively  clear,  but  rigorously  proved  in  [Goodman  (1981)],  that  Q  [%  a  convex 
polygon.    So,  we  can  state  the  problem  more  specifically: 

Problem  6.1.    Given  a  non-convex  polygon  P,  find  a  convex  polygon  Q  contained 
in  P  with  maximum  area. 


'n-I 


Figure  6.1. 


For  the  rest  of  the  chapter,  we  assume  that  we  are  given  a  fixed  but  arbitrary  n-s:ded 


polygon       P  =    (vf,, 


.  v„ 


has      k  >  Q      reflex       corners: 


,  w 


(0  s  i|  <  1,  <  ,  .  .  .  ,    <  (j  <  n).     Vi'e    will    write    u^   for    v.     For   0  s  i  <  n ,    let    e     be    the 
directed  edge  of  P  from  v,  to  v,_  |,^^^  „,.    For  simplicity,  we  assume  that  no  three  corners  are 
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collincar.  (Sec  Figure  6.1.)  As  for  the  solution,  we  assume  that  a  certain  convex  polygon  Q 
=  (**'o'**'i'  •  •  •  ■^m-\)  ^^^  ^^^  maximum  area  among  all  convex  polygons  included  in  P.  Let 
J,  denote  the  directed  edge  of  Q  from  w,  to  >*',-i(mod„,.  (0  s  ;  <  m).  Note  that  there  may  be 
more  than  one  maximum  area  included  polygon,  so  when  we  say  that  Q  has  a  certain 
property,  we  mean  that  at  least  one  such  polygon  has  the  property. 

6.2.    A  Fioiteness  Criterion 

First,  wc  give  some  definitions  which  arc  derived  from  the  previous  chapter: 

Definition.  Wc  define  a  chord  of  P  to  be  a  line  segment  satisfying  the  following:  (i)  It 
passes  through  at  least  one  corner  of  P.  (ii)  It  is  fully  contamed  in  P.  (iii)  Its  two  endpoints 
are  on  the  boundary  of  P.  An  extremal  chord  is  a  chord  that  passes  through  two  corners  of 
P.  In  particular,  an  edge  of  P  is  an  extremal  chord.  If  v^^  is  fully  contained  in  P,  then  let 
£■,  denote  the  extremal  chord  that  overlaps  v/v'.  and  is  directed  from  v^  to  v  Let  E,  be  the 
set  of  all  maximum  line  segment  contained  in  P  that  passes  through  v,  but  not  v,.,.  Let  E' 
be  the  set  {  £,  ^  |  £,  ^  is  defined  }  if  v,  is  a  reflex  corner,  and  the  set  {  £,  ,, ,  },  otherwise. 
Clearly  |E,*|  <  n.  It  is  easy  to  compute  all  the  sets  H,*  in  O(n')  time.  Let  E*  denote  the 
union  of  all  the  E,*'s. 

Intuitively,  we  must  cut  off  all  the  non-convex  corners  of  P  with  straight  lines  in  order 
to  get  a  convex  polygon  Q.  Furthermore,  each  of  these  lines  must  touch  one  or  two  reflex 
corners  of  P,  for  otherwise  we  can  move  the  line  away  from  the  center  of  Q  to  increase  its 
area.  So,  each  edge  of  Q  overlaps  cither  one  of  those  corners  or  an  edge  of  P .  In  other 
words,  Q  is  the  intersection  of  m  half-planes  defined  by  m  chords  (C,,  C-,,  .  .  .  ,C^)  of  P 
such  that  each  C,  i  E  ,  for  some  j,  0  ^  j  <  n.    Recall  that  for  a  directed  line  L .  L'  denotes 

the  closed  half-plane  to  the  left  of  the  line.    Thus,  Q  =   (Hc,'). 

I- 1 

Definition.  If  C  is  a  chord  of  P  thai  overlaps  an  edge  s  oi  Q  and  C  has  the  same  direction  as 
s,  then  we  say  that  C  determines  s.  We  call  an  edge  of  Q  extremal  if  it  is  determined  b\  an 
extremal  chord   of  P      If  for  0  s  <  <  m,  C,  determines  -f,- fin,„(j  „,  where   c   is  some   fixed 

number,  then  we  say  that  J  ^    (C^,C■^ C„)  determines  Q     Recall  the  definition  of  a 

sequence  of  directed  line  segments  being  upward  convex  and  enctosin/i  a  set  of  points  from 
the  previous  chapter.    For  two  chords  Cp,  and  C,.,  of  P  (wc  do  not  require  thai  C    precedes 

C,.|),   Cfi  (  E'   and   C,. ,  (.  E',   an   (C ,  C  )-chatn    C   is  a  sequence  of  chords   (C, C,) 

supported    by     (C .  C  )    which    is    upward    convex,    and    encloses    i  ^' .  i v,.,  )       For 

1  =  0 m,  let  X    be  the  intersection  of  C,  and  C,.,.    As  in  the  previous  chapter,  we  call 

the  X,  nodes  and  call  x  ,j:  an  edge  of  C  or  the  truncated  version  of  the  chord  C ,.  For  a  chord 
which  is  not  an  edge  of  P,  the  definition  for  balanced  chord  is  the  same  as  that  of  Chapter  5, 
with  reflex  corners  playing  ihc  role  of  the  point  set  R     Any  chord  which  is  also  an  edge  of  P 
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is  balanced.^   Balanced  chains  and  simple  balanced  chains  arc  similarly  defined. 


Figure  6.2.  A  V-chain  and  an  i4 -chain 
(showing  only  the  truncated  versions  of  the  chords) 


If  a  chain  is  supported  by  a  pair  of  V-shaped  lines,  then  it  is  a  V-chain.  Otherwise  it  is 
an  A-chain.    (See  Figure  6.2.) 

As  mentioned  before,  the  smallest  corner  problem  is  a  subproblem  in  the  potato-peeling 
problem  and  is  equivalent  to  the  balanced  chain  problem.  So,  it  is  not  surprising  that  for  the 
potato-peeling  problem,  we  derive  a  finitcness  criterion  based  mainly  on  the  notion  of 
balancing. 

In  the  following,  wc  give  three  characterizations  of  Q: 

Lemma  6.1.  Every  edge  of  Q  is  balanced. 
Proof.    By  arguments  similar  to  Lemma  5.1.  D 


This  is  consistent  with  the  original  definition  of  balancing  in  Chapter  2. 
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C(x,) 

Figure  6.3.  The  case  where  /  =  2; 
the  dotted  lines  represent  the  chain  after  perturbation 


Lemma  6.2.  Let  i,  and  s  be  two  non-consecutive  edges  of  Q  which  are  determined 
respectively  by  chords  C  and  C  of  P .  If  (C ,  C  )  is  a  V-pair,  then  there  exists  an  extremal 
edge  s^  in  Q ,  such  that  i  <  a  <  7  (mod  m). 

Proof.    For  the  purpose  of  contradiction,  let  us  assume  that  there  is  no  extremal  chord  in  the 

(C,  C  )-chain,  C.    In  other  words,  C  is  a  simple  balanced  chain.    Let  (xq,  jt, Xj)  be  the 

nodes  of  C,  and  (/>,,  p,,  .  .  .  . p ,)  be  the  single-pivots  of  C.    Clearly,  (x^.  or, x,)  is  a 

non-critical    (C,P|,  .        ,p)-paih      So,    we    can    infinitesimally    perturb    Xr    along    C    while 
keeping    C    upward   convex    and    enclosing    the    same   set  of  points.     Let   the    result  of   the 

perturbation   be   (xp,  x, x,.)   and  let   C    =  (C, C,)  where   C     be   the   chord   that 

overlaps  x,_  |X  ,,  1  =  1 /.    (Sec  Figure  6.3  ) 


Let  >■  =  C(x)nx^|t  .  If  C  and  C  are  not  parallel,  then  C(x)  *  C  and  y  ^  x  .  It  is 
easy  to  see  that  C  contributes  a  gain  in  the  area  of  Q  over  C  that  equals  the  area  of  Ax;X:V 
(see  Figure  6.3),  contradicting  the  optimality  of  Q  Finally,  if  the  two  lines  C  and  C  are 
parallel,  then  C(x)  overlaps  C  ,  v  =  x,,  and  the  area  of  Q  remains  the  same  after  the 
perturbation  Hence,  we  can  perturb  the  chain  without  decreasing  the  area  of  Q, 
contradicting  the  optimalit\  of  Q.  We  conclude  that  in  all  cases,  there  is  at  least  one 
extremal  chord  in  C.  Hence,  there  is  at  least  one  extremal  edge  s^  in  Q.  and  from  the 
convexity  of  Q .  1  <  a  <  y(mod  m).    0 
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Lemma  6.3.  There  are  at  least  two    extremal  edges  in  Q. 

Proof.    Let  C  =  (Cq,  C, C^.,)  be  the  chain  that  determines  Q.    Then  clearly,  the  chain 

C  =  (C|,  Cj C„.j)  is  a  V-chain  supported  by  (Cq,  C„.,).    By  Lemma  6.2,  there  is  an 

extremal  chord  in  C  ,  say  C,,  1  s  /  s  m-2.    Since  (C,.,,  C,.2.  •  •  •  .C",_2)  is  again  a  V-chain 
supported  by  (C,,  C,_,),  there  is  extremal  chord  other  than  C,.    □ 

To  summarize,  we  characterize  the  solution  polygon  Q  as  follows: 

(1)  Q  has  at  least  two  extremal  edges. 

(2)  Each  pair  of  (counterclockwise)  adjacent  extremal  edges  of  Q  is  A -shaped. 

(3)  The  chain  supported  by  a  pair  of  adjacent  extremal  edges  of  C  is  a  simple  balanced 
chain. 

6.3.    A  Potato-peeling  Algorithm 

The  lemmas  from  the  previous  section  provide  us  with  a  finiteness  criterion  in  the  sense 
that  there  are  finite  number  of  possibilities  for  the  extremal  chords  in  the  chain  that 
determines  Q.  We  can  have  an  exponential  time  algorithm  which  computes  Q  by  guessing 
these  extremal  chords,  provided  we  can  also  compute  the  simple  balanced  chain  supported  by 
two  consecutive  extremal  chords  in  exponential  time. 

In  fact,  we  can  do  better.  These  simple  balanced  chains  can  be  computed  in  polynomial 
time  by  (some  modification  of)  the  computation  for  the  smallest  corner  problem  described  in 
the  previous  chapter.  In  addition,  by  organizing  the  guesses  for  extremal  chords  present  in 
Q,  eliminating  redundant  guesses,  the  whole  problem  can  be  solved  in  polynomial  time.  This 
section  describes  such  an  algorithm. 

6.3.1.    The  Dynamic  Programming  Approach 

We  will  again  take  the  dynamic  programming  approach  to  organize  our  computation. 
In  this  subsection,  we  will  show  how  to  break  up  the  potato-peeling  problem  into 
subproblems,  and  how  to  divide  the  area  of  Q  arzcng  the  constituent  subproblems.  Since 
there  are  a  finite  number  of  extremal  chords,  it  is  only  natural  that  we  divide  the  solution 
polygon  Q  along  the  lines  of  extremal  edges  of  Q.  We  choose  a  certain  point  on  the  chord 
that  determines  an  extremal  edge  for  this  purpose: 

Dennition.    For  an  extremal  chord  C   determined   by   two  corners  v,  and  v    of  P,  call  the 
midpoint  between  v,  and  v  ,  the  reference  point  of  C. 

Let  Q  be  determined  by  a  sequence  of  extremal  chords  including  C  and  C  .  We  can  cut 
the  (C,  C  )-chain  of  Q  along  the  line  joining  the  respective  reference  points  c  and  r   of  C  and 


"  We  can  actually  prove  a  stronger  result  -  three  extremal  chords.    But  for  our  purpose,  two 
are  enough. 
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C  .  This  may  result  in  two  non-simple  polygons.  Therefore,  we  need  a  way  of  measuring 
area  such  that  (i)  nonsimple  polygons  can  be  accounted  for,  and  (ii)  when  we  "join"  two 
polygons,  the  areas  will  add  up  to  the  area  of  the  resulting  polygon. 

Derinitlon.    Let  T  =  (Iq,  f, 'n-i)  ^^  ^  polygon  that  is  not  necessarily  simple,  with  n 

corners  where  (x,,  y,)   is  the  coordinate  of  corner  r,,  x  =  0 n-\.    The  area^  of  T  is 

defined  as  follows: 

^  .-0 

If  T  is  simple  polygon,  then  this  definition  gives  the  expected  notion  of  area,  with  a 
positive  sign  if  the  corners  arc  given  in  counterclockwise  order,  and  negative  otherwise.  If  T 
is  not  simple,  then  7  defines  more  than  one  region  in  the  plane.  If  we  traverse 
counterclockwise  along  the  boundary  of  T,  Some  of  these  regions  lie  to  the  left  of  the 
traversal  and  others  to  the  right.  This  definition  of  area  measures  the  summation  of  the 
areas  of  these  regions  with  the  regions  to  the  left  carrying  a  positive  sign  and  the  regions  to 
right  a  negative  sign.  It  is  easy  to  see  when  we  "join"  two  nonsimple  polygons  along  a 
common  edge,  the  areas  add  up  nicely.    More  specifically,  by  "joining"  polygon 

A   =  (r,,  a,,  a, a„,  r-,) 

and  polygon 

B   =   {r,,b„  b. fc„,  r,) 

"along"  r^r^,  we  get  polygon 

C  =   (a,,  fl; a„.  fc,,  b. bj. 

And  the  definition  guarantees  that 

Area{C)  =   Area{A)  ~  Area{B). 

Dennition.  For  C  =  f  ^  and  C  =  E,  ,,  let  P(C .  C  )  denote  {  vv,  \  j  <  a  <  *(mod  n)  K  the 
set  of  points  counterclockwise  from  v.  to  v,  (exclusively).  A  pair  of  chords  (C,  ,  C^  )  is 
admissible  if  i  <  j  <  k  <  /(mod  n)  A  (C,  ,,  C,  ,)-chain  is  also  called  admissihU  if 
(C,  ,,  C ^  ,)  is  admissible  and  the  first  and  last  nodes,  x,~^  and  x„  (for  some  m),  lie  in  £,  and 
E^  ,  respectively.    A  (C,  C  )-chain  C  and  a  (C  ,  C  )-chain  C   are  compatible  if 

(1)      Both  J  and  J   are  admissible. 


'  This  is  called  the  signed  area  in  [Lopshits  (1963))  and  as  pointed  out  to  us  by  Professor 
Micha  Sharir  can  be  derived  by  summing  up  over  (,  %(-t_-i  .  |)(y,+ y,.  |).  the  area  of 
trapezoid  determined  by  (i,,  y,),  (x,.  ,,  y,. ,),  and  their  projections  on  the  x-axis. 
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(2)     The  last  node  x„  (for  some  m)  of  C,  the  first  node  Xq  of  C   and  the  pivots  {p ,  q)  of  C 
occur  in  the  following  order: 

6.3.2.    An  Algorithm 

The  following  algorithm  for  the  potato-peeling  problem  proceeds  in  three  main  steps. 
First,  all  the  simple  balanced  chains  supported  by  compatible  pairs  of  extremal  chords,  are 
computed.  Secondly,  a  dynamic  programming  process  is  performed  to  find  the  partial 
solution  to  the  subproblem  of  finding  the  "best"  chain  supported  by  a  pair  of  extremal 
chords.  And  finally,  a  maximization  process  over  the  possible  combinations  of  these  partial 
results  is  done  to  find  the  overall  solution.    The  algorithm"  goes  as  follows: 


Figure  6.4. 


"  We  only  describe  the  computation  of  the  area  of  Q  instead  of  computing  the  convex 
polygon  itself.  It  is  easy  to  modify  the  algorithm  to  compute  the  polygon  in  addition  to  its 
area. 
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MIC  ,  C)      - 


A  cc  c 


Figure  6.5. 


Algorithm  6.1.  Potato-Peeling  -  Area  Measure 

Step  1  To  compute  the  simple  balanced  (C,  C  ,  P{C,  C  ))-chain  (when  one  exists)  for  all 
compatible  pairs  of  extremal  chords  (C ,  C  ),  we  go  through  essentially  the  same 
balanced  chain  computation  as  in  the  previous  chapter.  The  area  of  the  simple 
balanced  (C,  C  )-chain  is  stored  in  A(C,  C  ).  -4  is  a  matrix  indexed  by  pairs  of 
chords.    The  differences  can  be  summarized  as  follows: 

(1)  In  the  place  of/?-,  we  now  have  5. 

(2)  When  computing  the  simple  balanced  (C,  C  )-chain,  we  check  the  possibility 
of  composing  the  simple  balanced  {C ,  C  )-  and  (C  ,  C  )-chains,  for  all 
C     €   =;,  /  <  *  <  ;. 

(3)  The  set  of  points  involved  when  we  talk  about  simple  balanced  (C,  C  )-chain 
is  P(C,C  )  instead  o{  R{C ,  C  )  =  R  r\C'  nC  '  . 

(4)  During  the  construction  of  simple  balanced  chains,  we  also  compute  their 
area^ 


Step 


For  a  pair  of  admissible  extremal  chords  (C,  C  ),  we  compute  the  maximum  area 
of  all  chains  supported  by   (C,  C  ) /"    This  maximum  area  is  stored  in  M(C,  C  ) 


where  M   can  be  computed  in  stages  where  in  stage  s   (s 
compute  the  entries  M(C,  C  )  with 


1, 


1),  we 


C  i  ~' .    C    €  E*    and   j-i  ^   5 (mod  n) 
Note  that  m  stage  j  =    1,  we  have  ;  =   /-  l(mod  n)  and  the  constraint  that  (C,  C  ) 


'  This  will  only  takes  0{n)  time  for  each  simple  balanced  chain  constructed     Thus  the  total 
running  time  for  computing  simple  balanced  chains  remains  the  same  -   0(n  ) 

^  M(C,  C  )  =   -<*>,  if  [C ,  C  )  is  not  an  admissible  pair  of  chords. 


77 


forms  an  admissible  pair  implies  that  C  must  be  the  chord  E , j.  Let  c  and  c  be  the 
reference  points  of  C  and  C  respectively.  Then  M{C ,  C  )  is  the  area  of  Acv^c  . 
(Sec  Figure  6.4.)    For  stage  j  >  1,  if  there  is  a  simple  balanced  (C,  C  )-chain  then 

M{C,  C  )  =  A{C,  C  ). 

Otherwise,  M{C ,  C  )  is  determined  by  a  chain  that  has  an  extremal  chord  and  can 
be  computed  as  follows:^ 

W(C,  C  )  =   max    {  M (C ,  C  )+  MiC  ,  C)  +  Area{Accc)  }, 

c' 

where  C  ranges  over  E^*,  /  <  *  <  ;'(mod  n)  and  c  is  the  reference  point  of  C  . 
(Sec  Figure  6.5.) 

Step  3       Since  there  are  at  least  2  extremal  edges  in  ^,  the  area  of  Q  is* 

max   {M(C,  C)  +  M{C  ,C)  \C  €  E,',  C    €  Ej ,  for  all  i  and;  }.  D 

6.3.3.    Analysis  of  Complexity 

From  the  previous  chapter,  we  know  that  Step  1  takes  Oin^)  time  and  0{n^)  space.  In 
Step  2,  we  spend  0{n')  time  to  carry  out  the  maximization  for  each  entry  of  M .  Thus  the 
whole  matrix  M  takes  0{n^)  time  in  total  to  compute.  Each  entry  of  M  requires  0{n)  space. 
In  total,  ©(n-)  space  is  needed  for  the  storage  of  the  matrix  M.  And  finally.  Step  3  takes 
0{n')  time  and  does  not  require  more  than  0(n)  space  to  keep  track  of  the  best  solution. 
Puning  all  the  costs  together,  we  have 

Theorem  6.4.    The  potato-peelmg  problem  can  be  solved  in  0{n  )  time  using  O(n-)  space. 

6.4.    Unsolved  Problems 

Consider  the  three-dimensional  version  of  the  potato-peeling  problem: 

Problem  6.2.    Given  a  non-convex  polytope  P,  find  a  convex  polytope  Q  contained 
in  P  with  the  maximum  volume. 

The  balancing  property  of  the  potato-peeling  problem  in  two  dimensions  generalizes  to 
three-dimensional  space.  More  specifically,  every  face  of  the  maximum  volume  convex 
polytope  contained  in  the  given  non-convex  polytope  P  is  balanced  in  the  sense  that  the 
geometric  center  of  the  face  is  in  P  (see  [Steward  (1942).  Levy  (1939)]  ).    Besides  that,  we  do 


^  Note  that  we  do  not  check  compatibility  between  the  chains  that  contribute  the  areas 
M{C,C  )  and  M(C  ,  C  ).  If  they  are  not  compatible,  then  they  cannot  form  an  optimal 
chain.    Thus,  it  is  not  necessary  to  exclude  them  from  the  maximization 

"  Notice  again  that  we  do  not  check  compatibility  beDAcen  A/(C,  C  )  and  M{C  ,  C)  for  exactly 
the  same  reason  as  in  step  2. 
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not  know  anything  else  that  will  help  us  in  solving  the  problem. 

Unlike  for  the  minimum  area  enclosing  polygon  problem,  in  the  formulation  of  the 
potato-peclmg  problem,  we  stated  that  the  largest  convex  included  polygon  is  to  be  found,  no 
matter  how  many  sides  it  has.  An  obvious  variation  of  the  potato-peeling  problem  is  to  find 
the  maximum  area  triangular  (or  quadrilateral)  potato.    In  general, 

Problem  6.3.    Given  an  n-sidcd  non-convex   polygon  P   in  the  plane  and  a  fixed 
number  t,  3  ^  *  s  n,  find  a  convex  i-gon  Q  included  in  P  with  the  maximum  area. 

In  the  next  chapter,  we  will  study  the  potato-peeling  problem  with  perimeter  as 
measure. 


CHAPTER  7 


Potato-Peeling  —  Perimeter  Measure 


In  the  previous  chapter,  we  studied  the  potato-peeling  problem,  the  problem  of  finding 
maximum  area  convex  polygon  included  in  a  given  non-convex  polygon.  In  this  chapter,  we 
look  at  the  same  problem  except  now  the  size  of  the  included  polygon  is  measured  by  its 
perimeter.  The  potato-peeling  problem  under  the  perimeter  measure  can  be  solved  by  using 
essentially  the  same  dynamic  programming  technique.  However,  a  stronger  finiteness 
criterion  than  the  one  we  have  with  area  measure  makes  possible  a  much  simpler  algorithm 
and  a  faster  running  time,  0(n*). 

7.1.    A  Finiteness  Criterion 

Lemma  7.1.  Lei  C  =  (C,,  C,  .  .  .  ,C^)  be  a  chain  that  determines  the  maximum  perimeter 
convex  polygon  Q  included  in  the  given  polygon  P  For  each  i,  \  S.  i  ^  m.  C ^  is  extremal. 
Furthermore ,  double-pivot  chords  must  be  balanced. 

Proof.    There  are  three  cases  that  need  to  be  considered: 

(1)  C,  is  an  edge  of  P.  Then  it  is  extremal  by  definition.  It  is  also  obvious  that  any  local 
perturbation  of  C,  that  maintains  Q  C  P,  would  result  in  a  smaller  Q  thus  decease  the 
perimeter  of  Q . 

(2)  C,  has  a  single-pivot  v,.  C,  is  embedded  in  a  butterfly  which  is  centered  al  Vj  and  is 
supported  by  the  two  neighboring  chords  C,.,  and  C,.,.  By  Butterfly  Lemma  II,  we 
know  the  perimeter  determined  by  a  chord  embedded  in  a  butterfly  is  at  a  minimum 
when  the  chord  is  balanced.  And  we  also  know  this  minimum  is  the  only  local 
extremum,  thus  we  can  always  perturb  a  chord  and  increase  the  perimeter  of  Q 
determined  by  the  chord  unless  the  chord  is  already  extremal. 

(3)  C,  =  ab  has  a  double-pivot  {p ,  q).  Then  by  definition,  C,  is  extremal.  Furthermore, 
by  an  argument  similar  to  (2),  we  can  show  that  R(C ,,  p)  <  1  and  >5(C,,  q)  >  1;  C,  is 
balanced.    D 

This  lemma  provides  us  with  a  finiteness  criterion  for  the  problem  in  question  in  the 
sense   that   we   can   guess   the   sequence   of   extremal  chords  that   determines   the   maximum 
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perimeter  potato  and  then  check  the  balancing  condition  of  double-pivot  chords.  Since  the 
possible  combinations  of  the  extremal  chords  are  finite,  the  number  of  guesses  is  finite  and 
we  have  a  finite  algorithm.  It  is  clear  that  this  naive  algorithm  is  exponential  in  time.  Again 
the  way  to  cut  the  complexity'  down  is  by  using  dynamic  programming  to  organize  guesses, 
subproblcms,  and  partial  solution.  Next  we  will  show  an  0(n*)  time  algorithm  for  the 
maximum  perimeter  potato-peeling  problem  using  the  dynamic  programming  approach. 

Notice  that  in  the  discussion  above,  we  did  not  exploit  the  balancing  condition  which  is 
more  complicated  and  more  difficult  to  compute.  We  do  not  know  of  any  way  of  using  the 
balancing  condition  in  a  faster  algorithm  at  this  point. 

7.2.    An  Solution  to  the  Maximum  Perimeter  Potato-peeling  Problem 

Wc  need  the  following  notation  and  definition  for  the  presentation  of  the  algorithm: 

Notation.  For  two  points  p ,  q  ox\  the  boundary  of  P,  let  D[p ,  q]  denote  the  length  of  the 
path  counterclockwise  from  p  to  q  along  the  boundary  of  P,  if  there  is  no  reflex  corner  along 
the  path.  Otherwise,  D[p ,  q]  =  -».  For  two  points  p,q  on  the  line  determined  by  a 
extremal  chord  C,  let  d{pq)  denote  the  length  of  pq  if  the  direction  of  C  is  from  p  to  ^. 
Otherwise,  d{pq)  denotes  the  negation  of  the  length  of  pq. 

Definition.  For  C,C  €  E,  define  the  size  of  a  chain  (C,  C  )-chain  C  with  nodes 
i'O'  'i ^m)  as  follows: 


X  =   I  m 

where  c  and  c  arc  the  respective  reference  points  of  C  and  C  .  When  we  patch  two  chains  to 
form  another  chain,  this  definition  has  the  advantage  of  only  requiring  an  addition  for 
computing  the  size  of  the  combined  chain.  And  when  two  (C,  C  )-  and  (C  ,  r)-chains,  C  and 
C  ,  combine  to  form  a  convex  polygon  Q,  the  perimeter  of  Q  is  exactly 

\c\-\c\. 

Let  M{C ,  C  )  denote  the  maximum  length  of  all  chains  supported  by  (C ,  C  ). 

7.2.1.    An  Algorithm 

In  the  following,  we  describe  an  algorithm  for  finding  the  maximum  perimeter  potato. 
The  algorithm  proceeds  in  three  main  steps.  First,  the  set  of  extremal  chords  E  is  computed. 
Secondly,  a  dynamic  programming  process  is  performed  to  find  the  partial  solution  to  the 
subproblcm  of  finding  the  best  chain  supported  by  a  pair  of  compatible  extremal  chords. 
And  finally,  a  maximization  process  over  possible  combinations  of  these  partial  results  is 
done  to  find  the  solution.  As  m  the  previous  chapter,  instead  of  computing  the  potato  itself, 
wc  will  only  compute  its  size.    The  algorithm  goes  as  follows: 
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Algorithm  7.1.  Potato-Peeling  -  Perimeter  Measure 

Step  J  First,  wc  compute  the  set  of  all  extremal  chords.  With  each  extremal  chord,  we 
will  keep  track  of  its  two  endpoints  and  the  edges  of  P  on  which  the  cndpoints  arc 
located. 

Step  2       M  can  be  computed  in  stages.  In  stage  s  {s  =   1,  2 n-l),  we  compute  the 

entries  Af(C,  C  )  where 

C  €  E,,     C    €  Ej,      and    j-i  =  s(mod  n). 

(1)  In    stage    5    =     1,    we    have  j  =  /+l(inodn).     For    C  €  H,    and   C   i  H^, 
M{C,  C)  is 

rf(cT,)  +  D[v„  vj  +  rf(c"7^), 

where  c  and  c  are  the  reference  points  of  C  and  C   respectively. 

(2)  For  stage  j  >  1,  M{C,  C  )  is  determined  by  a  chain  that  has  an  extremal 
chord.    Thus  wc  can  compute  M(C,  C  )  as  follows:' 

M{C,  C)  =  max    {  A/(C,  C')  +  A/(C',  C)  }, 
C 

where  C    ranges  over  E^,  ik  has  the  range  i  <  k  <  j{moi  n). 

Step  3       The  maximum  perimeter  is  given  by 

max     {  M(C,  C  )  +  A/(C  ,  C)  }.  n 
c.c'  i  =. 

1 .1.1.    Analysis  of  Complexity 

From  the  previous  chapter,  wc  know  that  Step  1  takes  0{n^)  time  and  space.  In  Step  2, 
we  spend  O(n-)  time  to  carry  out  the  maximization  for  each  entry  of  M.  Thus  the  whole 
matrix  M  takes  O(n^)  steps  in  total  to  compute.  Each  entry  of  M  requires  0(n)  storage.  The 
cost  of  computing  matrix  D  is  dominated  by  this.  In  total,  O(n-)  space  is  needed  for  the 
storage  of  matrix  M.  And  finally,  Step  3  takes  0{n~)  time  and  does  not  require  more  than 
0{n)  space  to  keep  track  of  the  best  solution.    Putting  all  the  costs  together,  wc  have 

Theorem  7.2.    The  potato-peeling  problem  can  be  solved  in  0{n^)  time  using  0{n-)  space. 


'  Note  that  we  do  not  check  compatibility  between  Af(C,  C  )  and  M (C  ,  C  ).  If  they  are  not 
compatible,  then  they  cannot  form  an  optimal  chain.  Thus  it  is  not  necessary  to  exclude  them 
from  the  maximization. 
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7.3.    Unsolved  Problems 

Obviously,  in  dimension  higher  than  two-dimension,  the  problem  is  even  more  difficult. 
Consider  the  three-dimensional  case: 

Problem  7.1.  Given  a  non-convex  polytopc  P ,  find  a  convex  polytope  Q  contained 
in  P  with  maximal  surface. 

Wc  do  not  even  know  what  the  balancing  property  is  in  higher  dimensions. 

As  with  the  area  measure,  in  the  formulation  of  the  problem,  we  stated  that  we  want  to 
find  the  convex  included  polygon  with  maximum  perimeter,  no  matter  how  many  sides  it  has. 
An  obvious  variation  of  the  potato-peeling  problem  would  be  to  find  the  maximum  perimeter 
triangular  (or  quadrilateral)  potato.    In  general. 

Problem  7.2.  Given  a  n-sidcd  nonconvex  polygon  P  in  the  plane  and  a  fixed 
number  i,  3  ^  t  s  n,  find  a  convex  Jk-gon  Q  included  in  P  with  the  maximum 
perimeter. 


CHAPTER  8 


Miscellaneous  Problems 


In   this  chapter,  wc  will   go  over   three   problems  that  arc  variants  of  the  enclosure 
problems  wc  have  studied  in  previous  chapters. 

8.1.    The  Non-coDTex  Enclosure  Problem 

If  in  the  enclosure  problem,  we  relax  the  condition  that  the  enclosing  polygon  has  to  be 
convex,  we  have  the  following  problem: 

Problem  8.1.  Given  a  simple  n-sided  nonconvex  polygon  P  and  a  number  k, 
3  <  k  <  n,  find  a  A-sided  (non-convex)  polygon  Q  enclosing  it  having  minimum 
area. 

First,  we  give  some  notation  and  definitions  needed  for  stating  a  finiteness  criterion: 
Recall  that  in  Chapter  3,  A-  and  V-pairs  are  defined  with  respect  to  a  convex  polygon  Q. 
Now,  we  generalize  the  notion  to  simple  polygons: 

Notation.  For  a  directed  line  segment  s,  let  7  denote  the  same  line  segment  in  the  opposite 
direction  to  s. 

Definition.  For  two  edges  s,  and  Sj  of  a  simple  polygon  Q,  let  s,  be  7,  if  w,.,  is  a  reflex 
corner  and  j,  otherwise.  Similarly,  let  s,  be  J,  if  w^  is  a  reflex  corner  and  s  otherwise.  Then 
J,  and  Sj  are  an  A-pair  if  s,  and  s^  are  an  A-piir  according  the  definition  in  Chapter  3. 
Otherwise  j,  and  Sj  are  a  V-pair.  An  edge  is  an  A -edge  or  a  V-edge  depending  on  whether 
the  supporting  pair  of  lines  are  an  A -pair  or  not. 

Lemma  8.1.  All  edges  in  Q  are  balanced.    Furthermore,  every  A -edge  in  Q  is  flush. 
Proof.    By  an  argument  similar  to  the  one  for  Butterfly  Lemma  I.  a 
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Lemma  8.2.   For  a   V-pair  (e^,  e),   the  optimal  k-sided  chain  is  either  a  simple  balanced 
chain  or  it  contains  an  flush  edge. 

Proof,    follows  from  Lemma  8.1.  D 

Lemma  8.3.  There  are  at  least  two  extremal  edges  in  Q. 

Proof.    We  can  at  least  find  two  different  A -pairs,  thus  two  flush  edges,  \wQ.U 

Notation.    Let  M J^i,  j)  denote  an  (i,  y)-chain  of  length  m  with  minimum  area. 

M  can  be  computed  as  follows:    If  («,,  e^)  is  an  A -pair  and  path  (v,,  v  )  is  convex,  then 
^mi''  J)  ~  "«('.  j)^    Otherwise,  if  (e,,  e^)  is  a  V-pair  and  there  is  a  simple  balanced  {e     e )- 

chain,  C  of  length  m  enclosing  {  v,_, v^_|  },  then  M^{i,j)  is  the  area  of  p(C)  minus  the 

area  of  Acr  o  where  c  and  c  arc  the  reference  points  of  e^  and  e  respectively  and  o  is  the 
intersection  of  e,  and  ej.'   If  none  of  above  is  true,  then  set  A/„(i,  j)  to 

min      M„{i,i)  +  M„-„-i(i,j). 
/-I- 1 j-\ 

""O.     .     .m-  1 

M  can  be  computed  in  0{lrn^)  time,  given  matrix  a*  and  the  extra  areas  of  admissible  chains, 
a*  matrix  takes  0(n^\og  niog  k)  time  to  compute.  The  areas  of  admissible  chains  below  a  set 
of  points  can  be  computed  the  same  way  as  for  the  potato  peeling  problem  with  a  slight 
modification  in  0{n')  time. 

Having  computed  M,„{i,  j).  we  can  now  easily  find  the  area  of  Q  in  0{kn-)  time  as 
follows: 

min      «,('.  ;)  +  A/i-:-„0'.'). 

1=0  .    k-2 

The  overall  complexity  of  the  algorithm  is  dominated  by  the  O(n^)  time  needed  for  the 
computation  of  balanced  chains  which  is  also  the  bottleneck  of  the  potato  peeling  algorithm. 
It  remains  an  open  problem  whether  we  can  compute  balanced  chains  more  efficiently. 

Previously,  most  problems  in  computational  geometry  assume  some  kind  of  convexity. 
The  problem  in  question  is  unusual  because  there  is  no  assumption  of  convexity  neither  on 
the  input  nor  on  the  output  The  solution  we  have  for  the  problem  is  a  direct  derivation  of 
the  two  problems  we  have  solved.  However,  it  still  represents  a  significant  result,  since  it 
does  not  rely  on  or  take;,  advantage  of  the  assumption  of  convexity. 


'     Before    computing    q'(i,  _/),    we     have    to    take    the    convex    hull    of    the    part    of    P 
counterclockwise  bcrwccn  v    and  v,. 

"  Notice  the  reversal  of  supporting  lines  and  the  area  taken  for  M„(i,  j)  is  not  the  area  of  C ■ 
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8.2.    The  Shape  Enclosure  Problem 

In  this  section,  wc  look  at  a  restricted  version  of  the  enclosure  problcna: 

Problem  8.2.  Given  a  n-sided  convex  polygon  P  and  a  convex  *-gon  5,  3  ^  *  <  n, 
find  a  minimum  area  *-gon  Q  such  that  Q  is  obtained  from  S  by  scaling,  translation 
and  rotation  and  Q  contains  P . 

Definition.    Given  a  convex  polygon  Q  =  {q^ Qn-\)  ^^^  ^  point  p  in  the  plane,  let /f  be 

a  polygon  congruent  to  Q.  Let  r,,  0  s  i  <  n,  be  the  corner  of  R  corresponding  to  corner  q^ 
oi  Q.  If  q^,  r,,  and  p  are  collinear  for  0  s  /  <  n  then  wc  say  i?  is  a  transformation  of  Q 
centered  zip.  We  call  p  the  center  of  the  transformation.  Furthermore  if  the  transformation 
results  in  a  smaller  polygon,  wc  call  it  a  centered  shrinkage.  All  the  polygon  obtained  by 
transformation  of  P  centered  at  a  fixed  point  p  form  the  family  of  P  centered  at  p . 

Lemma  8.4.  Q  intersects  at  least  4  corners  of  P . 

Proof.  Since  we  are  interested  in  finding  minimum  area  enclosing  polygon,  we  can  assume 
that  edge  e  ol  Q  touches  at  least  one  corner,  say  v,  of  P.  If  there  were  no  other  points  of 
contact,  wc  can  shrink  Q  centering  at  v,  until  an  edge  e  oi  Q  touches  another  corner  v^  of  f. 
And  then,  if  there  is  no  other  intersection  point  between  P  and  Q  other  than  v,  and  v^,  we  can 
do  the  following:  (i)  If  e  is  parallel  to  e  then  we  translate  Q  in  the  direction  of  e,  until  there 
is  another  contact  between  P  and  Q.  (ii)  Otherwise,  we  shrink  Q  centering  at  a  point  p, 
p  -  eC\e  %o  that  there  is  another  contact  between  P  and  Q.  Hence  Q  touches  at  least  three 
corners  of  P .  Extend  the  three  edges  of  Q  that  touch  P  to  form  a  triangle  T .  We  say  that 
the  polygon  Q  is  embedded  in  T.  Note  that  the  polygon  Q  with  the  same  3  contacts  is 
determined  by  the  triangle  in  which  it  is  embedded.  If  Q  touches  4  corners  of  P,  then  we  are 
done.  Otherwise,  since  there  are  only  3  points  of  contact  the  corresponding  triangle  T  is  free 
to  rotate  around  P  while  maintaining  the  same  3  contacts  with  P  and  the  same  shape.  It  is 
easy  to  see,  from  Lemma  4.2,  that  we  can  always  perturb  7"  to  decrease  its  area  (thus  the 
area  of  Q),  contradicting  the  optimality  of  Q.   The  contradiction  completes  the  proof.    D 

Now  let  us  consider  the  case  where  5  is  a  quadrilateral.  It  is  clear  that  given  four  points 
(corners  of  contact  of  P  with  Q),  there  is  at  most  one  quadrilateral,  each  edge  of  which 
passes  through  one  of  the  four  points.^  The  amount  of  work  required  to  find  such  a 
quadrilateral  is  independent  of  the  size  n  of  the  problems;  it  takes  0(1)  time.  Thus,  the 
finitcncss  criterion  immediately  gives  us  a  naive  algorithm.  By  "guessing"  the  four  corners  of 
contact  and  computing  the  quadrilateral  for  each  guess,  we  could  have  an  0{n')  time 
algorithm.    But,  since  all  quadrilaterals  enclosing  P  "intersperse"  with  one  another,  there  are 


'  There  are  four  degrees  of  freedom  for  the  solution.    But  it  needs  to  satisfy  four  constrants 
to  be  of  the  shape  of  S. 
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linear  number  of  choices  for  the  four  corners  of  contact  instead  of  0{n^). 

Theorem  8.5.    The  minimum  fixed-shape  quadrilaterals  can  be  found  in  0{n)  time. 

Next,  we  look  at  the  general  case  where  *  a  4.  Essentially,  we  guess  the  four  edges  of 
Q  that  touch  corners  of  P .  The  four  edges  of  Q  intersecting  P  extend  to  form  a  quadrilateral 
Q  which  has  the  shape  of  quadrilateral  S  formed  by  extending  4  edges  of  5.  Wc  say  Q  is 
the  extension  of  Q  from  5  to  S.  So  the  problem  reduces  to  finding  a  quadrilateral  Q  with 
the  shape  of  S  enclosing  P  and  verifying  that  P  is  enclosed  in  Q  which  is  the  extension  of  Q 
from  5  to  5. 

Lemma  8.6.  Let  Q    be  a  k-gon  which  is  an  extension  of  a  quadrilateral  Q  from  S   to  S  and 
Q    t   closes  a  ngon  P .    Then  it  takes  0{n)  time  to  verify  that  Q  encloses  P . 

Proof.  Let  Q  consists  of  four  edges  a  ,  b  ,c  ,  and  d  which  correspond  to  the  edges  a,  b,  c, 
and  d  of  S  respectively.  Take  any  two  edges  of  Q  which  are  not  parallel,  say  a  and  b  .  Let 
a  and  b  intersect  at  a  point,  say  p .  Consider  the  family  Tof  Q  centered  at  p  .  Consider  an  s 
in  S,  other  than  a  and  b,  corresponding  to  an  edge  e  oi  Q .  The  tangent  to  P  with  the 
direction  of  e,  together  with  a  and  b  ,  determine  a  polygon  in  7.  Wc  denote  the  area  of  such 
a  polygon  by  A(s).    Then,  it  is  clear  that  Q  encloses  P  if  and  only  if 

A{c)  =  A{d)  =  mm  {  A{s)  |  5  is  an  edge  of  S ,  s  ^  c,  s  *  d  ]. 

As  we  advance  the  edge  s  of  5,  the  corresponding  tangent  to  P  also  advances.  Thus,  all 
relevant  tangents  can  be  determined  in  time  proportion  to  n,  the  number  of  sides  of  P .  And 
each  area  can  be  computed  in  0(1)  time  with  0(*')  precomputation.'*  So,  it  takes  only  0(n) 
(  =  0{n  -■  *))  time  to  verify  the  containment  of  P  m  Q .    D 

There  are  0(k^^  different  possibilities  for  S    and  each  one  takes  0{n)  time  to  process, 
thus 

Theorem  8.7.    Given  a  n-sided  polygon  P .  and  k  gon  S.  the  minimum  area  (perimeter)  k- 
gons  that  enclose  P  and  are  congruent  to  Q .  can  be  found  m  0{nk')  lime. 

8.3.    Triangular  Enclosure  with  Angular  Restrictions 

For   the   enclosure   problem,    if   we    insist   that   the    solution    be    a   triangle   and   put  a 
restriction  on  its  angles,  then  we  would  have  the  following  problem: 


^  All  that  wc  need  is  the  ratio  of  the  area  of  S  to  the  area  of  the  triangle  formed  by  extending 
three  edges  of  S . 
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Problem  8.3.  Given  an  n-sided  convex  polygon  P  and  a  number  a,  0  <  a  <  it, 
find  a  minimum  area  triangle  Q  such  that  Q  encloses  P  and  one  of  the  interior 
angles  is  equal  to  a. 

We  call  such  a  triangle,  an  a-triangle,  and  call  the  two  edges  next  to  the  a  angle, 
adjacent  edges,  and  the  edge  opposite  to  a,  the  opposite  edge.  Without  proof,  we  state  the 
properties  that  minimum  area  a-trianglcs  share  with  minimum  area  enclosing  triangles: 

Lemma  8.8.  One  of  the  edges  in  the  smallest  a-triangle  is  flush  with  the  given  polygon  P 

Lemma  8.9.  Minimum  area  a-triangles  intersperse  with  one  another. 

In  the  following,  we  give  the  sketch  of  an  0{n  log  n)  time  algorithm  for  finding  the 
minimum  area  enclosing  a-triangle.  Essentially,  we  divide  the  solution  triangles  into  two 
types:  (i)  triangles  whose  opposite  edges  are  flush  with  P,  and  (ii)  triangles  with  one  of  their 
adjacent  edges  flush  with  P. 

(i)  Assume  the  flush  edge  is  adjacent  to  o.  Exploiting  the  interspersing  property,  we 
are  able  to  preprocess  P  in  linear  time,  so  that  it  only  takes  0(1)  to  find  the  other  adjacent 
edge  next  to  the  flush  edge  and  the  opposite  edge.  With  the  precomputation,  we  can  find 
minimum  area  a-triangles  with  an  adjacent  edge  flush  with  P  in  linear  time. 

(ii)  It  is  relatively  more  difficult  to  find  minimum  area  a-triangles.  For  each  e  in  P,  it 
takes  0(n)  time  to  find  minimum  area  a-triangle  whose  edge  opposite  the  a  angle  is  flush 
with  e,  since  we  have  to  consider  all  0{n)  possible  choices  for  the  other  two  points  of 
contact.  Using  the  interspersing  property,  we  can  repeatedly  halve  the  number  of  choices. 
Thus  a  divide-and-conquer  algorithm  can  be  devised  to  run  in  0{n  log  n)  time  to  go  through 
all  possibilities. 

Theorem  8.10.    The  smallest  a-triangle  can  be  found  in  0(n  log  n)  time. 

8.4.    Unsolved  Problems 

We  do  not  have  finiteness  criteria  or  algorithms  for  problems  with  arbitrary 
combination  of  angular  and  lateral  constraints.    For  example,  we  do  not  know  how  to  solve^ 

Problem  8.4.     Given    a    simple    polvgon    P,    find    the    smallest   hexagon    in   W:.  f, 
enclosing  P . 

Recall  that  W^a  „  is  a  restricted  iamily  of  polygons  consisting  of  hexagons  for  which 
there  exist  two  edges  e,  and  f,,,,^,,   f,   tha;  arc  parallel  and  equal  in  length. 


See  [Dori-Ben-Bassat  (1984)]  for  an  approximation  algorithm 
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The    inclusion   version   of   the    shape    problem    is    unsolved;  we   do   not   how   find   the 
maximum  subset  of  a  polygon  with  a  certain  shape: 

Problem  8.5.    Given  a  polygon  P  and  S,  find  the  largest  polygon  included  in  P  and 
is  congruent  to  S . 


CHAPTER  9 


Conclusions 


9.1.    Remarks 

In  this  paper,  we  give  the  first  polynomial-time  solution  to  the  potato-peeling  problem 
(both  the  area  and  perimeter  metrics);  a  problem  with  only  exponential-time  solutions  known 
previously.  We  also  solve  the  problem  of  finding  the  minimum  area  enclosing  i-gon. 
However,  only  the  special  case  Jk  =  3  of  the  perimeter  version  of  the  same  problem  is  solved; 
it  is  not  clear  how  our  solution  can  be  generalized  to  cover  the  general  problem.  We  also 
cover  several  restricted  versions  of  the  enclosure  problem.  The  restricted  versions  of  the 
potato-peeling  problem  are  still  open. 

In  solving  these  problems,  we  introduce  such  geometric  concepts  as  balanced  chains  and 
interspersing  between  chains  and  polygons  which  are  interesting  in  their  own  right  and  hold 
the  key  to  the  solution  of  the  problems.  Computing  these  chains  is  the  bottleneck  to  a  faster 
algorithm  for  some  of  the  problems  we  have  discussed.  Whether  these  chains  can  be 
computed  faster  is  still  an  open  problem. 

It  is  quite  likely  that  some  of  our  algorithms  are  not  optimal,  although  we  do  not  know 
the  lower  bounds  of  these  problems.  It  remains  an  open  problem  how  these  lower  bounds 
can  be  established. 

9.2.    Further  Work 

This  thesis  touches  on  only  the  simplest  cases  of  optimization  problems.  Some  solutions 
may  not  seem  practical,  for  we  are  just  beginning  to  study  these  kind  of  problems.  We 
consider  what  we  have  done  a  good  starting  point  for  further  investigation  of  different  kinds 
of  geometric  optimization  problems.  There  are  a  lot  of  interesting  and  practical  problems  yet 
to  be  solved  We  feel  that  work  on  optimization  problems  should  be  pursued  in  the 
following  dirccoons: 

(1)  Metrics  other  than   just  area  or   perimeter;  combinations  of  area   and   perimeter  (for 
example,  the  area  over  the  square  of  perimeter). 

(2)  Restriction  and/or  metric  independent  algorithms. 
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(3)  Higher  dimensions. 

(4)  Parallel  algorithms  (sec  [Aggarwal-Chazelle-Guibas-O  "Dunlaing-Yap  (1986)].) 

(5)  Other  kinds  of  optimization   problems  than  the  inclusion  and  enclosure  problems;  an 
example  is  the  division  problem  which  can  be  defined  as  follows: 

Problem  9.1.    Given  a  convex  polygon  P,  find  a  line  dividing  P  into  two  polygons 
Q  and  R,  minimizing  (or  maximizing)  the  sum  of  a  certain  function  on  P  and  g.' 


'  For  example,  in  the  finite  element  analysis,  one  is  motivaied  to  subdivide  a  polygonal  shape 
into  more  or  less  "regular"  shapes.  That  seems  to  suggest  that  the  sum  of  the  area  over 
square  of  the  perimeter  of  rwo  resulting  polygons  be  minimized. 
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